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1 .0 INTRODUCTION AND SUMMARY 


This report considers possible implementations of the RF links for 
command, telemetry, and voice communication between the Shuttle Orbiter (SO) and the 
Space Station (SS). For these links, there are three implementable and realistic candidates. 
These are: 

1) All S-band link 

2) All Ku-band link 
and 

3) Hydrid S-band/Ku-band Link 

Another possibility would involve the use of UHF for voice communication, but this 
approach was not considered here because of potential frequency utilization problem. 

Table 1-1 shows link margin summary for Shuttle/Orbiter RF link 
implementations at either S-band or Ku-band.* The table indicates that the Ku-band 
implementation has significantly higher margins than the S-band implementation. This is 
due primarily to the use of high gain (24 dB) antennas on the Space Station end of the link. 

The advantage of larger antenna appertures provided by S-band is offset by 
the fact that only 0 dB antenna gain is assumed at the Space Station end for the S-band 
operation. This assumption is very conservative and it results rather low margins for the S- 
band link which carries the telemetry and voice from the Space Station to the Orbiter. If the 
Space Station antenna gain at S-band is increased to 3 dB, then the margins improve 
accordingly. At this point, we do not see any reasons why such antenna gain increase 
could not be acceptable if the decision is made to operate at S-band. 

The unique feature of the S-band implementation considered in this report is 
technique for transmitting two-way 16 kbps delta-modulation voice between the Shutde and 
the Space Station at the maximum range of 37 km. This approach, as proposed by 
Axiomatix, uses the existing standard subcarriers for communicating digital voice. For the 

’Hybrid use of S- and Ku-bands would involve same link margins, therefore, no listing is 
given for a "hybrid" approach. 




2 



Shuttle/Space Station link, the subcarriers in the range from 65 KHz to 95 KHz are 
considered as possible candidates. These subcarriers are typical of CIE and/or CIU 
equipment and are readily handled by the payload interrogator (PI). For the Space 
Station/Shuttle link, the standard 1.7 MHz subcarrier can be used for digital (16 kbps) 
voice. Consequently, by modulating the appropriate subcarriers at both ends of the S-band 
link a viable method for a two-way voice communication can be provided. 

From the S-band and Ku-band link budgets considered, it appears that 
either implementation alone or the combination of the two can meet the requirements for the 
Shuttle Orbiter/Space Station link. Furthermore, as described in this report, there appears 
to be a possibility of adapting some of the existing S-band payload link equipment for 
providing two-way voice capability to the link. This capability will permit the Orbiter to 
satisfy the voice link requirements without actually using the Ku-band MA system. 

From the standpoint of the Orbiter, this appears to be relatively low cost 
solution, requiring only minor baseband modifications to CIE or similar equipment and the 
corresponding transponder. Neither RF nor any antenna modifications will be required. 
However, remaining at S-band beyond the IOC capability will impose dual-band (S- and 
Ku-bands) operation requirement on the Space Station. The impact of the latter 
requirement remains to be addressed. 

The reference made above to the CIE is not to the actual CIE which is a part 
of DOD's SGLS network. What is meant here is simply an equipment which is "CIE-like" 
in design, i.e., it has more than one subcarrier on the link to the payload. In fact, it could 
possible be a modified PSP. Whatever equipment it may be, the main requirement is that 
the 16 kbps digital voice can be "patched in" into it at the Orbiter and also be recovered and 
separated from telemetry for on-board use. Commands, however, could originate either on 
the ground or on-board. 
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2.0 REQUIREMENTS 


Table 2-1 shows the requirements for the Orbiter/Space Station links as 
defined in Reference [1]. The requirements listed in this table can be interpreted to mean 
that (1) the commands and telemetry link can be satisfied by the present S-band Orbiter 
capability and (2) the two-way voice capability will be added when the Shuttle Orbiter EVA 
communication system is upgraded to be compatible with Space Station's multiple-access 
(MA) system. The table can be also interpreted that the S-band and Ku-band 
communication links can be used simultaneously to satisfy both the commands/telemetry 
and the voice link requirements. 

As explained in the subsequent sections of this report, however, there is a 
possibility to implement the addition of voice capability at S-band with minor modifications 
of the existing Orbiter's S-band equipment. Specifically, Axiomatix has come up with 
techniques for such an S-band implementation which will provide simultaneous 
transmission of commands and voice on the Orbiter/Space Station link and telemetry and 
voice on the Space Station/Orbiter link. 

The advantage of this approach is that it provides the voice capability for the 
Oribter/Space Station link during the early stages of Space Station operations. In principle, 
this S-band capability could be retained indefinitely without ultimately converting to Ku- 
band MA operation. But, the continuing utilization of S-band will necessitate additional 
equipment and antennas on the Space Station. Thus, the current thinking among Phase B 
contractors is to satisfy most of the space-to-space link requirements with the Ku-band MA 
system including the Orbiter/Space Station link. If such an approach is adopted, then the 
Orbiter/Space Station and the Space Station/Orbiter link can be implemented with the Ku- 
band MA link. This will eliminate the requirements for separate commands and telemetry 
link at S-band. 
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3.0 IMPLEMENTATION POSSIBILITY AND TRADE-OFFS 

3. 1 Shuttle Orbiter/Space Station Link 

For this link there are following three (3) possible realizations: 

(1) All S-band link 

(2) All K-band link 

and 

(3) Hybrid S-band/K-band link. 

The pros and cons of each of these implementations are discussed below. 

3.1.1 Shuttle Orbiter/Space Station Link at S-Band 

3 . 1 . 1 . 1 Possible Implementation at S-Band 

Figure 3-1 shows functionally an all S-band implementation of 
Shuttle/Space Station link. The important feature of the implementation shown is that the 
CIE (Communication Inferface Equipment) may be used as an interface between the 
baseband signals and the PI at the Orbiter end of the link. The CIE equipment was 
developed by Rockwell for use with DOD payloads, but it has certain features which make 
it suitable, with some modifications, for this link.* 

Specifically, the existence of more than one subcarrier on the forward and 
return links provide for transmission of both commands and voice on the forward link and 
reception of telemetry and voice on the return link. 

The return link implementation is based on the availability of two 
simultaneous subcarries, i.e., 1,024 MHz and 1.7 MHz which can be phase-shift-keyed 
(PSK) at the Space Station transponder (Xponder) with telemetry (16 kbps) and digital 
voice (16 kbps) , respectively . 

The payload interrogation (PI) can demodulate these subcarries and the data 
streams can subsequently be recovered by CEE. Figure 3-2 shows the functional diagram 

*Here we do not imply the actual CIE equipment as part of SGLS network. What we mean 
here is simply an equipment which is based on CIE design with modification. 
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Figure 3- 1 . Orbiter/Space Station Link at S-Band 
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of this implementation. It is important to note that the return link does not require 
considerable CIE modification because the CIE already has the features which permit such 
operation. 

The forward link, however, requires modification of the frequency-shift 
keying (FSK) scheme used for DOD payloads. Specifically, three tones ,i.e., 65 KHz, 
76 KHz and 95 KHz are used to transmit symbol rates up to 2000 symbols/second. The 
low symbol rate of the current system is not a limitation of PI modulation capability. This 
is because the PI modulator on the forward link responds primarily to the tone frequency 
and only in a second order fashion to the tone’s fine modulation structure. Consequently, 
several methods for transmitting both 2 kbps command and 16 kbps voice are available. 
Figure 3-3 shows the modification for the CIE which will provide this simultaneous 
capability. Part (a) of Figure 3-3 shows the frequency division multiplexing implementation 
(FDM) of command (CMD) and voice data for input to the payload interrogator (PI). As 
shown in the figure, the two tones, i.e., 65 KHz and 95 KHz, have been replaced by two 
subcarriers at 64 KHz and 96 KHz, respectively. These subcarriers are bi-phase 
modulated with the respective binary data streams 2 kbps (CMD) and 16 kbps (VOICE). 

The reason for changing each of the S/C frequencies by 1 KHz (95 KHz to 
96 KHz and 65 KHz to 94 KHz) is to place the CMD subcairier in the null of modulation 
spectra of the 94 KHz subcarrier. To further prevent the frequency interaction between the 
two S/C modulation streams bandpass/ filters (BPF) may be required for each modulated 
subcarrier prior to summing. The detailed requirements for these filters have to be 
determined in the future if this approach is adopted. 

Another method of implementing combined transmission of commands and 
voice is shown in part (b) of Figure 3-3. In this approach straight forward time-division 
multiplexing can be used. The maximum combinded data rate will be 18 kbps, plus some 
overhead to accommodate frame synchronization. This combined data stream is then bi- 
phase modulated on a 96 KHz subcarrier and delivered to PI. 


9 



4>-Mod 


CMD 
2 kbps 



S/C I 64 KHz 
<(>-Mod 


VOICE _ 
16 kbps - 



S/C I 96 KHz 



(a) Frequency Division Multiplexing Implementation of CMD and 
VOICE Data for Forward Link 



| CLK 


TOPI 

(Forward Link 
toSS) 


(b) Time Division MUX of CMD and VOICE on a Bi-Phase Modulated S/C 


CMD 





2 kbps 


MUX 

18 kbps + OH 

FSK 





Mod. 


VOICE 






16 kbps 





S/C (65 to 95 MHz) 


TOPI 

(Forward Link 
toSS) 


(c) Time Division MUX of CMD and VOICE on an FSK S/C 


Figure 3-3. Modification of CIE to Multiplex Command Voice on Forward Link 
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The advantage of the TDM is that it conserves baseband spectrum. On the 
other hand, it requires multiplexing/demultiplexing equipment. 

Part (c) of Figure 3-3 shows another implementation of TDM, this time with 
a frequency-shift keyed (FSK) modulator. The detailed trade-offs of PSK versus FSK 
modulation are not critical at this point in time. Both approaches, i.e., PSK and FSK, are 
discussed here for the purpose of demonstrating the alternatives available. 

Moderate modifications of CIE may be required to implement either one of 
the schemes shown in Figure 3-3. However, these modifications are only at baseband and 
should affect the RF portions of the PI. Neither will they affect the antenna system. 

Baseband modifications of the transponder* are also required in order to 
work with the modified Orbiter equipments. Again, these modifications are not considered 
severe, compared to modifications which would be required if RF modification were found 
necessary. 

3. 1 . 1 .2 Link Budget for SO/SS Link at S-Band 

To get a better quantitative idea of the performance of an S-band 
communication between the Orbiter and the Space Station baseline link budgets have been 
derived. These are shown in Tables 3-1 and 3-2. 

The Shuttle Orbiter to Space Station budget is shown in Table 3-1. For this 
budget it is assumed that PI/CIE equiment is used such as shown in Figure 3-2. It is also 
assumed that a linear phase modulation of PI carrier is implemented such as shown in (a) of 
Figure 3-3. 

In one of the versions of the link budget, approximately equal modulation 
indexes are assumed for the command subcarrier (Bj = 1.1 radian) and the voice subcarrier 
(B 2 =1.0 radian). With such a power division, however, the command channel carries 


*Here we imply either a modified network transponder or a transponder similar to those 
used by the payloads. 
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Table 3-1. S-Band Link, Shuttle Orbiter to Space Station 
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about 10 dB more margin than the voice channel, although both channels show adequate 
margins. 

To provide for more even distribution of link margins, the modulation 
indexes can be readjusted as shown in Table 3-1. With such power division link margins 
are about 8 dB for both channels. With such margins a total pointing loss can increase by 
as much as 6 dB with the links still being operational. 

Table 3-2 shows the Space Station to Shuttle Orbiter link budget. From the 
table, it is evident that this link has much less margin than the forward link. This is due to 
the fact the data rate is 16 kbps on both links. Because of the relatively low margins, a total 
increase of 3dB in pointing loss may "wipe out" the voice link. Thus, a budget for a 
contingency operation link with telemetry only going from Space Station to Orbiter is 
shown. 

The consideration of such a contingency depends on the performance of the 
S-band antennas at the Space Station. If such antennas have adequate coverage, then our 
concern for the voice link failing at 37 km is not valid. Thus, a better understanding of the 
antenna coverage of the S-band at the Space Station is in order before the final conclusion 
can be drawn. Another possibility is to consider the use of steered beam S-band antennas 
on the Shuttle. Such possibility is being considered as a possible upgrade of Shuttle 
avionics. 

3.1.2 Shuttle Obiter/Space Station Link at Ku-Band 

3 . 1 .2 . 1 . Implementation Consideration 

Figure 3-4 shows the functional representation of an all Ku-band link 
between the Orbiter and the Space Station. The implementation shown is based on 
utilization of the Ku-band multiple access (MA) system which is planned for the Space 
Station. Thus, at the Space Station end the Orbiter is handled as one of the multiple users 
of the RF resources assigned to the Station. Such handling, therefore, does not present 
additional burded to the Station. 
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On the Orbiter end, however, Ku-band equipment identical, or similar to 
that employed by Station’s user space vehicles (S/V) must be utilized. For the Orbiter this 
may involve considerable amount of new equipment installation and integration. Of 
particular, concern may be the integration of Ku-band antennas into the Orbiter's antenna 
systems ensemble. The exact nature of such an impact requires further, more detailed 
study of all of the factors involved. More information for such a study will become 
available when the configuration for the Space Station is better defined. 

3.1 .2.2 Link Budget for SO/SS Link at Ku-Band 

Table 3-3 shows a link budget for a Ku-band implementation of the 
command and voice communication links between Shuttle Orbiter and the Space Station. 
Table 3-4 shows the corresponding return link between the SS and SO. The power 
-division is assumed to be proportional to the data rates. 

In all, it appears that within the assumptions used in Tables 3-3 and 3-4, the 
Ku-band implementation of the SO/SS link has a relatively good margin. 

3.1.3 Hybrid S-Band/Ku-Band Implementation 

Figure 3-5 shows a functional representation of a hybrid link which utilizes 
both the S-band and Ku-band equipment. As indicated in the figure, the S-band link is 
used only for the commands and telemetry interchange between the Orbiter and the Space 
Station. In this case, the Space Station is treated as a payload. Voice communication takes 
place via a Ku-band link, such as the MA link used by the Space Station to communicate 
with other space vehicles. 

These appears to be no obvious advantage to this type of operation for the 
following reason: If the Ku-band MA system is installed into the Orbiter, then this system 
should be utilizable for all of the communication interchange between the Orbiter and Space 
Station. The only possible utilization of the S-band link may be in the earlier phases of the 
assembly of the Space Station. Once the Ku-band system is installed on the Orbiter then S- 
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Table 3-3. Shuttle Orbiter/Space Station Link at Ku-band for Command and Voice 
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Figure 3-5. Orbiter/Space Station Links at S-Band and Ku-Band 







band system use to the communicate with the Station will only impose the requirement of 
carrying S-band equipment on the latter beyond the IOC phase. 
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4.0 DISCUSSION OF TRADE-OFFS 


From the S-band and Ku-band link budgets presented above, it appears that 
either implementation alone can meet the requirements for the Shuttle Orbiter/Space Station 
link. Specifically, as described in Section 3. 1.1.1, there is a possibility of adapting some 
of the exsiting S-band payload link equipment for providing a two-way voice capability to 
the link. This capability will permit the Orbiter to satisfy the voice link requirements 
without actually using the Ku-band MA system. 

From the standpoint of the Orbiter, this appears to be a relatively low-cost 
solution, requiring only minor baseband modifications to CIE or similar equipment and the 
corresponding transponder. Neither RF nor any antenna modifications will be required. 
However, remaining at S-band beyond the IOC capability will impose dual-band (S and 
Ku-bands) operation requirement on the Space Station. The impact of the latter 
requirement remains to be assessed. 

Among the factors to be used in the evaluation process are the following 

criteria: 

• Authorized Frequency Band(s) 

• Sufficient Bandwidth 

• Minimum Interference 

• Development Risk 

• User Impact 

• Growth Potential 

• IOC cost. 

Explanation of these criteria is presented below: 

Authorized Frequency Band(s) 

This refers to exclusion of UHF bands from use in future NASA space 
programs. In the future this way change, but for the purpose of this report, UHF was not 
considered. 


Sufficient Bandwidth 


To accommodating simultaneous users. Space Station communication link 
require considerable amount of RF bandwidth. This is particularly true when MSC 
operates with multiple video channels. Bandwidth at S-band is limited. Ku-band is much 
more acceptable from this standpoint. 

Minimum Interference 

This implies both the effect of our links (i.e., space links) on other users as 
well as the effect of other emitters on our links. There is considerable history of ground- 
based interference on UHF and S-band links. In comparison, Ku-band looks promising 
from the standpoint of minimum interference. 

Development Risk 

This criterion is based on the technology state and the amount of available 
design history. From this standpoint, an all S-band implementation rates high. 

User Impact 

The amount of modifications required by Orbiter equipment can be 
considered as user impact. For S-band operation, this impact is minimal. For Ku-band, 
the impact is high, but the fact that MA Ku-band system will be developed from other users 
can reduce the cost of the impact. 

Growth Potential 

From this standpoint of the Space Station, the Ku-band system appears to 
be most desirable because it can accommodate more users in the future. Use of specialized 
S-band equipment on the Space Station limits growth. 

IOC Cost 

This criterion needs analytical modeling based on the experience and 
projections of Phase B contractors for the Space Station. From the point of view of the 
Shuttle Orbiter, the S-band approach appears to be most cost effective for the IOC. For the 
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long term considerations, this may not be true and the Ku-band approach may be more 
desirable. 


The criteria listed above indicate that the process of evaluating various 
approaches is rather complex. The current, on-going Space Station programs and design 
definitions will provide more data for carrying out a more accurate assessment 
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OVERVIEW AND SUMMARY 


The problem of radar interference arises when the shuttle rendezvous with the 
Space Station using its Ku-band radar which freqeuncy hops between five frequencies in 
the range from about 13.8 GHz to 14.0 GHz. 

Of particular concern is the radar frequency (out of five) at 13.779 GHz. This 
frequency is only 4 MHz away from the 13.775 GHz Space Ground Link (SGL) received 
by the Space Station from TDRS. The Ku-band radar signal enters Space Station’s SGL 
receiver via sidelobes of the high gain antenna pointed to TDRS. 

The estimates indicate that the peak Ku-band radar power intercepted by sidelobes 
of SGL receiver when shuttle is 10 km away may be in the order of - 40 dBm (peak) for a 
sidelobe level of - 30 dB. In contrast, the nominal value of SGL signal at the Space 
Station receiver may be about - 80 dBm. This presents interference-to-signal ratio of about 
40 dB. 

Several technical as well operational "work-arounds" to solve the interference 
problem have been considered by Rockwell Space Station study team. In principle, we 
agree with these findings within the scope of "pros and cons" defined by Rockwell team. 
In addition, however, we are considering such possibilities as: 

1) Operational fix 

Use only a single radar frequency for Space Station rendezvous taking advantage of 
station's large radar cross section which may be also enchanced by a suitable reflector. In 
this case the radar frequency will be in the middle of the radar band, i.e., around 13.88 
GHz, which is the frequency for beacon mode. The reflector must be optimized for 
returning circularly -polarzied (CP) RF pulses such as used in beacon mode. 

2) RF and Receiver Fixes 

Use of cancellation techniques such as: 
a) antenna sidelobe cancellation 
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b) precision pulse cancellation based on accurate amplitude and phase 
information with respect to the received radar pulse. 

The operational fix appears to be most desirable from standpoint of requiring 
minimum amount of new hardware for either Space Station or TDRSS. However, 
feasibility of precision approach and rendezvous using the beacon mode for skin return 
from Space Station remains to be established. 

RF and receiver fixes using cancellation techniques are still available if operational 
fixes are not accepted. The advantage of these technquies for solving SGL interference is 
that they do not require any operational coordination. In other words, these techniques 
belong into "set them and forget them" category. Moderate hardware development is 
required, however, at the Space Station receiver. 

We have also examined the Ku-band radar interference level to the reception by the 
Space Station MA receiver of the signals transmitted by users such FF, OMV, EMU and 
the NSTS itself. These signals arrive at SS in 13.64 - 13.70 GHz NB return band and in 
14.00 - 14.30 GHz WB return band. These signals may potentially be interfered with the 
extreme lower and upper radar transmission frequencies, respectively. 

Our analysis, however, indicates that with the exception of OMV and FF wideband 
return at maximum ranges, the radar interference is below thermal noise even when the 
users are at their respective maximum ranges. 

Link margin degradation to OMV and FF video WB return links is only about 1.5 
dB leaving margins of 2.6 dB and 4.3 dB, respectively. Furthermore, the estimated 1.5 
dB degradation is based on a worst case assumption that the Orbiter is approaching the 
Space Station along the axis of the high gain (46.3 dB) antennas used for communication 
with FF and OMV at long ranges. Operationally, this may not be the case. 

Thus, it appears that the major remaining problem is to deal with Ku radar 
interference to the reception of the SGL. 

The material that follows represents our analysis and other support data dealing 
with Ku radar interference to Space Station communications. 
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TABLE THAT FOLLOWS USES THESE GAINS FOR RADAR POWER CALCULATIONS 







Radar PRF = 3 KHz 

HGA = MA High Gain Antenna (46.3 dB) 

MGA = MA Medium Gain Antenna (24dB) 

HGA/SL = MA High Gain Antenna Sidelobe (- 30 dB) 





















































RELATIVE SPECTRAL DENSITY OF KU-RADAR 
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Plot of corresponding spectral envelope is shown on next page. 





SPECTRAL REPRESENTATION OF RADAR INTERFERENCE 
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AT ALL RANGES. 
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Means are available to do so manually on Shuttle. 
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= COHERENT AMPLITUDE DETECTOR 
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UP WEAK SIGNAL TO LEVEL REQUIRED BY DATA DEMODUALTOR. 
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CANCELLATION MAY PROVIDE FOR UNRESTRICTED APPROACH DIRECTIONS. 
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1 7) Margin Loss Due to Data - 1 .5 dB Radar interference adds to thermal noise 

18) Modified Margin 4.3 dB 
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* Rockwell International, Space Station Work Package 2, Definition and Preliminary Design Phase, 
Communication and Tracking Report, SSS-85191, Volume III, p. B-14, Table B-l. 1.2-1. 
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1.0 INTRODUCTION 


This interim report includes some preliminary considerations for antenna switching, 
power control, and AGC functions for the Space Station’s multiple-access (MA) system. 
The background for the general requirements for these functions is given below. 

Antenna Switching 

Antenna switching will be required to provide for optimum link conditions despite 
the relative movement of the MA system users with respect to the Space Station (SS). For 
movement within the proximity operations zone, several omni atnennas located at various 
points of the Space Station will provide for near- spherical coverage, but some criteria for 
selecting the proper antenna are required. Also, the quality of the received signal must be 
sampled to determine when to switch from an omni to a medium gain antenna and vice 
versa, the latter scenario being a part of power control. 

Power Control 

The users of the Space Station's MA system transmit signals at two widely 
separated rates. For example, the telemetry rate is about 100 Kbps and the video rate is 
about 22 Mbps. Also, the range to a user may be from about a few meters (EMU) to about 
37 km (OMV, NSTS). This creates a potential problem of interchannel interference and 
receiver overloading. Thus, means of controlling EIRP of users to reduce the Min/Max 
signal differential at Space Station’s receivers is required. 

AGC Function 

The multiple access system of the Space Station receives Ku-band signals from 
several users. These signals are picked up by antennas placed at various locations of the 
SS. The signals are amplified, filtered and supplied via a bus (IF or Fiber-Optical) to the 
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programmable signal processors (PSP) for data demodulation. The type and distribution of 
the AGC function has to be identified. 

The role of Axiomatix in dealing with these requirements can be summarized as 
follows: 

1) Identify key issue in each of these tasks, 

2) Analyze tradeoffs between available alternatives, 
and 

3) Recommend baseline approaches. 

As the title of this report indicates, the material presented here is only of a preliminary 
nature. The report addresses primarily the issues as we see them now. It also lists 
additional tasks to be carried out by us towards the recommendation of a baseline design 
for each of the functions specified. In some cases we also identify the approaches which 
appear promising at this point. In arriving at the material presented in this report we had to 
examine the technical material contained in the reports prepared by the Phase B contractors. 
Not all of the aspects of the MA system design are finalized in these reports. Thus, we had 
to assume certain baseline approaches which were rated in these reports as most promising 
candidates. In other cases, we had to apply our best engineering judgement and expertise 
to identify the issues. 

The material which follows is presented in a vu-graph form to simplify the process 
of communicating our ideas. As concepts and baselines will evolve and grow in 
complexity, a conventional format of the report will be utilized. 
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ANTENNA SWITCHING SCENARIOS | ^3 Avi«mot!v 

TASK REQUIREMENTS I MXIUll IcllIX 



switch from an omni to a medium gain antenna and vice versa, the latter scenario being a part of 
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preset level, the receiver is switched to a different 
antenna. (See next chart for comments regarding 
BER measurement). 
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Adv. BER detection response may be faster than that using sync words 

Disadv. Not all links are coded to allow for this feature. 
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Interfaces with other related functions at the Space Station 
(i.e., antenna selection, AGC) 






POWER | 
CONTROL COMMAND 










METHODS FOR USER EIRP CONTROL I AxiOITIcltiX 



(b) Amplifier Bypass and Antenna Switchin 










POWER CONTROL— REMAINING ISSUES I ^^§1 AxiOITIcltiX 
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Consider the interaction between the power control function and the related functions of 
antenna switching and AGC. 
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distribution of the AGC function has to be identified. 
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Complicates PSP design by requiring an AGC function within the PSP units. 






and antenna switching must be defined in quantitative terms. 
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EXECUTIVE SUMMARY 


In this report various technical aspects of implementing the power control for the 
Space Station Multiple Access (MA) system are addressed. Particular emphasis is placed 
on the RF/IF processing of the Ku-Band MA signals received by the station. A fiber optic 
(FO) cable link is assumed to be carrying the IF signals between the boom locations of the 
receiver/transmitter (R/T) units and the central processing unit located in the habitat module. 
Capabilities and limitations of using an FO cable link as a subunit of the power control 
function are examined. 

The motivation for the power control arises due to the following three factors: 
(1) the FDMA configuration of the return links, (2) the near/far problem, and (3) the 
dynamic range limitation of the Space Station's RF receiving and the IF distribution 
components. Figure 1 shows the power control function scenario. 

Two receiver architectures were examined for their ability to handle the FDMA 
signals. These architectures are: (1) a wideband receiver proposed by the McDonnel 
Douglas team and (2) a channelized receiver proposed by the Rockwell International team. 

The salient feature of the first approach is a single wideband downconversion of the 
entire return link band (300 MHz wide) to a wideband IF. Thus, each antenna of an 
antenna group is connected to a separate receiver which then connects to its own FO cable 
driver. The main advantage of this approach is the simplicity of the Ku-Band receivers. 
With a channelized approach a group of antenna receivers have a capability to amplify 
selectively the individual channels of the return link band. Although this approach is more 
complicated than the wideband approach it provides for a channel-dedicated AGC function 
within each receiver. 

Figure 2 presents a side by side comparison of our estimates for the third order 
intermodulation behavior of the two receiver architectures. The bandwidth assumed is that 
of a digital video link, i.e., 25 MHz. As indicated in the figure, the origins of the third 
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Figure 1 . Power Control Functional Scenario. 




Note : All power levels are referred to the system input. 


LNA 

Third Order 



Figure 2. Third Order Intermods for (a) Wideband Receiver and 
(b) Channelized Receiver. 
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order intercepts are different for the two receiver architectures. It can also be seen from this 
figure that a channelized receiver can, at least in principle, accommodate a 20 dB higher 
level of the input signals before the internal generation of the third-order intermods. This 
capability is due to the channelized AGC available with a channelized receiver. 

Because of the potential advantages of the channelized receiver, this configuration 
was adapted as a baseline for subsequent analysis to determine the power control window 
capabilities and limitations of the Space Station receiving system. Of particular significance 
was the analysis of the FO cable link capabilities. 

Based on a practical model of an FO cable link, it was established that the basic 
limitation for the power control window is the FO cable loss which causes the degradation 
in the system SNR at the optical receiver. In view of this degradation, it is estimated that a 
power window for the high SNR users, such as EMU and MSC-3,* may be limited to 
10 dB. For the lower SNR users, such as FF and OMV, the power control window of up 
to 15 dB appears feasible. Figures 3 and 4 show how the FO cable loss affects the system 
SNR transfer curve. 

However, a conservative recommendation is a 6 dB ± 3 dB power control window 
with a channelized receiver. The "flattening" of the curve at the higher input SNR means 
that it is more difficult to estimate with accuracy the large SNR values which are the result 
of the increasing signals due to the range closure. This is based on a baseline assumption 
that the power control signal is developed at the coherent demodulators located at the 
receiving end of the FO cable link. 

The possibility of using a PIN diode to control the power input to the user's 
transmitter, and thus control the output EIRP, was examined. It was determined that 
60 dB of power control can be obtained in this manner without antenna switching. 
Figure 5 shows a functional block diagram for this viable concept. 


* 


A three channel 8-PSK system. 
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CONDITIONS: 

N = 5 (five carriers, equal power) 
m = 0.5 (total mod index) 

P b = 0.00225 W (2.25 mW total average optical power at transmitter) 
B = 25 MHz (system bandwidth) 

RIN = -131 dB/Hz (transmitter noise) 


L = FO Cable Loss 
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Figure 3. System SNR Transfer Curve as Function of FO Cable Loss 
(EMU and MSC-3 Power Control Window of 10 dB). 
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N = 5 (five carriers, equal power) 
m = 0.5 (total mod index) 

P b = 0.00225 W (2.25 mW total average optical power at transmitter) 
B = 25 MHz (system bandwidth) 

RIN = -131 dB/Hz (transmitter noise) 
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Figure 4. System SNR Transfer Curve as Function of FO Cable Loss 
(FF and OMV Power Control Window of 15 dB). 
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Figure 5. Functional Block Diagram for User Power Control. 



1.0 INTRODUCTION 


The subject of this report is the return link power control for the user segment of the 
Space Station's Ku-band multiple access (MA) system. The motivation for the power 
control is due to the following factors: 

(1) the baseline configuration for the return link(s) is the frequency division 
multiplexing (FDMA), 

(2) there is a large spread between the signal levels received by the Space Station 
from the users at the maximum and the minimum ranges, 

and (3) the limitations of the dynamic range capabilities of the Space Station's RF 
receiving and the IF distribution components. 

Figure 1-1 illustrates a functional scenario for the Ku-band MA system's return link power 
control. As shown in the figure, a number of users may be transmitting to the Space 
Station simultaneously from different ranges. In some cases, they may be transmitting to a 
common antenna. In other cases, they may be transmitting to different antennas and yet 
their signals may be combined in either a common RF or a common IF equipment. 

Figure 1-2 shows the signal power received at the Space Station from various users 
as a function of user range and signal modulation. From the figure it is evident that a 
power differences of as much as 90 dB may exist. For example, a telemetry received 
signal from a free flyer (FF) at 2000 km may be at about -105 dBm while a video signals 
received from another FF at 1 km may be at about -14 dBm. 

It must be pointed out, however, that just because one signal is far above the other, 
it does not mean that there will be an interference to the weaker signal. A great deal 
depends on such factors as the return link frequency plan, adjacent channel filtering at the 
user transmitter and the receiver configuration at the Space Station. Of more direct concern 
is the near/far problem which may exist for the users within the same transmission band 
and serviced by the same antenna. 
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The near/far problem, if not alleviated by the user power control, may lead to a 
performance degradation of the Space Station’s receiving and signal processing equipment. 
This degradation may be in the form of: 

(1) intermodulation products generated in such receiver elements as the front end, 
the down converter and the cable driver. 

(2) interchannel interference in both the receiving and the signal processing 
equipment, 

and (3) signal limiting in various element of the receiving equipment, particularly the 
drivers which feed either the RF/IF or the optical signal distribution cables. 

The power control function can be implemented by sensing either the power levels of the 
received signals or their signal-to-noise ratio (SNR) and then using this information to send 
the power control commands to the users via the forward link. The required information 
can be obtained from the received signals on either an absolute basis, or a relative basis, or 
some combination of both. In this report we baseline a method for measuring the SNR of 
the received signals and using this information to implement the power control. 

The degree to which the power, or more generally the EIRP, is to be controlled 
depends on a tolerable "power window" at the Space Station. Different criteria influence 
the power window in a different way. The criteria for the interchannel interference may be 
different from the criterion dictated by a limitation in the dynamic range of an optical cable 
driver. For example, a 2 or 3 dB power window may satisfy most of the criteria and thus 
eliminate most of the performance degraders. On the other hand, imposing on the users a 
power control within such narrow limits will increase the complexity of user's transmitting 
equipment. 

In the material which follows various criteria for the power window are examined 
and various implementation tradeoffs are discussed. 
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2.0 TECHNICAL CONSIDERATIONS 

2.1 Received Signal Spread 

In this section we take a closer look at the significance of the signal level spread due 
to the near/far problem. For the purpose of discussion, we will consider primarily the 
video link signals and we will assume that a similar argument holds for the narrowband 
signals, i.e., the telemetry and the digital audio. 

Table 2.1-1 summarizes the signal level characteristics of the Ku-band video return 
links. Of particular note in this table are the minimum signal-to-noise ratios of the different 
users. It must be noted that the spread in SNRs ranges from 15.6 dB (FF) to 27.8 dB 
(MSC-3). This represents a spread of about 12 dB. The spread of the margins, however, 
extends only from 4.8 dB (MSC-3) to 9.6 dB (EMU), i.e., a difference of less than 5 dB. 
If the margins are reduced to zero, then there is a spread in SNR of 13.2 dB. 

The point that we are making here is that because of the nature of the links, it is not 
possible to equalize the working SNRs of these different signals. Neither is it possible to 
equalize the margins. This means that when a combination of signals from various users 
appear at the same antenna, there will be an unavoidable spread in signal levels which is not 
a function of range. Also, the criteria for power control may have to be different for 
different links because of this intrinsic difference. As a result, the power control algorithm 
may have to be tailored to each user. 

Figure 2.1-1 shows video link signal levels and the corresponding margins as the 
function of the range. 
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2.2 The Effect of Receiver Implementation on System Intermods 
2.2.1 Wideband Ku-Band Receiver 

Figure 2.2. 1 - 1 shows a wideband Ku-band receiver implementation. The salient 

feature of this implementation is a single wideband downconversion of the entire return link 
band (300 MHz wide) to a wideband IF. As shown, the output of the downconverter 
extends over a bandwidth of 450 to 750 MHz. Each antenna of an antenna group is 
connected to a separate receiver which in turn connects to its own FO cable driver. 

The main advantage of this approach is the simplicity of the Ku-band receivers. 
Each receiver consists of a wideband (300 Mhz) bandpass filter, followed by a low noise 
amplifier (LNA), which in turn is followed by a wideband downconverter. The function of 
the bandpass filter is to restrict the signals applied to the LNA to only those within the MA 
return link band and to provide the image rejection. The low noise amplifier provides the 
gain required to maintain a low noise figure despite the downconverter loss. The output of 
the downconverter is then increased by a fixed gain which is tailored to provide the 
required input level to the cable driver. 

At the other end of the FO cables, the signals are detected by FO receivers and 
applied to an IF switch matrix. This matrix permits any antenna to be connected to any 
PSP unit. Channel separation and bandpass filtering takes place in the PSP units. Thus, 
the intermodulation can occur at any point within the receiver/FO cable driver chain. 
Because there may be several signals received simultaneously by any one antenna at any 
time, it is important to maintain a high degree of linearity in this receiver in order to keep 
the IM interference at tolerable levels. 

Of specific concern is the effect of the third order intercept which is due to the FO 
cable driver assembly. This intercept is reflected to the system input at a level of about -36 
dBm. The effect is shown in Figure 2.2. 1-2. In this figure all of the levels have been 


This approach is being considered by the McDonnell Douglas team. 
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Figure 2.2. 1-2. Wideband Receiver Third Order Intermodulation 
as Function of Input Signal. 
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normalized to a bandwidth of 25 MHz which is typical of a digital video link. From the 
figure it is evident that as the input signals increase above -60 dBm, the intermodulation 
products begin to rise above the thermal noise in the 25 MHz bandwidth. 

Figure 2.2. 1-3 shows the channel degradation due to the third order 
intermodulation as the function of the input signal level. This figure shows that to keep the 
degradation below 0.5 dB, the received input signal level should not exceed -59 to -60 
dBm. This requirement places a rather stringent demand on the user power control. 


18 



Figure 2.2. 1-3 . Channel Degradation Due to Third Order Intermodulation in 
the Wideband Receiver. 
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2.2.2 Channelized Ku-Band Receiver Implementation 

Figure 2.2.2- 1 shows a group of antenna receivers which have a capability to 
amplify selectively the individual channels of the return link band. This capability is 
provided by implementing each antenna receiver with several parallel receivers. The 
number of parallel channel receivers in each antenna receiver is equal to the total number of 
channels in the return link band (or bands). 

As shown in the figure, each antenna is connected to a wideband bandpass filter 
followed by an LNA. The LNA is connected to an N-way power divider. The N output 
ports of the divider feed N channel receivers. Each channel receiver has its own 
downconverter (not shown) and at least one, or possible two, IF amplifiers.* ** Thus, the 
channel receivers can be tuned to the signals picked up by a particular antenna. Similarly, 
the return channels picked up by other antennas are applied to their respective channel 
receivers. The outputs of all M antenna receivers are then combined and applied to the 
cable driver. 

The advantage of this receiver implementation is that once the received signals are 
separated into the individual frequency-selective channels, the levels of these signals can be 
controlled by AGC prior to the summation at the input to the cable driver. Such control 
provides two following features: (1) the levels of different received signals can be 
"equalized" at the input to FO cable driver and (2) the signal levels can be held at a 
relatively constant level despite signal increase due to range closure. 

The availability of a channelized AGC prior to the FO cable driver implies that, at 
least in principle, the input signals can vary over a wide range of power levels, before the 
intermodulation becomes a problem. 


* 

This approach is being considered by the Rockwell International team. 

||c ]|c 

The detailed architecture of the channel receivers is not addressed here. These receivers may be either 
single conversion or double conversion type. 
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Figure 2.2.2- 1. Channelized AGC Equalizes Signal Levels into the Cable Driver. 




It is important, however, that the circuitry preceding the line driver is not 
overloaded by the received signals prior to the frequency selective channelization. This 
requirement is due to the fact that the LNA and the downconverters of the channel receivers 
are wideband and thus have a limited dynamic range compared to the channel receivers. 
Consequently, the power control is still required but its function is primarily to reduce the 
interchannel interference due to spectrum sidelobes and to prevent intermodulation in the 
wideband LNA and the downconverters. 

With the channelized receiver, it is reasonable to assume that the primary source of 
the intermodulation distortion is the LNA unit where all the channels are amplified on a 
broadband basis. Figure 2 . 22-2 shows the third order intermodulation level originating in 
the front end (LNA) portion of a channelized Ku-band receiver. Figure 2.2. 2-3 shows the 
degradation, on a channel basis as function of the input signal level. Comparing this figure 
to 2.2. 1-3, we can see that the channelized receiver can accommodate about 20 dB more 
input signal level than the wideband receiver. This means that, at least in principle, a wider 
power control window can be tolerated with a channelized receiver. This, of course, is 
based on the assumption that a channelized power level equalization takes place ahead of 
the FO line driver. 

The disadvantage of the channelized approach is the relative complexity of the 
antenna receivers. This disadvantage may be offset by the intrinsic component redundancy 
provided by the availability of several channel receivers within each antenna receiver. 
Thus, if one channel receiver fails, another may be used in its place. 

Although a noncoherent AGC is shown in Figure 2.2.2-1 for the purpose of 
functional clarity, such AGC may not necessarily be the final choice, and a coherent AGC 
may be called for this application. 
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Figure 2 . 22 - 2 . Channelized Receiver Front End Third Order 
Intermodulation as Function of Input Signal. 
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Figure 2.2.2-3. Channel Degradation Due to Third Order Intermodulation in Front End 
LNA as Function of Signal Input to a Channelized Receiver. 
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2.3 Cable Driver Considerations 

The dynamic range of signals applied to either an RF/IF or a FO cable driver is 
among the factors, which must be taken into account when specifying the degree of power 
control. The specific effects to consider in this case are the intermodulation products and a 
limited dynamic range of practical cable drivers. However, as explained in the preceding 
sections, the architecture of the Ku-band receiver ahead of the cable driver influences to a 
great degree the interaction between the capabilities of the driver and the power control 
requirement. For example, if channel separation is provided prior to the cable driver, then 
such means as AGC can be used to equalize the levels of different channels at the input to 
the driver. With the spread in signal levels reduced in this manner, the requirements for the 
power control can be relaxed somewhat without seriously impacting the performance of the 
cable driver unit. 

In this report we are baselining an FO cable driver such as General Optronics Corp. 
Model GO DIP 1000 or equivalent. Presented below are the factors which determine the 
tradeoffs between the drivers capability and the power control requirements. 
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2.3.1 FO Cable Driver Functional Diagram Description 

Figure 2.3.1- 1 shows a functional diagram of an FO cable driver. We constructed 
this diagram based on the information available to us from the data sheets on the GO DIP 
1000 laser module. As shown in the figure, an RF/IF amplifier supplies RF power to an 
RF to optical pow'er converter which in this case is an injection laser diode. The laser diode 
is biased into the conduction region by the bias current. In the absence of the RF 
modulation, the laser diode emits an optical power which is proportional to the bias current. 
When the RF power developed by the driver amplifier is applied to the laser diode, the bias 
current following through the diode is amplitude-modulated at the RF/IF rate. This results 
in an intensity modulation which is propagated along the FO cable and is detected at the 
other end by an appropriate photo-detector. The optical limitations of the FO cable driver 
shown arise from the following: 

(1) Third order intermodulation 

(2) Maximum power clipping 

and (3) Optical noise. 

In addition to these optical limitations, there is a limitation contributed by the 3rd order 
intermodulation generated in the RF driver amplifier. 

The optical limitations are functions of the bias current. If the bias current is 
increased, the allowable optical modulation is reduced, but the optical noise is also reduced, 
which is a desirable factor. Also, when the allowable modulation is decreased, the 
maximum clipping can be reached with a smaller RF signal, a result which is undesirable 
from the standpoint of obtaining a maximum dynamic range out of the optical power 
converter (i.e., laser diode). Thus, finding an optimum operating print for a laser diode 
becomes an important consideration. 

Furthermore, the intermodulation contributed by the driver amplifier is an important 
factor. In most cases it is this intermodulation which limits to the dynamic range of an FO 
driver. For example, if the third order intercept of the driver amplifier is at 24 dBm and the 
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Figure 2. 3. 1-1. FO Cable Driver Functional Block Diagram. 
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system gain, including the gain of the driver amplifier, is about 60 dB, this third order 
intercept is reflected to the system input as a -36 dBm intercept. Consequently, it becomes 
evident that such factors as laser diode bias current and the system gain are some of the 
parameters which have to be optimized before the actual capabilities of an FO line driver are 
established. 
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2.3.2 Optimizing System Gain Ahead of the FO Line Driver 

Figure 2.3.2- 1 shows the optical driver output as the function of the system gain 
preceding the optical driver. As is shown in the figure, there is a certain fixed level of the 
optical noise which is delivered by the optical driver to the optical cable. This level 
(-51 dBm) is based on the optical noise density of -125 dBm/Hz in a 25 MHz bandwidth. 
This noise is the relative intensity noise (RIN) and is a function of the laser diode. 

On the upper end, the maximum modulated signal out of the optical element is 
limited at about -2.7 dBm which corresponds to 50% modulation with respect to the 3mw 
peak unmodulated output of the laser diode. The maximum modulated signal is the 
function of the bias current, which in this case is assumed to be set to allow for 50% optical 
modulation. 

For low values of system gain, the desired signals are degraded by the optical noise 
(i.e., RIN) of the laser diode. As the system gain is increased, the received signals and the 
thermal noise are amplified until the thermal noise exceeds the RIN. At about 60 dB of 
system gain, the degradation due to the RIN is only 0.1 dB. From this point on any 
additional gain reduces the dynamic range of the optical driver. 

Thus, from Figure 2.3.2- 1, it is evident that the dynamic range of the optical driver 
can be traded off for some degradation in the link margin. This degradation would come 
from operating the system at some lower gain, say 50 dB, where the degradation due to the 
contribution from the RIN is about 1 dB. At system gain below 40 dB, the margin 
degradation increased inversely proportional to system gain as the received signals get 
buried in the RIN. This is shown in Figure 2.3.2-2. 

The SNR of the received signal must be also considered when determining the 
system gain. As shown in the figure, the SNR of 15.6 dB (FF) can be handled by a gain 
of about 76 dB while that of 27.8 dB (MSC-3) needs a reduced gain of about 64 dB to 
avoid clipping. But, as has been said earlier, there is a need to increase the system gain 
above the point at which the optical noise contribution is negligible. Thus, the discussions 
that follow will assume a baseline minimum value of system gain of about 60 dB. 
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Figure 2. 3. 2-1. Optical Driver Output as Function of System Gain 
Preceding the Driver. 
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Figure 2. 3. 2-2. Link Degradation Due to FO Driver Optical Noise (RIN) 
as Function of System Gain Preceding the Driver. 
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2.4 Power Control 

In this section we explore the various aspects of the power control. Specifically, 
we consider how the receiver implementation at the Space Station impacts the power 
control "window". 

2.4.1 Power Control for a Wideband Receiver 

In this section we will consider a baseline power control window for a wideband 
receiver such as shown in Figure 2.2. 1-1. We will consider the effect of a particular power 
control window on the operation of the FO cable driver as it is influenced by the wideband 
nature of the receiver preceding the driver. 

As shown in Figure 2.2. 1-1, each wideband receiver is dedicated to a single 
antenna and it feeds a single cable driver. Thus, only the signals picked up by a particular 
antenna are applied to a related FO cable driver. However, an antenna may be receiving 
signals from different users at different ranges. For example, signals from OMV and FF 
may be picked up by the same high gain (46.3 dB) antenna at their extreme ranges. Thus, 
they will appear together at the output of the same FO line driver. Also, signals from an 
MSC-3 and an EMU may be picked up by the same hemi antenna. This may happen when 
the reception from an EMU is switched from the medium gain (24.3 dB) antenna to the 
reception by a hemi (3 dB) to permit a close-in operation. If this switch-over occurs at a 
range of about 60 meters, the margin for the reception from the EMU will be about the 
same as the margin for the reception at 1 km. 

At this point let us baseline a 3 dB power control window and determine what effect 
this window has on the system operation. First, we must consider the third order 
intermodulation intercept referenced to the system input. As shown in Figure 2.2.1-2, this 
intercept is at -36 dBm, and the performance degradation resulting from this intercept is 
shown in Figure 2.2. 1-3. From that figure we conclude that to keep the performance 
degradation at less than 0.5 dB the input signals should not exceed about -59 dBm. 
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Figure 2.4. 1-1 shows the effect of a 3 dB power control on the video links received 
from the FF, OMV, EMU and MSC-3. As can be seen from this figure, the power control 
keeps all signals below the -59 dBm level with a margin of about 5 dB. Thus, there are no 
significant intermods due to the RF/IF circuitry of the wideband receiver. The signal level 
spread is 15.2 dB between the minimum (FF at 2000 km) and the maximum (MSC-3 at 
1 meter) signals. This corresponds to the existing link budget spread in SNRs of these 
different users, i.e., 12.2 dB plus the power control window of 3 dB 

Figure 2.4. 1-2 shows the four received signals as they may appear at the output of 
the FO cable driver. In this figure all of the signals pertaining to the output of the cable 
driver have been referred to the system input by subtracting 60.5 dB of the system gain.* 
Thus, the +4.8 dBm (3 mw) value of the peak modulated optical power appears as -55.7 
dBm in Figure 2.4. 1-2. Other output-referenced signals are also reduced accordingly. 
Both the thermal and the optical noise are based on a channel equivalent noise bandwidth of 
25 MHz (i.e., 74 dB-Hz). 

From Figure 2.4. 1-2, it is evident that even with a relatively narrow power control 
window of 3 dB, there is not much operating margin at the output of the optical cable 
driver. This is due to the fact that the SNRs of the four signals differ by about 12 dB in the 
25 MHz bandwidth. 

Of particular concern is the fact that the maximum signal of the MSC-3 is only 
one dB away from the level of the maximum signal modulation. This means that any drift 
in gain in the portion of the receiver which precedes the optical driver may drive the 
MSC-3 signal into limiting along with the other signals which may be present at the same 
antenna output. In fact, the RMS addition of the highest MSC-3 signal (-64.2 dBm) and 
the highest EMU signal (-67.2 dBm) will make the combined signal exceed the -63.2 dBm 


* 


At this system gain the degradation due to the optical noise does not exceed 0.1 dB 
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Figure 2.4.1 -1 . Min/Max Received Power Spread for a 3 dB Window Power Control 
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Figure 2.4. 1-2. Optical Cable Driver Output Levels for 3 dB Window 
Power Control (Wideband Receiver). 









threshold by about 0.5 dB. And this does not even include the secondary effect of adding 
other lower level signals. 

One of the solutions to this problem is to narrow the power control window to 2 dB 
or less. However, in view of the added complexity to the power control function 
necessitated by such a "tight" window this approach does not look attractive, particularly 
within the user segment. Instead, one may consider operating the system with a reduced 
margin by lowering the gain prior to the optical driver. The margin degradation in this case 
will then by due the increased level of RIN relative to the system thermal noise. Another 
alternative would be to lower the bias on the FO driver and thus allow more modulation 
range on the cable driver. 

Perhaps the most technically workable alternative, although a more expensive one, 
may be to increase the power handling capability of the laser diode used for the FO cable 
driver. The trade-offs in this area are, however, beyond the scope of this report. 
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2.4.2 Power Control with a Channelized Receiver and AGC 
2.4.2. 1 General Consideration 

As was explained in Section 2.2.2, the advantage of a channelized receiver is that 
the received signals can be separated (i.e., bandpass filtered) according to their channel 
frequencies prior to being applied to an FO cable driver. Such separation permits AGC to 
be applied to the individual channels to achieve some degree of signal level equalization at 
the input to the FO cable driver. In this section we consider various aspects of the 
interaction between the power control requirement and the channelized AGC function. For 
the purpose of comparison with the broadband receiver, it is initially assumed here that a 
3 dB power control window is implemented. The effect of using wider power windows is 
then considered in view of the AGC capability inherent with a channelized receiver 
approach. 

Before considering quantitatively a specific situation, the effect of the channelized 
AGC prior to the FO cable driver needs to be discussed. Figure 2.4.2. 1-1 shows, in a 
qualitative way, the effect of a channelized AGC on the signal and noise levels at the output 
of the FO cable driver. As shown in that figure, the dynamic range of the FO cable driver 
is limited at the upper end (2) by the maximum signal modulation (MSM) level and at the 
lower end by the optical noise (6). 

Consider first what happens to the signal (1) when the input signal increases due to 
a decrease in range. At the initial acquisition, the signal will be at an SNR determined by 
the thermal noise. Without the AGC the signal will increase due to a closer range as shown 
by the dotted line (1) until it reaches the maximum signal modulation level (2). At this 
point the clipping will occur in the optical transmitter and the signal will be distorted. The 
level of the thermal noise, however, will remain at a constant level (4) regardless of the 
increase in the signal level. 

When the AGC is activated, the signal will be kept at a nearly constant level (3) 
despite an increase in signal strength. At the same time the AGC action will, due to the 
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RELATIVE SIGNAL AND NOISE LEVELS— dB 



Figure 2.4.2. 1-1. Effect of Channelized AGC on Relative Signal and Noise 
Levels at the Output of FO Cable Driver. 
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reduction in gain, lower the level of the thermal noise reaching the FO cable driver. This is 
shown by (5). At some point, however, the front end thermal noise will be at the same 
level as the optical noise, and thus both noises will be contributing to the total noise at the 
output of the FO cable driver. As the signal strength continues to increase further, the 
AGC action will push the thermal noise below the level of the optical noise. Beyond this 
point the output signal-to-noise ratio will remain constant and it will be limited by the 
optical noise. 

Figure 2.4.2. 1-2 shows FO cable driver output SNR vs. system input thermal 
SNR. The curves presented there demonstrate the effect of the AGC upon pushing the 
thermal noise below the optical noise. In this case the optical output noise consists of both 
the random noise and of third-order intermods due to the optical driver per se, i.e., 
exclusive of the RF/IF driver amplifier circuitry. The initial ratio of the thermal to the 
combined optical noise is 14 dB for the assumed model of the cable driver. 

The curves show different behavior as the function of the AGC threshold above the 
thermal noise. The important common feature of the three curves shown is the leveling off 
of the output SNR as function of the increasing input SNR. Curve (1) is for the AGC 
threshold set at 26.7 dB above their thermal noise. This AGC setting keeps the maximum 
level of any signal at a level 6 dB below the peak optical modulation level of -2.7 dBm. 
Curve (2) corresponds to an AGC threshold of 19.3 dB which is the minimum SNR for the 
OMV signal. Curve (3) corresponds to a threshold setting of 10 dB and its purpose is to 
provide additional information on the behavior of the output SNR. 

From Figure 2.4.2. 1-2, it is evident that the output SNR ratios level off at their 
respective asymptotes which are equal to the respective AGC threshold setting plus 14 dB. 
The effective degradation in the output SNR w/r to the input SNR is also represented by a 
set of respective curves. Although not necessary detrimental to the link margins, the 
leveling of the output SNR curves presents a limitation on the accuracy with which the 
output SNR can be measured. Specifically, if the AGC threshold is located beyond the FO 
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driver, the leveling reduces the accuracy with which the power control decision can be 
made. We will address this subject in more detail elsewhere. 
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♦Operating Curve for data in Figure 2.4. 2.2-1 


Figure 2.4.2. 1-2. FO Cable Driver Output SNR vs. System Input Thermal Noise SNR. 
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24.2.2 Example of a 3 db Power Window 

Figure 2. 4.2. 2-1 illustrates the effect of a 3 dB power control window on the FO 
line driver output. This figure is similar to Figure 2.4. 1-2 except that all signal levels are 
referenced to the FO cable driver output instead of the system input. 

For the example shown in Figure 2.4.2.2-1, it is assumed that the AGC threshold 
is set at a level equal to the signal received from an OMV. In other words, the AGC action 
prevents all other signals from exceeding the nominal level of an OMV signal. As shown 
in the figure, this level is 12.5 dB below the maximum signal modulation level (-2.7 dBm) 
of the FO driver. The absolute level of this AGC threshold is -15.2 dBm. 

From Figure 2.4.2.2-1, the following can be observed: 


(1) Signals from an FF are below the AGC threshold level and, therefore, they are 
not affected by the AGC action. 

(2) Maximum signals from an OMV remain at the level of the AGC threshold. 
However, as the received OMV signal increases within a 3 dB power window 
the AGC reduces the channel gain accordingly. This pushes the thermal noise 
down by 3 dB. 

(3) Maximum signals from an EMU and the MSC-3 are affected by the AGC even 
outside of the power control range. Specifically, this means that the AGC 
reduces the channel gain for these signals even at their nominal link budget 
values. AGC then provides further gain reduction as these signals increase 
within their respective 3 dB power windows. Furthermore, because the MSC- 
3 and the EMU signals have relatively high initial SNRs the AGC action forces 
the thermal noise levels closer to the optical noise, thus degrading the margins 
of these signals by 0.8 and 0.2 dB, respectively. 


The scenario represented by Figure 2.4.2.2-1 is conservative from the standpoint of 
providing 12.5 dB of margin between the AGC threshold and the maximum signal 
modulation level (-2.7 dBm) of the FO driver. In an actual situation the 12.5 dB margin 
can probably be reduced to 10 dB or perhaps even to lower value. The actual margin value 
will depend on the number of simultaneous signals handled by the FO cable driver and on 
the relative magnitudes of these signals. 
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Figure 2.4.2.2-1. Optical Driver Output Levels for 3 dB Power Control Window 
(channelized received with AGC). 
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24.2.3 Possibilities and Limitations of Wider Power Windows 

In the preceding section, a 3 dB power window was assumed and the effect of such 
power window on a channelized receiver was considered. It was shown that with the AGC 
threshold set at 12.5 dB below the MSM level the FF and the OMV signals can be received 
without any margin degradation. It was also shown that the EMU and the MSC-3 signals 
can be received with 0.8 dB and 0.2 dB margin degradation, respectively. 

A power window wider than 3 dB can alleviate the power control requirements 
imposed on the user. It can also loosen the tolerances required of the AGC at the Space 
Station end of the link. Thus, the channelized architecture of the receiver and the 
availability of a channelized AGC naturally raise a question of whether a power window 
wider than 3 dB can be specified for this case. 

Providing for an operation with a wide power window can be obtained by a 
combination of the two following conditions: 

( 1 ) moving the AGC threshold level closer to the MSM level 
and 

(2) allowing the AGC action to move the thermal noise level closer to the level of 
the optical noise. 

Neither one of these conditions, however, can be realized without paying a penalty. 
Raising the AGC threshold level may result in amplitude clipping, particularly with multiple 
signals being applied to the FO cable driver. On the other, hand, allowing more gain 
reduction with the AGC decreases system output signal-to-noise ratio due to a higher level 
of the optical noise. Consequently, a wider power window must be a compromise between 
these two limitations. 

In considering the possibilities of using wider power control windows, we must 
also take closer look at the effect of the FO cable loss on the system output SNR. 
Depending on the magnitude of the FO cable loss various noise sources play the 
predominant roles. 
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As shown in Figure 2.4.2.3-1, there are three basic sources of the optical noise. 
One source is at the transmitter and two are at the receiver. At the transmitter the source is 
due to the relative intensity noise (RIN). RIN is inversely proportional to the cube of the 

transmitted power. However, once a given optical transmitted power is selected, RIN does 
not affect the C/N 0 at the receiving end of the link. This noise was already mentioned in 

Section 2.3.2. 

The noise sources which do affect the C/N 0 value at the receiver are the quantum 

noise and the thermal noise of the receiving device. The quantum noise degrades the 
receive C/N 0 ratio at the rate of 1 dB/1 dB of received power reduction, but the thermal 

noise provides a corresponding C/N 0 reduction at the rate of 2 dB/1 dB. The role which 
these noise sources play in degrading the C/N 0 ratio depends on the level of the optical 

signal impinging on the optical receiver. 

Figure 2.4.2.3-2 shows a typical behavior of the C/N 0 ratio for an optical link as 

the function of the received optical power (References [1] and [2]). As shown in this 
figure, for the received power levels in the 0 to -10 dBm, the RIN dominates the C/N 0 

ratio. In the region extending from -10 dBm to -20 dBm, the quantum noise determines 
the received C/N 0 . Below -20 dBm the optical receiver’s thermal noise determines the 

performance of an FO link. Thus, depending on the transmitted optical power and the FO 
cable loss different noise sources determine the system performance. 

It must be pointed out that the trend shown in Figure 2.4.23-2 is only a generalized 
trend, and that for each particular link, an analysis has to be performed to determine the 
actual performance. Based on the data obtained from the manufacturers of the optical and 
FO equipment, we have performed an analysis using a specific model of our link. 
Appendix A contains this analysis. Based on this analysis we have determined that at the 
maximum estimated FO cable loss of about 20 dB, the receiver thermal noise will dominate 
the link performance, but that at losses of 15 dB or less, the link performance will be 
dominated by RIN. 
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Figure 2A.2.3-2. Carrier/Noise versus Received Optical Power. 
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Figure 2.4. 2.3-3 shows the system SNR transfer curve as a function of the FO 
cable loss. From this curve we can see that for the maximum loss of 20 dB the "flattening" 
occurs at output SNR of about 35 dB. This, however, appears to be sufficient to receive 
the MSC-3 and the EMU signals with some link margins. It also appears that on top of this 
a 10 dB power control window may be implemented with some reasonable degree of 
accuracy for these high SNR users. 

Figure 2.4. 2.3-4 shows the same system SNR transfer curve in conjunction with 
the low SNR users such as an FF and an OMV. From this figure it is evident that for these 
low SNR users a power control window of up to 15 dB can be implemented with a 
reasonable accuracy and for the maximum cable loss of 20 dB. Also, these users are 
accommodated with some margin by the FO link. 

Consider now some of the limitations which may arise as the result of striving for 
wide power control windows. These limitations are particularly pronounced for the high 
SNR users such as the EMU and the MSC-3. From Figure 2.4.2.3-3 it evident that a 10 
dB increase in the received power moves the output SNR curve into a relatively flat, 
compressed region. Also, it can be seen that the flatness of the SNR curve is more 
pronounced for the case of the maximum FO cable loss of 20 dB. Thus, as the cable loss 
increases with time from say 10 dB to 20 dB a power increase of 10 dB at the system input 
will be detected as an increase in output SNR of only 5.9 dB instead of 9 dB. This fact is 
illustrated in Figure 2.4.2. 3-5. As shown there, the SNR threshold Th(10) set at a level 
when the cable loss is 10 dB will not be exceed when the cable loss increases with time. 

To alleviate this problem the algorithm for the power control must include 
provisions for compensating for such variation. The data for the compensation can be 
derived from an AGC amplifier placed after the photo diode at the receiver end of the FO 
link. The purpose for this amplifier is strictly to compensate for the gradual losses in the 
FO cable. Thus, this amplifier is not to be confused with other AGC amplifiers which 
form the conventional (i.e., the RF and the IF) portions of each Ku-Band channel receiver. 
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Figure 2.4.2. 3-3. System SNR Transfer Curve as Function of FO Cable Loss 
(EMU and MSC-3 Power Control Window of 10 dB). 
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Figure 2.4.2.3-4. System SNR Transfer Curve as Function of FO Cable Loss 
(FF and OMV Power Control Window of 15 dB). 
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Another significant issue encountered with large power control windows arises due 
to finite tolerances on the power control decision threshold. This is illustrated in 
Figure 2.4. 2. 3-6. As shown there, the normal operating point for an MSC-3 user output 
SNR of 27.0 dB. The corresponding margin at this SNR is 4 dB. As the signal power 
increases, the decision to command a 10 dB power reduction may occur at either a 6 dB 
power increase or at a 10 dB power increase depending on whether the threshold has 
moved by -2 dB or has remained at its proper value. If the decision is at the correct value, 
the power reduction will bring the link to its original operating point, i.e., SNR (in) = 27.8 
dB and SNR (out) = 27 dB. If, however, the decision is made at the Th(-2 dB) point then 
the link is brought back to a much lower operating point which is still within the margin 
region but is dangerously close to the lower end of the margin zone. 

The situation is worse if the power control window for the MSC-3 is 15 dB. 
Figure 2.4. 2. 3-7 illustrates this point. Here we see that if the decision is made at the 
Th(-2 dB) point, the power control command may force the link to operate below the 
margin. 

Based on the examples presented, one can, therefore, conclude that power control 
windows over 10 dB should be used with caution for the high SNR users which operate in 
the "flattened" portion of the output SNR curve. This, however, does not preclude the 
lower SNR users (see Figure 2.4.2.3-4) from utilizing 15 dB windows provided that the 
decision threshold tolerances are set properly. 
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Figure 2.4.23-6. A Loose Power Control Decision Tolerance May Push Control 
into the Margin Region 




User: MSC-3 
FO Cable Loss = 20 dB 
Power Window = 10 dB 



Figure 2.4.2.3-7. Excessive Power Window is Vulnerable to Decision Tolerance. 
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2.4.24 Power Control Algorithm for a Channelized Receiver 

The proposed algorithm for the power control is based on the measurement of the 
output signal-to-noise ratio (SNR) of the controlled channel. Such measurement does not 
depend on the absolute signal level and thus it is less sensitive to the residual inaccuracies 
of the AGC. 

Figure 2.4. 2. 4-1 shows a flow chart for implementing such algorithm for a 
channelized receiver. As shown in the figure, the first step is to decode the user's ID, i.e., 
FF, OMV, EMU, etc. The user ID permits the algorithm program to look up such user- 
related parameters as the SNR(min), the margin and the power control (PC) window. This 
data may be stored in an EPROM to permit changes in the user-related parameters. The 
assigned user PC threshold, i.e., the value of the SNR(out) at which the power control is 
activated is also obtained from the EPROM. Based on the measured FO cable loss the PC 
threshold is adjusted accordingly as described in Section 2.4.2.3. 

The steps described above should take place prior to the actual user acquisition by 
the Space Station. After these steps are completed, the PC program awaits the carrier lock 
flag. 

Once the carrier lock flag is received, the program begins to test the link SNR of the 
received signal. The test is based on comparing the pertinent signals at the output of the 
link demodulator. If the SNR falls below the SNR(min) value, the program alerts the 
antenna switching algorithm. But, if the SNR is within the acceptable limits, the program 
continues to measure the link SNR. 

When the link SNR exceeds the preset PC threshold, i.e., the SNR becomes high, 
a command is sent to the user to reduce the power level by the appropriate power window. 
The program then returns to the testing for the carrier flag. 
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Figure 2.4.2.4-1. Power Control Algorithm for a Channelized Receiver. 
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2.4.25 The SNR Estimation 

* 

Figure 2.4.2.5-1 shows a circuit for deriving the signals for the estimation of the 
channel SNR. The particular circuit shown is for a QPSK Costas loop, but a similar 
technique can be applied to higher order modulations such as 8-PSK. As shown, the actual 
measurement performed is that of (S+N)/N which is a close approximation to S/N for large 
values of SNR. 

The (S+N)/N component is derived by averaging the absolute values of the a(t) and 
b(t) data streams. The estimate of the N-component is obtained by taking the RMS value of 
the link noise in the carrier error track channel. The bandwidth of the latter channel is 
approximately equal to the bandwidth of the "sidearm" lowpass (LPF) filters. 
Furthermore, the highpass filtering eliminates the carrier tracking error component from the 
N-term. The SNR estimation can be performed by digital signal processing circuitry after 
quantizing the S+N and N signals to the required degree of accuracy. 
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Figure 2.4.2.5-1. QPSK Costas Loop Signals for SNR Estimation. 



2.4.3 User Power Control 

The commands sent to the user can alter user’s EIRP by either reducing the 
transmitter power, or the antenna gain or both. The simplest method is to control the 
transmitter power over a wide range of levels. 

The power amplifier baselined for the MA Ku-Band system users may be a 1 watt 
cw solid state amplifier. The power output of such amplifier may be controlled by varying 
the RF input drive by means of a PIN diode. As shown in Appendix C, the input power, 
and hence the transmitter output of the user, can be controlled over a range of 60 dB by 
means of a PIN diode. Appendix C contains the support material for defining the user 
transmitter configuration. 
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3.0 CONCLUSIONS 

In this report we have examined various aspects of establishing the requirements 
and the limitations for the power control for the Space Station Multiple Access System 
Users. Particular emphasis was placed on the performance of the fiber optic (FO) cable 
link which is currently baselined for connecting the Ku-Band receiver/transmitter units with 
the modem/processing equipment. Two receiver architectures were examined. These are: 
(1) a wideband receiver and (2) a channelized receiver. It was determined that a 
channelized receiver can, at least in principle, accommodate a 20 dB higher level of the 
input signals before the internal generation of the third-order intermodulation products. 
This capability is due to the channelized AGC available with the channelized receiver. 
Thus, we have assumed a channelized receiver as a desirable baseline for our system 
model. 

It was also determined that, based on practical model of an FO cable link, the basic 
limitation for the power control window is the FO cable loss which causes the degradation 
in the system SNR at the optical receiver. In view of this degradation, it is estimated that a 
power window for the high SNR users, such as EMU and MSC-3, may be limited to 10 
dB, at best. For the lower SNR users, such as FF and OMV, the power control window 
of up to 15 dB appears feasible. A prudent approach, however, is to set a 6 dB power 
control window as a design goal and allow a ±3 dB tolerance. 

The possibility of using a PIN diode to control the power input to the user's 
transmitter, and thus control the output power, was also examined. It was determined that 
60 dB of power control can be obtained in this manner without antenna switching. 
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4.0 RECOMMENDATIONS 

Based on the findings reported in this document, it appears that a channelized 
receiver system with channelized AGC can provide significant advantages over a wideband 
receiver which amplifies all of the received signals by means of a single IF. Thus, a 
channelized receiver is recommended. 

To allow for implementation tolerances, it is recommended that a power control 
window be kept below 10 dB, although some users can work with wider windows. A 
conservative baseline is 6 dB ± 3 dB window. Because the SNR performance of the FO 
cable link is limited by the optical receiver sensitivity when large (= 20 dB) FO cable losses 
are encountered a low noise optical receiver design is recommended. The receiver model 
used in this report is based on a PIN optical diode and a transimpedance amplifier. It is 
suggested that the use of an avalanche photo-diode (APD) receiver be investigated during 
the design phases of the MA link system. The use of an APD receiver may improve the 
noise performance at the receiving end of an optical FO cable link. 
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FIBER OPTIC CABLE LINK PERFORMANCE ESTIMATES 
1.0 INTRODUCTION 

In this appendix we provide a quantitative estimate of the FO cable link 
performance. The particular emphasis is placed on determining the effect of the noise 
sources which affect the performance of the link. 

As shown in Figure A-l, there are three basic sources of the optical noise. One 
source is at the transmitter and two are at the receiver. At the transmitter the source is due 
to the relative intensity noise (RIN). RIN is inversely proportional to the cube of the 
transmitted power. However, once a given optical transmitted power is selected, RIN does 

not change at the receiving end of the link. 

The noise sources which do affect the C/N fl value at the receiver are the quantum 

noise and the thermal noise of the receiving device. The quantum noise degrades the 
receive C/N 0 ratio at the rate of 1 dB/1 dB of the received power reduction, but the thermal 

noise provides a corresponding C/N 0 reduction at the rate of 2 dB/1 dB. The role which 
these noise sources play in degrading the C/N 0 ratio depends on the level of the optical 
signal impinging on the optical receiver. 

Figure A-2 shows a typical behavior of the C/N 0 ratio for an optical link as the 

function of the received optical power (References [1] and [2]). As shown in this figure, 
for the received power levels in the 0 to -10 dBm, the RIN dominates the C/N 0 ratio. In 

the region extending from -10 dBm to -20 dBm, the quantum noise determines the 
received C/N 0 . Below -20 dBm the optical receiver’s thermal noise determines the 

performance of an FO link. Thus, depending on the transmitted optical power and the FO 
cable loss different noise sources determine the system performance. 
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^General Optronics Corp. GO DIP 1000 Laser Module 
and/or GO ANA Laser Diode Transmitter are baselined 
for our analysis. 




CARRIER/NOISE (dB/Hz) 



The trend shown in Figure A-2 is only a generalized trend. To obtain a more 
accurate estimate of the performance of our link, we have to make some reasonable 
quantitative assumptions and perform the required calculations. 

For the calculations that follow, we assume that a General Optronics Corp. 
GO DIP 1000 laser module is used at the transmitting end of the link. This module may be 
incorporated into the GO ANA analog laser diode transmitter. We also assume that the 
laser diode, which is rated at 3 mw cw power, is biased at 2.25 mw (i.e. P b = 2.25 mw) 

and that the total modulation index is 50% for a 5-channel link. We also assume a PIN 
diode and a transimpedance amplifier with an equivalent resistance of IkO (R eq = lkQ) and 
a noise figure of 3 dB (F t = 2) is used at the receiving end. 

Appendix B contains the data sheets on the General Optronics Corp. laser 
transmitter units. 
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2.0 CALCULATIONS 


The C/N 0 value at the receiving end of an FO cable link can be expressed by the 

(mR»P b ) 2 


following equation [1] 
1 


C ^o = (W) 


(RIN Rj P 2 b B) + 2q(R 0 P b + I d )B + (4KTB/R eq )F t 
source quantum receiver 


(A-l) 


where 

N = number of independent channels 

R 0 = diode responsibility (amperes/watt) 

m = modulation depth (total, all carriers) 

P b = average optical power received (w) 
q * electron charge (1.6 x 10 -19 coulombs) 

K = Boltzman constant (138 x 1(T 23 watt-sec/°K) 

I d = diode dark current (amperes) 

B = bandwidth of receiver (Hz) 

T = temperature (°K) 

R eq = equivalent resistance of photodiode load and amplifier (ohms) 

F t = noise factor of preamplifier 

RIN = source relative intensity noise (dB/Hz in numeric value) 

The three terms in the denominator of equation A-l identify the aforementioned 
contributions of the source, the quantum and the receiver noise terms. If we let B = 1 Hz 
and neglect the effect of the diode dark current (i.e., let I d = 0 amp) we can separate 
equation A-l into three separate C/N 0 terms as shown below: 


Source Noise: 


C/N 


1 m 2 
0 _ 2N RIN 


(A-2) 
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Quantum Noise: 


Receiver Noise: 


1 ( R o P b> 

C/N 0 = — (m2)— — 
0 2N 2q 


i , ( R o p b) 2 K, 

C/N n = — (m2) 

0 2N (4KT) F t 


For our baseline model we assume the following values: 

N =5 (Five channels) 
m = 0.5 (total for five channels) 

RIN = -131 dB/Hz (7.94 x 10~ 14 numeric) 

R 0 = 0.5 A/w 

P b = 0.00225 w (2.25 mw at transmitter) 

T = 290°K 

R = 1KQ 

eq 

F t = 2 (NF = 3 dB) 

I d = negligible 

Using these values we proceed to calculate the various C/N 0 terms. 


Source Noise C/N 0 : 


C/N„ = T^77*<0.5)V 


1 


0.25 


(2) (5) ~ 7.97 x 10“ 14 10 x 7.94 x 10" 14 

= 3.147 x 10 11 or 115 dB/Hz 
Quantum Noise C/N fl : 

q = 1.6 x 10 -19 coulombs 

2.81 x IQ" 4 


C/N 0 = — i — x (0.5) 2 x 

U \ /c\ K 7 / 1 n\ „ i n-1" 


(2) (5) 


(2) (1.6) x 10 y (10) (3.2) x 10 


-19 


= 8.19 x 10 or 139.4 dB/Hz 
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(A-3) 


(A-4) 



The value above is for no cable loss. For a 20 dB loss the C/N 0 value becomes: 


C/N 0 (L = 20 dB) = 139.4 dB/Hz - 20 dB = 119.4 dB/Hz. 


Receiver Noise C/N 0 : 

Assume 

r = 1000 Q. 

eq 

T = 290°K 

K = 138 x 10~ 23 watt-sec/°K 

(KT = 4 x 10" 21 ) 

F = 2 (NF = 3 dB) 


Then for zero cable loss (i.e., L = 0 dB): 


c/N o = TN (m2 ) 


(RpRb) 2 r 
4 K T F 


1 

" (2) (5) 


2 (0 5) 2 (2.25 x IQ' 3 ) 2 
(4) (4 x 10" 21 ) 


10 3 


(0.25) (0.25) (5.1) (IQ"*) 10 3 
(10) (4) (2) (4 x 10 _il ) 


= 9.96 x 10 14 or 150 dB/Hz 


For L = 20 dB this value will decrease by 40 dB: 

C/n q = (L = 40 dB) = 150 dB/Hz - (2 x 20 dB) = 110 dB/Hz 

Figure A-3 shows the three respective C/N 0 terms as the function of the cable loss. The 
combined value of the three C/N 0 terms is also plotted in Figure A-3 as the function of the 
cable loss. 

Using the combined C/N 0 value for the optical link, the SNR transfer curves can be 
plotted as the function of the system input SNR and of the FO cable loss. For this the 
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Figure A-3. Carrier/Noise Ratio per Channel as Function of FO Cable Loss. 


C/Nq value is converted to the nominal video link bandwidth of 25 MHz. Figure 2.4.2.3-3 
and 2.4.2. 3-4 of the report show this conversion. 
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APPENDIX B 


DATA SHEETS FOR GENERAL OPTRONICS CORP. 
ANALOG LASER DIODE TRANSMITTER UNITS 




820 nm 
MULTIMODE 
LASER MODULE 
GO DIP lOOO 

FEATURES: 

• High power/efficient coupling into fiber 

• Temperature and power monitoring 

• Low modal noise 

• Low self-pulsation noise 

• High linearity 


DESCRIPTION: The GO DIP 1000 is a 14 PIN dual- 
in-line package incorporating our GOLS 1000 injec- 
tion laser diode, coupled to a 50^m core fiber. A 
photodetector and a thermistor are included for 
power and laser temperature monitoring. An optional 
thermo-electric cooler can be used in conjunction 


with the thermistor for temperature stabilization and 
control. The entire assembly is hermetically sealed, 
containing no organic material Thousands of units 
are now operating worldwide in analog CATV 
distribution, digital telephone transmission and high 
speed data link systems. 


SPECIFICATIONS (@ 25°C): 


Parameter 

Symbol 

Min. 

Threshold Current 

Uh 


Light Output 

L 

3 

Quantum Efficiency 

. n 

0.08 

Spontaneous Emission 

Lth 


Wavelength 

X 

800 

Spectral Width 

AJt 

1.0 

Rise/Fall Time 

T 


Forward Voltage 

V, 


Monitor Current 

|d 

0.4 

Thermistor Resistance 

nt 

9.5 

Linearity 

2nd harmonic distortion 



3rd harmonic distortion • 

Fiber Pigtail (1 meter) 

* 


core diameter 


47 

cladding diameter 


122 

jacket diameter 



numerical aperture 


0.18 

Cooler Capacity 

A ^max 

40 

Cooler Current 

lc 



Typ. Max. Unit Remarks 

80 100 mA 

mW out of fiber core 
0.20 mW/mA in lasing region 


250 


fj W 

820 

850 

nm 

1.5 


nm 

0.5 


ns 

1.8 

2.0 

V 

12 

25 

mA 

10.0 

10.5 

k Q 

-35 


dB 

—45 


dB 


L at l th 

other X available 
30% CW 0 1 5mW 

9 L=1.5mW 
L * 1 .5mW 


50% modulation depth 
50% modulation depth 


50 

125 

950 

0.20 

0.9 


53 

128 

0.22 

1.1 


fim customer supplied fiber can be pigtailed on request 

f<m 

°C optional 

A at AT max 


ABSOLUTE MAXIMUM RATINGS: 

• Storage Temperature: — 20°C to +60°C 

• Operating Temperature: 

— 20°C to +60°C 

• Soldering Time: 1 0 seconds at 260°C 

• Reverse Voltage: 5V 

• Photodetector Reverse Voltage: 30V 

• Light Output: 5mW 


ORDERING INFORMATION: 

Model Number: GO DIP 1000 
Wavelength Specifications: Typical range is 800nm to 
850nm. Specify if you require: a) wavelengths outside 
this range; b) tight tolerance on wavelength. 

Options: a) customer supplied fiber; b) thermo-electric 
cooler (TEC). 


B-l 


ORIGINAL PAGE IS 
OF POOR QUALITY 


OflBBOBBBBBBBBBIU 




general optroniacorp. 

Two Olsen Avenue. Edisoa NJ 08820 
(201) 549-9000 TWX 710-997-9556 

ORIGINAL PAGE S3 B-2 
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SE qeneral optroriic/ corp. 



ANALOG 
LASER DIODE 
TRANSMITTER 
GO ANA 

FEATURES: 

• Wide electrical Bandwidth 

• Power and Temperature Control 

• Alarm Circuitry 

• Uses Hermetic DIP Package, GO DIP Series 


DESCRIPTION: This general purpose analog driver 
can be ordered with any of our DIP packaged laser 
diodes. The circuit incorporates all functions in order 
to monitor the optical power bias and the operating 
temperature. The transmitter can accommodate all 


types of analog or digital signals (with 50% duty cycle). 
Power and temperature are factory set depending on 
the laser module used inside the transmitter. 

A detector circuit enables an alarm and monitors the 
operation of the device. 


ABSOLUTE MAXIMUM RATINGS: 

• Input Signal: IVp-p. 

• Supply Voltage: -7V 

• Storage Temperature: — 50°C to +80°C 

• Operating Temperature: -40°C to +65°C 
SPECIFICATIONS (@ 20°C): 

• Power Supply: -5V, +/- 5%, 1.2A max. 

• Input Impedance: 50-75Q 

• Stability = 1% at 25°C, 3% over temperature 
range 

• RF Input: SMA 


ORDERING INFORMATION: 

GO ANA XXXX 

XXXX suffix from GO DIP light source 1000, 1300, 
3000, 4000, 6300 

Temperature setting is 25°C for transmitters 
incorporating GO DIP 1000, 3000 and 4000, and 
20°C for transmitters incorporating GO DIP 1300 
and 6300. Power bias is 1.5mW for multimode fiber 
pigtail, 500pW for 9pm single mode pigtail, 300pW to 
ImW for 5^m single mode pigtails. 
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APPENDIX C 


POWER/I ERP CONTROL FOR KU-BAND MA SYSTEM 
USER S TRANSMITTERS 



TECHNICAL MEMORANDUM 


TM 8706-2 


TO: 

FROM: 

SUBJECT: 

DATE: 


Sergei Udalov 
Richard Austin 

Power/IERP Control for Ku-Band MA System User's Transmitters 
June 9, 1987 


1.0 INTRODUCTION 

The purpose of this memo is address the issue of power/EIRP control for the Space 
Station (SS) Ku-Band multiple access (MA) system users. Due to wide variations in user 
ranges and antenna gains, the dynamic range of the received signal strength could be quite 
large. This may necessitate the implementation of a system for controlling the users 
transmitted power or EIRP. This memo examines an approach to the problem from both a 
system and hardware perspective. 

2.0 PROBLEM STATEMENT 

Users of the Ku-Band MA system fall into one of four categories. They are listed 
here along with the approximate gain of their antennas: 

( 1 ) Free Fliers (FF) — 30.3 dB 

(2) Orbital Maneuvering Vehicle (OMV) — 13.3 dB 

(3) Extravehicular Mobility Unit (EMU) — 3 dB 

(4) Mobile Service Center (MSC) — 3 dB 




Additionally, they will have a transmitted output power baselined at 1 watt cw (BPSK, 
QPSK, or 8-PSK*) and an operating frequency range of 14.5 - 14.76 or 14.6 - 14.9 
GHz. The difference in the transmitted EIRP of these users along with the wide variation 
in range will result in a large dynamic range of received signal strength at the SS receiver 
front end. Since a single receiver must accommodate all MA users, gain control at the 
receiver is ineffective when signals are received from more than one user at the same 
time.** As a result, some means must be employed to control user EIRP. This may take 
the form of a reduction in transmitter output power, switching to a lower gain antenna, or 
both. A quantitative analysis of the problem has shown that at least 30 - 40 dB of 
power/EIRP control may be required to maintain received signal strength levels within an 
acceptable range. 


3.0 SYSTEM APPROACH 

Switching between antennas as a means of reducing EIRP results in abrupt changes 
in the received signal power. Additionally, if only two antennas are employed, only two 
EIRP values are obtainable resulting in coarse power control. Coupled with the expense of 
engineering multiple antennas, it makes this approach the least attractive. The 30 - 40 dB 
reduction in EIRP is most easily obtained by varying the transmitter output power. 

An approach to this solution is shown in Figure C-l. The user transmitter output 
power is controlled by varying the driver signal to the solid state amplifier with a PIN diode 
attenuator. A measurement of the received signal strength for a particular user is made by 
the SS signal processor. If the signal strength falls below or rises above some 
predetermined value, a command to change the transmitter output power accordingly is 
generated and modulated onto the SS to user command link. The signal processor in the 
user decoders the command and sends the appropriate control signals to the attenuator. The 


* May require 2.5 watts 

This statement applies to a non-channelized, wideband receiver. 
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result is similar, a closed loop AGC system. How often the attenuator must be updated and 
with what granularity would have to be determined based on a study of the rate of change 
of signal strength due to user motion. 

4.0 HARDWARE CONSIDERATIONS 

It is worthwhile to examine current device technology and identify possible 
components that would allow control of the user transmitter output power. A quick review 
of the requirements shows: 

(1) User Output Power = 1 watt cw 

(2) Operating Frequency Range = 14.5 — 14.76 or 14.6 - 14.9 GHz 

(3) 30 to 40 dB of EIRP/power control 

In addition, we require that the 30 - 40 dB of EIRP/power control be achieved by varying 
the input power to the power amplifier in a linear or quasi-linear manner. 

The input-output characteristics of current solid state amplifiers follow the 
traditional "S" curve found with tube and other types of amplifier, i.e,. a thermal noise 
floor transitioning to a linear region transitioning into a saturated region. The linear region 
extends from the noise floor to the 1 dB compression point. This is the point where a 10 
dB increase in input signal power results in only a 9 dB increase in output signal power. 
The dynamic range of the linear region can be quite large — often greater than 70 dB. 

It's instructive to look at a specific example. Avantek makes a series of 
communication satellite uplink driver amplifiers (see attachment) that are typical of devices 
on the market whose performance characteristics closely match the stated requirements. 
Model #AWP 145101 has the following characteristics: 

Frequency Range: 14.0 - 14.5 GHz 

Output Power: 30 dBm (1 watt) 

Gain: 30 dB 
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Noise Figure: 


10 dB 


The thermal noise floor of this device is determined by the following: 

P T = NF - 114 dBm + 10 log Bw 

where 

P T = noise power in dBm 
NF = noise Figure =10dB 

Bw = operating bandwidth in MHz = 300 (14.9 - 14.6 GHz) 

P T = 10 dB- 114 dBm + 10 log 300 

= 10 dB - 1 14 dBm + 24.8 dB 
P T = -79.2 dBm 

If the 1 dB compressed gain maximum output power of the device is 30 dBm and the gain 
30 dB, then the maximum input signal power is 0 dBm. Hence, the dynamic operating 
range of the linear region is approximately 79 dB (0 dBm - 79.2 dBm) 

While step attenuators could be used to vary the input (and hence the output) power 
of the signal, given the linear characteristics of the amplifier, the use of continuous 
attenuators such as PIN diodes is more appropriate. The use of step attenuators can result 
in phase discontinuities in the signal when switched. Avantek packages many of their 
amplifiers with a PIN diode attenuator as part of the first stage of the amplifier. Control of 
the amount of attenuation is obtained by varying the bias voltage to the diode. This can be 
done mechanically by means of a potential meter mounted on the case, or electrically by 
means of a gain control circuit. Depending on the frequency, a single PIN diode stage can 
yield 30 dB of attenuation. If more attenuation is required, additional stages may be 
connected in series. For example. General Microwave Corporation’s Model 3250 digitally 
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controlled PIN diode attenuator provides up to 60 dB variable attenuation within the band 
of interest (see attachment). 

An alternative, but generally less desirable method of changing the output power is 
to vary the DC bias voltage to the output stage of the amplifier. This has the adverse effect 
of changing many of the parameters of the amplifier which, depending on the specific 
application, may or may not be a problem. Typical changes that result in varying the DC 
bias voltage are an increase in the noise figure, lowered gain, reduced linear operating 
region, and degradation of the AM-PM conversion characteristics. 

It will be noted that the communications satellite uplink driver amplifiers listed in 
the attached sheet and used as an example in the section have an operating frequency range 
of 14.0 to 14.5 GHz while the requirements specify 14.5 to 14.9 GHz operating range. 
Discussions with manufacturer representatives have verified that the performance 
characteristics listed would just as easily be obtained for devices in the required operating 
band. As the attached sheet shows, the amplifies listed are all designed to meet an existing 
application. The lack of devices in the 14.5 to 14.9 GHz range is only a reflection of the 
lack of a demand for such a product and not of any technical difficulties in achieving them. 

5.0 CONCLUSIONS 

We have presented a variable approach to the user power/EIRP control problem. 
While effective, it has the drawback of adding additional complexity to the system design 
of both the user and the SS. A simpler alternative might be to have the users control their 
output power based solely on their range to the SS. This would result in a simplified open 
loop system not involving the SS. 
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MEDIUM POWER, COMMERCIAL COMMUNICATIONS BAND AMPLIFIERS 

] Guaranteed Specifications @ 25° C Case Temperature pcs 

| Power Output Typical 

| Frequency toe 1 dB Cain Gain Intercept Nolae VSWq 

j Retponae Compre»»lon Gain Flatncaa Point lor Flyura (50 ohms) 

j (MHi) (dBm/Watts) (dB) (idB) IM Products (dB) Maximum 

! Model Minimum Minimum Minimum Maximum (dBm) Maximum In Out 


3700 to 4200 MHz Point-to-Point Communications Amplifiers 


AWP-42107 _7- 

. 3700-4200 

7: +25/3 

26 _ ; 


^+35_ 

10 

71-2 . 


soo;2aJ 

AWP-42108 

3700-4200 

+32/1.5 

33 

2 

+40 

10 

1.2 

1.2 

+15 

1200 

AWP-42109 

3700-4200 

+40/10 

41 

2 

+47 

10 

1.2 

12 

+ 15 

6000 


5925 to 6425 MHz Point-to-Point Communications or Satellite Uplink Amplifiers 

AWP-64107 5925-6425 +25/3 26 2 +34 10 1.2 1.2 +15 500 

AWP-64108 ^175925-6425 7^32/1 5 23^^:. 2 "£2^*40 77 1 2 

AWP-64109 5925-6425 +40/10 41 2 +46 10 1.2 1.2 +15 6000 


6425 to 7125 MHz Point-to-Point Communications Amplifiers 


AWP-71107 

AWP-71108 

6425-7125 

6425-7125' 

+2S/.3 

+32/1.5 

26 

33 

2 

2 

+34 

+40 

10 

10 

1.2 

1.2 

1.2 

1.2 

+15 

+15 

500 

1200 

6000^1 

AWP-71109 

6425-7125 

+40/1 0" 1 "" 3 


' , ""~+46 x ' rr *~' 

“10 ^ 


^1.2 


7125 to 7725 MHz Point-to-Point Communications Amplifiers 

AWP-77107 

7125-7725 

+25/. 3 

26 

2 

+34 

10 

1.2 

1.2 

+15 

500 

AWP-77108 7125-7725 “7 

r::+32/i .5~~~7' 33 ~ n ^ rz " r ‘. t~7 +40 77T7' 

jo 71 


37^.2 


12003 

AWP-77109 

7125-7725 

+40/10 

41 

2 

+46 

10 

^1.2 

1.2 

+ 15 

6000 


7725 to 8500 MHz Point-to-Point Communications Amplifiers 


AWP-85107 

7725-8500 

_ +25/.3 

^ 26 

_2 ^ 

*^♦34 _ 

__ IQ....-;, „1.2 ^ 


+15 - ~ 500 

AWP-85108 ~ 

7725-8500 

+32/1.5 ~ 

33 

2 

+40 

~~ 10~ 1.2 ~ 

1.2 

+15 1200 

AWP-85109 

7725-8500 

+40/10 

41 

2 

+47 

10 1.2 

1.2 

+15 6000 


7725 to 8275 MHz Point-to-Point Communications Amplifiers 


AWP-83107 

7725-8275 

. +25/3 

26 2 

+34 

10 

1.2 

1.2 

+ 15 500 

AWP -831 08 7^7“ 

'7725-8275 

7j7+32/i.s;^ 


-t40 •” 

10 *7^7' 

1 2™T 

,, T 1 .2~T 

TT;+ 15 ^77;; i 2 00 1^7; 

AWP-83109 

7725-8275 

+ 40 / 10 "'^ 

*"*41 - 

+47 

Yo~ 

1.2 ' 1.2 ' 

" '%15 ***“"* 6000 ^ 


10.7 to 1 1.7 GHz Point-to-Point Communications Amplifiers 


AWP-1 171 07 

10.7-11.7 

+25/.3 

26 

2 

+32 

10 

1.2 

1.2 

+15 

500 

AWP-1 17108 

10.7-11.7 

+32/1.5 

33 

2 

+39 

10 

1.2 

1.2 

+ 15 

1200 

AWP-1 17109 

"’10.7-1 17^"* 

~ +37/5 

^ 38 r 

"2 ” 

r ~*~'+43 

^10 ' 

1.2 

" tr! ~1.2 

^ rc ^+1 5 

^6000' 


14.0 to 14.5 GHz Communications Satellite Uplink Drivers 


In pul Power 

Current 

Voltage (mA) 

(VDCi3%> Typical 


AWP- 145001 14.0-14.5 +24Z.25 20 .25/40 MHz +30 10 1.2 • - 1.2 ' ' +15 400 

AWP-1 45001 _ . 1 4.0-14.5 +27/.S 20 . .25/40 MHz +33 „7-, - 10 W7, 1.2 j. 7 1.2 . ^ ^ . ..+15 ..^7 900 -3 

; AWP-145101 14.0-14.5 +30/1.0 * 30 .25/40 MHz +36 10 1.2 i’.2 +15 **1000 

AWP-145102 14.0-14.5 +33/2 30 .25/40 MHz +39 10 1.2 1.2 +15 4500 

[, AWP-145103' '“"14.0-14.5 ’^" '♦34/3 '”'^ 30 ' .25/40 MHz 7+40 10“'"T r "'l.2 4500'7 r3 

j AWP-1 45503 " 14 0-14.5 "“7' +34/3 ,17. 40 ^ .25/40 MHz .77 r +40 10 .TT',. 1.2 1.2 77 7*15 '77 4500 7 • 

AWP-1 45504 14.0-14.5' +36/4 ' 40 .25/40 MHz ’ - ‘^“+4l“‘ " 10 "'“'** *1.2 ’ 1.2 ‘ ~ +15" '“6000 

AWP-1 45505 14,0-14 5 +37/5 40 .25/40 MHz +42 10 T2 V2 +15 6000 

| For case drawings, contact Avantek, Inc., 180 Blue Ravine Road. Folsom, CA 95630-Phone (916) 985-1201. 


Note 1: Units are available with or without top and/or bottom heatsinks. 
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PHASE SHIFTERS • 


Model 3250 Ultra-Broadband 
Digitally-Programmable PIN Diode Attenuator 


• Frequency range: 0.2 to 18 GHz 

• Attenuation range: Up to 60 dB 

• Binary or BCD programming 

• Absorptive 

• Guaranteed Monotonicity 



The Model 3250 digitally programmable attenuator 
provides excellent performance characteristics over 
the frequency range of 0.2 to 18 GHz. Attenuation 
levels up to 60 dB are programmable in increments 
of 1 dB. 

The unit is an integrated assembly of a dual T-pad 
PIN diode attenuator and a driver consisting of a 
D/A and a V/l Converter. See figures 1 and 2. 

The Model 3250 operates as a bilaterally-matched 
device at all attenuation levels. It is supplied in a 
compact rugged package well-suited to military 
applications. 
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Fig . 1 -Model 3250, block diagram 
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Fig. 2-Model 3250, rf schematic diagram (unit 
consists of two such sections) 
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Model 3250 
Specifications 


PERFORMANCE CHARACTERISTICS 


Frequency Range 

.0.2 to 18 GHz 

Mean Attenuation Range 

0.2 to 12.4 GHz 

.60 dB 

12.4 to 18 GHz 

.50 dB 

Insertion Loss (max.) 

0.2 to 8 GHz 

.2.5 dB 

8 to 12.4 GHz 

.3.0 dB 

12.4 to 18 GHz 

.5.0 dB 

VSWR (max.) 

0.2 to 8 GHz 

. .1.75 

8 to 12.4 GHz 

. .2.0 

12.4 to 18 GHz 

. .3.0 


Accuracy of Attenuation (1) . ±0.3 dB 

Temperature Coefficient. . .0.02 dB/°C max. 

Power Handling Capability 

Without Performance 

Degradation From 0.4 to 100 mW cw or 

peak (see Figure 3 below) 


Power Handling Capability (con't) 

Survival Power (from -65°C to +25°C; see power 
derating curve, Figure 4 below, for higher 

temperatures) 2W average or peak (1 psec 

max pulse width) 

Switching Time 20 psec max. 

Programming Positive true binary (stan- 

dard) or BCD 
(option 1). For 
complementary code, 
specify Option 2. 

Minimum Attenuation Step . 1 .0 dB 

Logic Input 

Logic 0 -0.3 to + 0.7V 

Logic 1 +2.5 to +5.0V 


Power Supply 
Requirements . 


Power Supply Rejection. 


. +5V ±5%, 310 mA 
+ 15V ±5%, 100 mA 
-15V ±5%, 100 mA 

.Less than 0.1 dB/volt 
change in any supply 


FLATNESS 



Frequency (GHz) 

ATTEN. 

(dB) 

0.2 to 8.0 

0.2 to 12.4 

12.4 to 18.0 

10 

0.5 

0.7 

1.0 

20 

0.5 

1.0 

1.0 

30 

0.7 

1.5 

1.5 

40 

1.0 

1.5 

1.5 

50 

1.0 

1.5 

1.5 

60 

1.0 

1.5 

1.5 


( 1 ) At calibration frequency. All units are calibrated at 4 GHz. 
Calibration at other frequencies within the band is available upon 
request. 
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RATING 
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♦25°C 

TEMPERATURE 


Fig. 3-Modet 3250, maximum peak and average 
operating power without performance 
degradation. 
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Model 3250 


ENVIRONMENTAL RATINGS 

Operating Temperature 
Range — 54°C to +110°C 

Non-Operating 

Temperature Range. . . . -65°C to +125°C 

Humidity M1L-STD-202F, Method 103B, 

Cond. B (96 hrs. at 95%) 

Shock MIL-STD-202F, Method 213B, 

Cond. B (75G t 6 msec) 

Vibration MIL-STD-202F, Method 204D, 

Cond. B (.06" double amplitude 
or 15G* whichever is less) 

Altitude MIL-STD-202F, Method 105C, 

Cond. B (50,000 ft.) 

Temp. Cycling MIL-STD-202F, Method 107D, 

Cond. A, 5 cycles 


AVAILABLE OPTIONS 


PIN FUNCTIONS 


Option No. Description 

1 BCD programming (Binary is 

standard 

2 Complementary programm- 
ing (positive true is 
standard) 

7 Two SMA male rf 

connectors 

10 One SMA male and one 

SMA female rf connector 


PIN NO. 

BINARY 

BCD (Opt. 1) 

_ 1 

SPARE 

SPARE , 

2 

SPARE 

■'“SPARE 

3 

+ 5V V 

W+sv ;^T 

4 

DIGITAL & 

DIGITAL & : 

' Jr . 

POWER GND 

POWER GND 

' 5 

• GND V 

-r- 1 dB ■':£* 

6 . 

. GND ..v 

2 dB ~>r 

7 

■*v- i dB . 

:.--StS 4 dB -V 

8 

V- 2 dB J • 

•r?M : e dB y*? 

9 

. 4 dB " 

10 dB 

10 

8 dB 

■ 20 dB *: 

11 

16 dB ' 

40 dB 

12 

32 dB 

OPEN (NO 

. . ... 


CONNECTION) 

13 

+ 15V 

_,+ 15V 

14 

-15V 

. -i5v 

15 

SPARE 

SPARE . . 


ACCESSORY FURNISHED 

Mating power/logic connector 


DIMENSIONS AND WEIGHTS 


ITT CANNON OA15POH EOUIV 
WITH 01 10551 JACKPOSTS 
MATING CONNECTOR FURNISHED 


136 13 .51 DIA THRU 
4 HOLES 


" — i so n ' *o / 

110 , <2.5)/ I*. 

•—I4S.7I— —7 — <»•* ^ 

\ / i "1 

m in / i 111 


MTG 

SURF ACS 


SMA MALE 
FOR OPTION 10 


L 1 


J6 19 71 FOR SMA FEMALE RF CONNECTOR 
.» (12.7) FOR SMA MALE SMA FEMALE (3) 


MODEL 3250 
Wt: 4 oz. (113 gm) approx. 


Dimensional Tolerances, unless otherwise indicated: .XX ± .02; .XXX ±£05 
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1.0 INTRODUCTION AND OVERVIEW 

This report presents a baseline design of a universal 4/8-PSK modem for potential 
utilization by the Space Station Ku-band, multiple-access system users. Although the 
present concept of the return links can be satisfied by a conventional 4-PSK (i.e., QPSK) 
modulation providing either two independent channels, or one channel coded at rate 1/2, 
the motivation to consider the 8-PSK mode results from the MSC requirement for 3 digital 
television channels. Thus, the idea expressed in this report is that a 4-PSK modem design 
baselined for the major portion of users can be expanded to include an 8-PSK capability at 
a reasonable cost in increased complexity. In other words, the intent here is to assume a 
4-PSK design as a baseline which is easily expanded to include an 8-PSK capability. 
Ideally, the component partitioning would be such that a major portion of the modem will 
consist of the 4-PSK capability, with the 8-PSK capability provided in the form of plug-in 
modules for the MSC. user as well as for other users which may require 3 digital TV 
channels in the future. 

The design considerations presented in this report are driven by the requirement to 
transmit and receive digital TV signals at rates up to 22 Mbps - 25 Mbps per channel. 
Thus, circuit configurations which can be implemented by high speed circuitry are given 
primary consideration. 

The requirements for the modem in terms of data rates are summarized briefly in 
Section 2. Modem implementation is outlined in Section 3. In Section 3.1 modulator 
configurations are presented. Section 3.2 then describes a design for a corresponding 
demodulator. Section 4 presents conclusions. Appendix A contains a simplified analysis 
of modem operation in the 4-PSK mode. 
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2.0 MODEM REQUIREMENTS 

The return link requirements for Ku-band MA users are summarized in Table 2-1. 
From the table, it can be seen that there are basically two return rates; one at about 
100 kbps and the other at about 22 Mbps. Although, as shown in the table, the two 
Phase B contractors established slightly different requirements for these two rates, the 
implementation driver for the modem is still the upper rate which may be as high as 
25 Mbps according the RI estimate. 

The functional goals which are the motivators for considering an 8-PSK modem to 
meet the MSC requirement for 3 simultaneous digital TV channels are as follows: 

1) Addition of the three-channel capability must not increase RF channel 
bandwidth significantly beyond that for the two-channel mode. 

2) Modulation method should not change drastically when the modem 
switches from 2-channel to 3-channel mode. 

3) The modem must have as much hardware commonality as possible 
between the 2-channel and the 3-channel modes. 

Figure 2-1 shows a functional diagram of a 4/8-PSK multi-channel modem 
utilization within the Ku-band MA system. As indicated, the 4-PSK (i.e., QPSK) mode is 
the baseline, and the 8-PSK mode is considered as a multi-channel capability. This 
diagram is responsive to the goals stated above. One of the salient features of the concept 
shown in Figure 2-1 is that a constant envelope signal is provided by the 8-PSK operation 
making it fully compatible with the RF amplification equipment used with the 4-PSK 
mode. 

Figure 2-2 shows a functional hardware partitioning for the proposed 4/8-PSK 
modulator and the corresponding demodulator. The main idea expressed in Figure 2-2 is 
that the 4-PSK mode is the baseline mode for the modem and that the 8-PSK mode is a 
hardware "add-on" to be utilized by such users as the MSC. Such partitioning requirement 
determines the modem implementation/configuration design described in the subsequent 
section of this report. 
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Requirement Rates 


User 

TL M/Voice (kbps) 

Video (Mbps) 

Comments 


MCDD 

RI 

MCDD 

RI 


MSC 

160 

128 

22 

25 

Three simultaneous video 
signals are required. 

EVA 

100 

128 

22 

25 


NSTS 

100 

128 

N/A 

N/A 


OMV/OTV 

100 

128 

22 

25 


FF/COP 

100 

128 

22 

25 



MCCD = Requirement established by McDonnell Douglas Corp. team [1]. 
RI = Requirement established by Rockwell International team [2], 


Table 2-1. Return Link Requirements for Ku-band MA-Users. 
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Figure 2-2. Functional Hardware Partitioning for the 4/8 PSK 
Modulator (a) and Demodulator (b). 
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Because the subject matter of this report deals mainly with the implementation of the 
modulation/demodulation functions of the proposed modem, we do not address the issues 
of an IF frequency at which the actual modem implementation should take place. It suffices 
to state, however, than an IF frequency in the range of 150 MHz to 700 MHz is 
envisioned. The final selection of the IF frequency will be determined by the frequency 
plan of a particular Ku-band MA system selected for the Space Station/user application. 
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3.0 MODEM IMPLEMENTATION 

3.1 Modulator Configurations 

3.1.1 Configuration for 4/8-PSK with Two Quadrature Modulators 

Because the main idea for a universal modem is to provide a maximum amount of 
component commonality between the 4-PSK and 8-PSK modes, let us consider first an 
implementation which uses two quadrature modulators to extend a 4-PSK capability to that 
of an 8-PSK. For this we first consider a single 4-PSK (i.e., QPSK) modulator 
implemented with two orthogonal bi-phase modulators. 

Part (a) of Figure 3. 1.1-1 shows the configuration for such a quadrature modulator 
implemented with two bi-phase modulators whose phases are orthogonal to each other. As 
shown in part (a), the bi-phase modulator #1 is driven by baseband data stream a(t), and 
bi-phase modulator #2 is driven by baseband data b(t). The resultant signal constellation is 
shown in part (b) of Figure 3. 1.1-1. 

The positions of the resultant QPSK vectors shown in (b) have been arranged 
according to the Gray coding technique. With this technique a 90° phase error in the 
demodulator output results in only a one-bit output error at the time of the error. With other 
coding schemes both output bits will be in error when a 90° wrong decision is made at the 
demodulator. This is an important point, because in the subsequent implementation of an 
8-PSK modulator, the phase constellation representing positions of 3-bit symbols will be 
also encoded with Gray code. 

Figure 3. 1.1 -2 shows a functional diagram of 4/8-PSK modem implemented with 
two quadrature modulators. Each of these modulators is like the quadrature modulator 
shown in (a) of Figure 3. 1.1-1. Thus when the modem is in the 4-PSK mode, only 
modulator #1 is operational. In this modp it accomodates two baseband channels a(t) and 
b(t). These can be either two independent channels or two components of rate 1/2 encoded 
single data stream. 
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Carrier Phase 
to Modulator #2 



Figure 3. 1 . 1 - 1 . Quadrature Modulator Implementation with Two 
Bi-Phase Modulators 
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Reference Carrier 



Figure 3. 1.1 -2. Functional Block Diagram of 4/8-PSK Modem Implemented 
with Two Quadrature Modulators. 



For the 8-PSK mode quadrature modulator #2 is activated in addition to modulator 
#1 and the outputs of the two modulators are added together to provide an 8-PSK signal. 
The resultant signal constellation is shown in Figure 3. 1.1-3. As can be seen from this 
constellation, as well as from Figure 3. 1.1 -2, the outputs of the two quadrature modulators 
are shifted by 45° with respect to each other. This phase shift permits the two modulators 
to generate eight vectors which are positioned at multiples of 45° along the constellation 
locus circle. These vectors, labeled 0 through 7, are encoded in Gray code as indicated in 
the corresponding parethesis. 

It is important to note that use of two quadrature modulators requires two pairs of 
binary numbers each number consisting of two bits. These bit pairs are labeled a,b and c,d 
for modulators #1 and #2, respectively. The data inputs a,b and c,d are also indicated in 
.the block diagram of Figure 3. 1.1-2. 

However, the maximum number of baseband data streams applied to the modem in 
8-PSK mode is three. Thus, a method for generating four data streams to drive two 
quadrature modulators is required. This method is provided by an octal mapper shown in 
Figure 3. 1.1-4. In that Figure the logic equations which relate the three binary input data 
streams dj(t), d 2 (t), and d 3 (t) to the required binary pairs a,b and c,d are indicated. These 

equations provide for Gray code mapping of the three input data streams by means of two 
quadrature generators. 

Figure 3. 1.1-5 shows logic implementation for the octal mapper. Because of 
relatively high clock rates used by the modem (=25 MHz), it is important that proper timing 
is maintained between signals a, b, c, and d. As can be seen in Figure 3. 1.1-5, signals c 
and d have no logic elements in their paths, while signals a and b have traveled through 
logic gates. Thus, a simple way to compensate for this would be insert similar gates into 
paths dj to c and d 2 to d. These gates should be identical to the type of gates used for 

generating signals a and b. Typically, these may be spare gates in the same IC package as 
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a, b = Quadrature Modulator # 1 Code 
c, d = Quadrature Modulator # 2 Code 


Figure 3. 1.1 -3. Generation of an 8-PSK Signal by Means of Two Quadrature 
Modulators Shifted by 45 Degrees. 



di(t) a(t) = di (t) • d 3 (t) + d 2 (t) - d 3 (t) 


d 2 (t) b(t) = di (t) • d3(t) + d 2 (t)«d 3 (t) 


d 3 (t) c(t) = dj(t) 


d(t) = d 2 (t) 


Figure 3. 1.1-4. Octal Mapper for 3-Channel 8 PSK Modulator Implementation 
with Two Quadrature Modulators. 
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Inputs 


Outputs 



Figure 3. 1.1-5. Octal Mapper Logic Implementation. 
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the "functional” gates. It is also recommended that emitter-coupled logic (ECL) elements 
be considered for modem implementation. 

3.1.2 Configuration for 4/8-PSK with Two AM Modulators 

An alternate approach to using two quadrature modulators for generation of 8-PSK 
signal would be to synthesize the eight phases of the constellation by using the appropriate 
in-phase (I) and quadrature (Q) components of each constellation point. This is illustrated 
in part (a) of Figure 3. 1.2-1. As shown there, signal 0 (Gray code 0, 0, 0) consists of an 
I-component of 0.924 and a Q-component of 0.383 for a normalized RF envelope of unity. 
Part (b) of Figure 3. 1.2-1 gives the I and Q values of all of the sixteen components 
required to form the 8-PSK constellation. Figure 3. 1.2-2 shows the block diagram for 
implementing this approach. As shown in the Figure the Quad/Octal Mapper selects the 
proper digital equivalents of the I and Q components based on the required mapping of the 
input data bits d^t), d 2 (t), and d 3 (t) into the Gray code. 

These digital components are applied to two digital/analog (D/A) converters which, 
in turn, feed the analog I and Q signals to the corresponding amplitude (AM) modulators. 
The AM modulators are driven by RF carriers which are orthogonal to each other and, 
consequently, the summed outputs of these modulators represents the correct phase of 
resultant vector. 

The number of bits (n) used for D/A conversion must be a compromise between 
accuracy and the speed, the latter being determined by the settling time of the D/A 
converters. Our preliminary analysis indicates that an 8-bit accuracy is sufficient for this 
application. Fast settling D/A converters are available with current ECL technology which 
is recommended for this application. 

For the 4-PSK mode, the mapper puts out digital words which result in I and Q 
components of ±1.000. Based on combination of input bits d 3 (t) and d 2 (t), the appropriate 

quadrature coding is generated as was shown in Figure 3. 1.1-1, part (b). 
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(a) 8-PSK Signal Constellation with Gray Code 
Mapping of Coded Binary 3-Tuple 


Octal 

Symbol 

Phase 

(deg.) 

Magnitude of 
Quadrature Components 

I 

Q 

0 

+22.5 

+0.924 

+0.383 

1 

+67.5 

+0.383 

+0.924 

2 

+ 112.5 

-0.383 

+0.924 

3 

+ 157.5 

-0.924 

+0.383 

4 

-157.5 

-0.924 

-0.383 

5 

-112.5 

-0.383 

-0.924 

6 

-67.5 

+0.383 

-0.924 

7 

-22.5 

+0.924 

-0.383 


(b) Signal Definition for 8-PSK 


Figure 3. 1.2-1. Signal Constellation (a) and Coordinate Definition (b) 
for 8-PSK Modulation 
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Figure 3. 1.2-2. Functional Block Diagram of 4/8 PSK Modem Implemented 
with Two AM Modulators in Quadrature 




It is important to note that only one of the AM modulators is actuated at any one 
time when the 4-PSK mode is being used and each of the AM modulators is putting out a 
carrier of unity magnitude and of phase either 0° or 1 80°. Thus, in designing the balanced 
modulators attention must be paid to switching speeds and other factors which determine 
the transient behavior of the modulators. An alternate approach to implementing the 4-PSK 
mode with two AM modulators is generate a set of I and Q components corresponding to 
90° constellation spacing (i.e., I, Q = ±0.707) and use two modulators simultaneously, 
with only the RF phase changing to generate the required QPSK vectors. There may be 
some implementation advantage with this alternate approach. 

At this point, however, we believe that the details of the implementation and 
hardware trade-offs are better left to the modem designers. Also, some breadboard tests 
. may be required to experimentally determine the finer points of the design trade-offs. 
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3.2 Demodulator 

3.2.1 Block Diagram Description 

The block diagram for a 4/8-PSK demodulator is shown in Figure 3.2. 1-1. As 
shown in the figure, there are four phase detectors. Phase detectors PD-A and PD-B 
comprise one orthogonal pair. Detectors PD-C and PD-D comprise the second orthogonal 
pair. Both detector pairs are driven by the local reference carrier VCO, but the VCO signal 
to the second pair (PD-C and PD-D) is shifted by 45 degrees with respect to the first pair 
(PD-A and PD-B). This relationship is shown in Figure 3.2. 1-2. 

It must be noted that for demodulation of a 4-PSK signal only phase detectors A 
and B are sufficient. However, for an 8-PSK demodulation scheme described here, all 
four phase detectors are required to develop the carrier tracking error and to recover the 
data. Thus, in all subsequent discussions, the 8-PSK mode will be emphasized due to the 
fact that the 4-PSK mode is only a smaller subset of the 8-PSK mode. 

As indicated in Figure 3.2. 1-1, the baseband outputs of the phase detectors A 
through D are applied to a unit labled Data Demodulation and Carrier Error Detector 
Circuitry. In this unit the outputs of the phase detectors are processed to generate the 
carrier phase tracking error 0. Also, the output data streams d^t), d 2 (t), and d 3 (t) are 

recovered and applied to their respective destinations. 

The bandwidth of the lowpass filters (LPFs) must be optimized to the incoming 
data rate. Thus, a bandwidth selection feature is also indicated in the block diagram. For 
example, for telemetry only mode the "low" rate is about 100 kbps. On the other hand, for 
digital TV and telemetry, the "high" data rate is in the order of 22 Mbps to 25 Mbps. The 
requirement for the modem to accommodate the high data rate is the primary driver which 
determines the implementation tradeoffs at the circuit design level. 
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Figure 3.2. 1-1. Block Diagram for an 4/8-PSK Demodulator 






(dj , dj , dj ) 


0, 1,2, 7 are octal symbols for phase cod 


Figure 3.2. 1-2. Phase Relationship 
Pairs A/B and C/D a 


2 \ 



3.2.2 Carrier Tracking Error Circuit 

Figure 3.2.2- 1 depicts a block diagram of the carrier tracking error circuit. The 
circuit consists of four baseband multipliers, BM-A through BM-D, four hard limiter 
circuits and several analog adders/subtractors. Similar to the demodulator, only the upper 
portion of the circuit is required for the 4-PSK mode. However, since the uniqueness of 
this modem is defined by its 8-PSK capability, we will concentrate on that aspect of the 
circuit. 

The specific requirement of the carrier tracking circuit, whether it is in the 4-PSK or 
in the 8-PSK modes, is to develop some sort of an "S-curve" which will provide a stable 
lock point, i.e., a point at which the carrier tracking error 0 goes to zero. Providing such a 
curve for suppressed carrier, data-modulated systems is complicated by the fact that the 
carrier tracking signal, in addition to being proportional to the carrier, tracking error also 
contains data modulation. This modulation has to be removed from the carrier error signal 
to make it useable for its intended function. For the antipodal biphase modulation the "third 
multiplier" in the Costas loop provides the function of removing the data modulation from 
the carrier tracking error developed in the "Q-channel" of the loop. As the modulation 
complexity expands to 4-PSK and 8-PSK, the implementation of the carrier tracking loop 
also increases in complexity. 

Using the nomenclature of Figure 3.2.2- 1, the carier tracking error e 0 for the 8- 
PSK mode is: 

e 0 = -A sgn B + B sgn A - C sgn D + D sgn C 
= -Vs cos[0 + <j> m (t)] sgn Vs sin [6 + <f> m (t)] 

+VS sin[0 + <{> m (t)] sgn Vs cos[0 + <J> m (t)] (1) 

-VS cos[0 + 4> m (t) - 45°] sgn Vs sin[0 + <J> m (t) - 45°] 

+V§ sin[0 + <J> m (t) - 45°] sgn Vs cos[0 + <t» m (t) - 45°] 
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BM-A 

= Vscos[e + <i> m (t)] /C>, A sgn B 



al 
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Figure 3.2.2- 1 . Carrier Tracking Error Circuit 



Part (a) of Figure 32 . 2-2 shows the plot of the carrier error voltage e e versus carrier 

tracking error in degrees for an 8-PSK demodulator. The plot shown there is the plot of 
Equation (1) above. Note that this plot is independent of the phase modulation function 
<|> m (t) as long as <J> m (t) is defined as (2k + l)7t/8 where k=0,l,2,...,7 for the 8-PSK 

modulation. It must be also noted that the S-curve "tracks" linearly only within an interval 
defined by ±22.5° which is the spacing between the phases of an 8-PSK signal (see Figure 
3.2. 1-2). 

Part (b) of Figure 3. 2.2-2 shows an S-curve for the 4/8-PSK demodulator in the 4- 
PSK mode.* As can be seen from this curve, the stable lock points for the 4-PSK mode 
are spaced by 90° and the linear tracking range extends over the range of ±45°. This is 
consistent with the QPSK modulation. 

3.2.3 Data Recovery Circuitry 

3.2.3.1 Data Recovery for 8-PSK Mode 

The function of the data recovery circuitry in the 8-PSK mode is to convert the 
four outputs (A, B, C, and D) of the four respective phase detectors (PD-A, PD-B, PD-C, 
and PD-D) into three independent data streams defined by djft), d 2 (t), and d 3 (t). Part (a) 

of Figure 3.2.3- 1 provides a functional representation of this requirement. This "four into 
three" conversion is an inverse operation of the octal mapper shown in Figure 3. 1.1 -4. 

To provide for better physical understanding of the demodulator action in the 
8-PSK mode we have identified axes A and B as I and Q, and axes C and D as I' and Q', 
respectively. Based on the polarity of the I, Q, I', and Q' signals, the bits dj, d 2 , and d 3 is 
0 if Q-output is positive. Conversely, dj = 1 if Q = output is negative. Similarly, from 
Figure 3.2. 1-2, we note that d 2 =1 when I-output is negative and d 2 = 0 when it is 
positive. 


*For derivation of the 4-PSK carrier phase error term, see Appendix A. 
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Figure 3. 2.2-2. Carrier Tracking Error (a) for 8-PSK and (b) for 4-PSK Modulation. 
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+ d,(t) 


Q(B) 
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d 2 (t) 

d 3 (t) 


(a) Functional diagram of data recovery circuit. 
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(b) Data recovery algorithm for an 8-PSK demodulator. 


Figure 3.2.3- 1. Data Recovery Circuitry (a) and Decoding Algorithm (b) 
for an 8-PSK Demodulator. 
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The decoding of the last bit, d 3 , is somewhat more complicated. Observation of 
the phase diagram in Figure 3.2. 1-2 indicates that d 3 = 1 for two cases. One is when both 
I' and Q’ are positive, thus representing octal numbers 1 and 2. The second case for d 3 = 1 

is when both I' and Q' are both negative, i.e., octal numbers 5 and 6 are received. The 
condition d 3 = 0 occurs when I' and Q' have different polarities. This occurs for octal 

symbols 3, 4, 0, and 7. Table 3.2.3- 1 supports the above-presented arguments in terms of 
the actual magnitudes and polarities of the outputs of an 8-PSK demodulator. 

Figure 3. 2.3-2 shows an 8-PSK data demodulation circuit which implements 
the data recovery algorithm given in part (b) of Figure 3.2.3-1. As shown in Figure 3.2.3- 
2, the I/Q and I’/Q* signals are applied to four threshold circuits. The threshold circuits 
convert the analog signals into the digital (binary 1 or 0) formats according to the decoding 
algorithm chosen. The operation of convertin the four baseband signal outputs of an 8- 
PSK demodulator into three data streams, i.e., dj, d 2 , and d 3 , is thus accomplished. 

3.2.3.2 Data Recovery for the 4-PSK Mode 

For the 4-PSK mode, i.e., a conventional QPSK modulation, the demodulator is 
reconfigured to use only the upper half of the circuitry shown in Figures 3.2. 1-1 and 
3.2.2- 1. The details of the demodulator operation in the QPSK mode is given in 
Appendix A. A different labeling is used in the Appendix, but functionally a(t) in the 
Appendix is d t and b(t) is d 2 . 
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Output Axis* 

Octal 

Symbol 

Modulation 
Angle M 
Degrees 

I = cos M 

Q = sin M 

I'=cos(M - 45°) 

Q'=sin(M -45°) 

0 

22.5 

0.924 

0.383 

0.924 

-0.383 

1 

67.5 

0.383 

0.924 

0.924 

0.383 

2 

112.5 

- 0.383 

0.924 

0.383 

0.924 

3 

157.5 

- 0.924 

0.383 

- 0.383 

0.924 

4 

202.5 

- 0.924 

- 0.383 

- 0.924 

0.383 

5 

247.5 

-0.383 

- 0.924 

- 0.924 

- 0.383 

6 

292.5 

0.383 

- 0.924 

-0.383 

- 0.924 

7 

337.5 

0.924 

- 0.383 

0.383 

- 0.924 


M = (2d -1 )tc/ 8 d = 0,1,2 7 

♦Output signals normalized to a maximum value of unity. 


Table 3. 2.3-1. Demodulator Output Signals for an 8-PSK Demodulator. 
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3.3 Design Considerations 

The addition of an 8-PSK capability to the baseline 4-PSK mode imposes certain 
design considerations not only on the modulator and the demodulator units but also on the 
transmitter and the receiver units of the MA link as well. In this section, we address some 
of these considerations particularly from the standpoint of how they affect the performance 
of the 8-PSK mode. 

3.3.1 Channel BER 

3.3.1. 1 General 

Throughout this report we have been referring either to two data streams, 
identified as dj and d 2 , for the 4-PSK mode or to three data streams, identified as d lt d 2 , 
and d 3 , for the 8-PSK mode. Operationally, this means that we have either a two-channel 
capability (dj and d 2 in the 4-PSK mode) or a three-channel capability (dj, d 2 , and d 3 in 
the 8-PSK mode). Because the on-going discussion concerns the 8-PSK mode, we will 
refer to the independent channels available in this mode as channels 1, 2, and 3. 

Consider first the BER behavior of the 8-PSK mode. Figure 3. 3. 1.1-1 shows 
the BER performance comparison of the 8-PSK mode with that of the baseline 4-PSK and 
the 2-PSK modes which are included in the figure for purpose. An interesting feature of 
the 8-PSK curve is that channels 1 and 2 provide a slightly better BER performance than 
channel 3. This is due to the fact that, as was shown in Figure 3.2.3-1, the recovery of the 
data streams d 2 and d 2 is based on the hard discisions in I and Q half-planes, but the 
decision for the d 3 data stream is based on narrower regions of several quarter-planes. In 
terms of the penalty in the E S /E N ratio, this difference between channels 1, 2 and channel 3 
is only about 0.3 dB at BER of 10 5 , which is the specified BER for the digital video 
transmission. The advantage, however, is the simplicity of the data demodulation circuitry 
as was illustrated in Figure 3.2.3-2. 
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Channel BER 



Figure 3. 3. 1.1-1. BER vs. PSK Modulation Level and Symbol -to-noise 

Ratio E s /Efj- 
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An alternate method of recovering the data streams dj, d 2 , and d 3 would be to 
use only the I and Q axes and solve for the modulation angle M using the I and Q 
components recovered by the I and Q-axis demodulators. This would be an optimum 
symbol detection method, but for the Gray incoding, it would still make channel 3 BER 
somewhat worse than channels 1 and 2. 

We at Axiomatix have also explored a data decision algorithm which equalizes 
the BER in all three channels. However, such an equalization results in a BER which is 
equal to that of channel 3 as shown in Figure 3. 3. 1-1. Thus, the conclusion is that the RF 
link which uses an 8-PSK modulation would be designed to achieve BER of 10‘ 5 in 
channel 3 based on curves in Figure 3.3. 1-1. In this case channels 1 and 2 would simply 
have an additional margin of 0.3 dB. 

3.3.1.2 Effect of Data Angle Imbalance on Channel BER 

As explained in Section 3.2, the data demodulation boundaries are defined by 
two sets of orthogonal axes I/Q and I'/Q’. An examination of the phase relationship shown 
in Figure 3.2. 1-2 indicates that for an ideal implementation all eight phases of the 8-PSK 
signal are never more than 22.5 degrees from either one of the orthogonal axes defined by 
the I/Q and I'/Q' pairs. With a practical implementation of the modulator/demodulator, a 
question arises in regards to the angular (i.e., phase) tolerance, or imbalance, with which 
the 8-PSK data is generated at the transmitter and demodulated at the receiver. Imbalance is 
defined as per Figure 3.3. 1.2-2. Figure 3. 3. 1.2-1 shows the effect of imbalance on the 
BER of the channel 1 , 2, and the channel 3 data outputs. From this figure it is evident that 
at Es/E n of 18 dB, there is very little angular tolerance allowed for imbalance. This, in 

turn, translates into a very "tight" angular tolerance on the angular relationship between the 
I, Q, I’, and Q' axes. For Es/E N = 18.5 dB the allowed tolerance increases, and for Es/E n 
= 19 dB the allowed tolerance falls within the limits (90° ± 2°) of the practical quadrature 
hybrids. Such hybrids would be used for implementing both the modulator and the 
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1GLE 

3.3. 1.2-1. Effect of Imbalance on 

Channel BER vs. E s /E^ Ratio 
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demodulator units. This increase of Ej/En from 18 to 19 dB can be attributed in a link 
budget to what is generally referred to as an "implementation loss" of 1 dB. 

It must be noted in Figure 3.3. 1.2-1 that the equal BER values for channel 1,2 
and channel 3 appear at an angle which is approximately 22.2° rater than the "nominal" 
angle of 22.5°. This is due to the fact that we are considering each channel individually 
rather than treating the average rate of all three channels. If we consider the BER value 
averaged over all of the three channels, we obtain curves shown in Figure 3. 3. 1.2-3. 
From these curves it is evident that the minimum in average BER value does occur at the 
nominal value of 22.5°. 

The significance of this discussion is that for application on hand, i.e., providing 
a capability to transmit three independent channels, one must analyze the performance of 
each channel individually rather than treat the performance of an 8-PSK modem on an 
"average" basis. 

3.3.2 Tolerance to RF Circuitry Impairments 

When utilizing higher order modulation techniques, such as 8-PSK, 16-PSK, 
16 QAM, etc., one must pay attention to the effect of the RF circuitry impairments on the 
performance of the modulation technique. This, of course, is also true for the bi-phase 
(2-PSK) and the quadriphase (4-PSK) modulations. But, because we are introducing an 
8-PSK technique here, our primary concern is the relative sensitivity of the 8-PSK 
technique to such impairments as compared to a 4-PSK modulation. 

Our consideration, therefore, focuses specifically on such impairments as 
(1) amplitude distortion, (2) delay distortion, and (3) carrier phase jitter. Although, at this 
point in time, it is difficult to determine the exact location of these impairments in the RF 
chain extending from the transmitter to the receiver, the purpose of this discussion is to 
make the potential users of the 8-PSK aware of the effect of these impairments. 
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3.3.2.1 Amplitude and Delay 


Figure 3.3.2. 1-1 shows the power penalties due to amplitude and delay 
distortion for several high order modulation schemes. For both the amplitude and the 
delay, there are two terms: (1) a symmetrical term and (2) an asymmetrical term. These 
terms are referred to, respectively, as "linear" and "quadratic" terms. Also, these terms are 
defined by their respective coefficients, i.e., Aj and A 2 for the amplitude response and dj 
and d 2 for the delay response. 

An example of using the information in Figure 3. 3. 2.1-1 may be as follows. 
Suppose we are considering a symbol rate of 22 Msps, the duration of each symbol is 
about 45 nsec. If the attenuation of the bandpass filters proceeding the demodulator is 1 dB 
at ±1 1 MHz (i.e., A 2 = ldB), the power penalty will be about 0.7 dB for the 8-PSK mode 
.and about 0.3 dB for the 4-PSK mode. Similarly, if d x = 9 nsec then the power penalty 

paid for this distortion alone will be about 0.9 dB for the 8-PSK mode and about 0.2 dB 
for the 4-PSK mode. These penalties define an increase in EJE N required to maintain BER 

of 10‘ 6 if each of these distortions is taken separately. 

In an actual system, there will be all four distortions present, each claiming its 
own power penalty. Thus, from the standpoint of the RF system design, the impact of 
using the 8-PSK has to be reflected in the "baseline" specification (i.e., the 4-PSK spec). 
If this is not done, then the 8-PSK mode has to be provided with sufficient amount of 
additional RF power to work over the channels baselined for the 4-PSK. 

3.3.2.2 Carrier Phase Jitter 

The effect of carrier phase jitter on the power penalty is shown in Figure 
3.3.2.2-1. For example, a phase jitter of 2 degrees RMS would result in power penalty of 
about 0.4 dB for the 8-PSK system and about 0.1 dB for the 4-PSK system. 

The information provided in Figure 3.3.2.2-1 is useful for determining whether 
the phase noise characteristics of the local oscillators, synthesizers and other frequency 
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1-1. Power Penalties Due to Amplitude and Delay Distortion 
at BER = 10- 6 (Ref. [3]). 
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CARRIER PHASE JITTER (DEGREES RMS) 


Figure 3.3.2.2-1. Power Penalties Due to Carrier Phase Jitter at BER = 10 -6 
(Ref. [3]). 
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sources used by the Ku-band MA system have a significant impact on the power required 
to establish an 8-MSK link. Conversely, if the phase noise level is established by factors 
other than the 8-MSK operation, then the designers of the 8-PSK link must consider means 
within the framework of the modulator/demodulator design to reduce the effect of the phase 
noise on link performance. One parameter, for example, which determines the carrier 
phase jitter is the bandwidth of the carrier tracking loop for the 8-MSK demodulator. For 
an MSC application the doppler shift may not be a problem because of small relative motion 
between the transmitter and the receiver units. Thus, the trade off is then between making 
the carrier tracking loop wide enough to maintain lock despite the carrier phase noise or to 
make it narrow enough to minimize the effect of the thermal noise. 

It must be also stated here that the power penalties indicated for the impairments 
are for BER of 10' 6 . For BER of 10' 5 , which is a baseline requirement for the digital TV 
links, these penalties will be less. 

3.3.3 Implementation Philosophy 

Digital implementation is recommended, wherever applicable, for the design of both 
the modulator and the demodulator units. Digital implementation of the modulator unit was 
discussed briefly in Section 3.1.2. For the demodulator design, the digital implementation 
will particularly be of great benefit. Figure 3.3.3- 1 shows a block diagram of a digitally 
implemented demodulator. As shown, the baseband signals from phase detectors A, B, C, 
and D are converted to digital format and applied to a digital processor. 

The processor performs the functions of (1) carrier tracking, (2) data demodulation, 
and (3) data clock recovery. For the low data rate (= 100 kps), several microprocessors 
may be used to achieve to desired processing speed. However, for the video data rates 
(=25 Mbps) special purpose, high speed processor circuitry has to be used. 
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Figure 3.3.3- 1. Block Diagram of a Digitally Implemented Demodulator. 









4.0 CONCLUSIONS 

It has been shown that the baseline 4-PSK capability of the Space Station's Ku- 
band MA system can be expanded for some users to provide three simultaneous digital TV 
channels by using 8-PSK modulation. The encorporation of the 8-PSK capability can be 
provided by designing what is referred to in this report as a "multi-channel" 
modulator/demodulator. The design philosophy is to provide the 8-PSK capability as a 
selectable "add-on" capability to the baseline 4-PSK mode. In this manner only those users 
who require the three channel capability will have to carry the additional circuitry. This 
circuitry may be either on a "plug in" card or in some other form of a plug-in module. 

This report presented a conceptual design for a 4/8-PSK modulator/demodulator 
unit. The design concept described here can be used as a starting point for a more detailed 
engineering design. The selection of an IF frequency still remains an open issue. But, the 
impact of the IF frequency value will be primarily on the circuit design of the 
modulator/demodulator rather than on the conceptual design. 

The subject of the clock recovery is not addressed in this report. We feel that it is a 
separate topic. However, there are techniques reported in literature which provide for a 
relatively easily implementable means of recovering the clock information from an 8-PSK 
signal [4]. For example, the technique described in [4] uses narrow-band filtering of the 
8-PSK IF signal to recover the clock component which is present in the RF envelope of an 
8-PSK modulated signal. 
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DEMODULATOR FUNCTION IN THE 4-PSK MODE 
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DEMODULATOR FUNCTION IN 4-PSK MODE 


A-l INTRODUCTION 

This appendix presents a simple analysis of the demodulator operation in the 4-PSK 
mode, which is generally referred to as the QPSK mode. 

Figure A-l shows the block of a QPSK demodulator. As shown, the demodulator 
consists of two orthogonal phase detectors PD-A and PD-B. For carrier error generation, 
the outputs of these detectors are cross-coupled via two hard limiters and baseband 
multipliers BM-A and BM-B. 

An analysis of the method of the carrier error generation is given below. 
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Figure A-l. Demodulator for 4-PSK (QPSK) Mode. 



A-2 ANALYSIS 

To simplify analysis, we will represent the input signal to the QPSK demodulator in 
terms of two in-phase and quadrature components. These components are labeled a(t) and 
b(t), respectively. Figure A-2 shows the idea. The input signal to the demodulator can, 
therefore, be written as: 

e in = VS [a(t) cos co 0 t + b(t) sin co 0 t ] , (1) 

where a(t) = ±1 and b(t) = ±1, depending on the data patterns of the respective data 
streams. The output of phase detector PD-A is 

e PDA = Vs [a(t) cos co 0 t + b(t) sin co 0 t] V2 cos(co 0 t + 6) 

=Vs/2 [a(t) cos 0 - b(t) sin 0] . (2a) 

Similarly, the output of phase detector PD-B can be obtained 

e PD B = Vs [a(t) cos fo 0 t + b(t) sin co 0 t] V2 sin(ca 0 t + 0) 

=Vs/2 [a(t) sin 0 - b(t) cos 0] . (2b) 

As shown in Figure A-l, the outputs of the phase detectors are, after low pass 
filtering, applied to baseband multipliers BM-A and BM-B. The outputs also are applied to 
two threshold detectors which not only determine the polarity of the appropriate signals, 
but also "clean up" the signals a(t) and b(t) making them look more or less like square 
waves. Thus, for a resonably high SNR in the LPFs, we can assume that 

sgna(t) = a(t) , 
sgn b(t) = b(t) . 

Using the relationship in (3) above, we proceed to find the expressions for the outputs of 
the baseband multipliers BM-A and BM-B: 
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e BMA = ^S/2 [a(t) cos 0 - b(t) sin 0] b(t) 

= Vs/2 [a(t) b(t) cos 0 - b(t) b(t) sin 0] (4a) 

= VS/2 [a(t) b(t) cos 0 - sin 0] 
and 

e BMB = ^S/2 [a(t) sin 0 + b(t) cos 0] a(t) 

= Vs/2 [a(t) a(t) cos 0 + a(t) b(t) cos 0] (4b) 

= VS/2 [sin 0 + a(t) b(t) cos 0] . 

It must be noted that in (4a) and (4b), we have made a valid assumption that a(t) a(t) = 1 
and b(t) b(t) = 1. For the next step, we will make a further resonable assumption, specially 
that 

a(t) b(t) = a(t) b(t) . (5) 

Using this assumption we will show that the carrier tracking error out of the summer is 


e 0 ~ e BMB ~ e BMA 
= 'FIs sin 0 , 


( 6 ) 


which for 0 = small becomes 

e e = "FlS 0 . 

The tracking capability of the QPSK demodulator extends only to ±45° because, as can be 
seen from Figure A-2, when 0 > 45° the output of PD-A becomes b(t) and the output of 
PD-B becomes a(t). Thus, the data streams "flip." 
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Figure A-2. Phase Diagram for DPSK Demodulation. 
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1.0 INTRODUCTION 

The Space Station Information System (SSIS) is intended to provide reliable 
information flow between SSIS users. In order to provide the quality of service required, 
forward error correction (FEC) and forward error detection (FED) codes will be used. 
Axiomatix was tasked to examine several candidate codes, and propose alternatives as 
deemed necessary. In this memo, we discuss the performance of the candidate block codes 
for the SSIS transfer frame and transfer frame header. Performance is determined by 
analysis or simulation, as appropriate. Details of the performance assessment are covered 
in this memo. 

1.1 Transfer Frame Requirements 

The SSIS transfer frame will consist of a header field and a data field. An error 
detection code may be included in the header field to provide the required header reliability, 
which we understand to be on the order of 10 -12 bit error rate, as per [Ref 1). In addition, 
an error correcting code will be added to the transfer frame to mitigate the channel errors 
and improve message throughput. The required output bit error rate of the decoded transfer 
frame depends on the grade of service: grade 1 is 1(T 12 , grade 2 is 10 -8 and grade 3 is 
10~ 5 . However, a more relevant measure of message fidelity is the overall probability of 
transfer frame error or header error. In terms of random bit error rate, the probability of 
correct message is the probability that all the constituent bits in the message are correct, and 
the probability of message error is given as: 

P ME = i-0-ber) n , 

with BER the bit error rate, P ME the message (codeword) error rate, and N the length of the 

j 2 

coded message, in bits. Thus, for a 64 bit transfer frame header, a bit error rate of 10 


corresponds to a header error rate of 6.4 x 10 -11 , and for a 2040 bit transfer frame 
codeword, a bit error rate of 10 -12 corresponds to a frame codeword error rate of 2 x 10 9 . 
If the transfer frame is encoded using a Reed-Solomon code, as proposed in [Ref 2], and 
error correction is performed, the resultant output will either contain no errors, or will 
contain several bursts of errors since the received sequence will be transformed to a 
codeword differing in at least 33 symbols of 8 bits each. The point behind the preceding 
discussion in that average bit error rate out is not the appropriate statistic to specify, but 
rather the average frame error rate should be considered. In the following discussion, the 
error detection and correction performance of the codes will be presented in terms of frame 
or frame header error rate, as the transformation to bit error rate is somewhat of an artifact. 

1.2 General Code Consideration 

We have examined the performance of the two codes proposed in [Ref 2], The 
performance of the (255, 223) Reed-Solomon code proposed for transfer frame error 
correction was investigated using analytical techniques. The codeword weight structure of 
Reed-Solomon codes is well known, and the message error rate and throughput can be 
calculated explicitly, based on the channel symbol error rate. At this point, we need to 
make a distinction between channel symbol error rate and channel bit error rate. The 
weight structure of the Reed-Solomon codes is based on symbols, with each non-zero 
symbol consisting of one or more non-zero bits. Thus, the weights given by the formula 
of Section 2 are independent of the number of non-zero bits, so the the weight in bits 
cannot be derived from this result. As a consequence, the probability of undetected error as 
a function of bit error rate, rather than symbol error rate, cannot be derived in closed form 
from the symbol weight structure. Even though symbols for the code can be reconstructed 
from bits transmitted over a binary channel, the codeword errors induced by bit errors have 
a different probability measure than the same errors induced by symbol errors. That is, in 
an M-ary symbol channel, with n = log 2 M bits , a received sequence with w errors has 



in a code of length N. However, 


maximum probability of ^1 - ^ — j in a code of length N. Howevei 

the same codeword, if it were to consist of symbols with one bit error each, would have a 


maximum probability of j^l - Jjj W . The latter probability is considerably larger 

than the former, and the potential for undetected errors using a binary channel could be 
higher than with a non-binary channel using the same Reed-Solomon code. Simulation 
with a shorter code, discussed in Section 3, indicates that there are not enough codewords 
with low binary weight symbols to cause this to be a problem for the Reed-Solomon code 
tested. The (255,223) code proposed for the transfer frame has enough distance, even after 
correcting sixteen errors, to provide a probability of undetected error less than 10" using 
a non-binary channel. However, simulation will probably be needed to verify that this 
code will perform as well with a binary channel. In fact, a rather lengthy simulation will be 
required to verify a 10 -13 error rate, but it will be worthwhile if this one code can provide 
both error correction and error detection with the required probability of false message. 
The simulation may be expedited by testing codewords with a low binary weight. We 
discuss this issue in more detail in Section 2. 

The frame header code candidates are examined using a simulation, thereby making 
it possible to verify the code performance with a binary channel. Since the error rates 
involved are higher, this simulation is possible on a personal computer. We have run an 
extensive simulation involving three potential header codes, a Hamming code, a binary 
BCH code, and a Reed-Solomon Code. All three codes are 64 bits long with 16 check 
bits, which gives an undetected error rate bound of 2 -16 , or 1.5 x 10 -5 at an input error rate 
of 0.5. Results of this simulation are discussed in Section 3. 


2.0 TRANSFER FRAME ERROR CORRECTION CODE 

We have evaluated the performance of the (255,223) 8-bit per symbol Reed- 
Solomon code proposed by Miller and Morakis [Ref 2]. Results from this paper, 
reproduced as Figure 2.0-1, show a probability of undetected bit error exceeding 10' at 
low signal-to-noise ratio. We believe that these results are based on a weak upper bound, 
or based on a maximum likelihood decoding algorithm, as our results indicate a probability 
of undetected error of less than 10' 13 at all signal-to-noise ratios for a non-binary channel, 
and we will argue that the same result may apply to a binary channel. Our results assume 
that we correct a maximum of 16 errors, the guaranteed capability of the code. 

We examine the error detection performance of the code using two techniques. 
First, we calculate the total number of codewords in the received signal space. Then, 
assuming we are going to correct a maximum of 16 errors, we calculate the number of 
codeword neighbors within distance 16. An erroneous codeword will be accepted if more 
than 16 errors occur, and the received sequence lands within 16 symbols of another 
codeword, which we call a neighbor. Thus, under the condition that landing on a 
codeword or a neighbor is no more or less likely than landing on any other possible 
received sequence at the same distance, the probability of error is upper bounded by the 
ratio of codewords and nearest neighbors to the total possible number of received 
sequences. 

We can calculate the number of codewords in an (N,k) code over GF(q) as q k , 
which in this case is 25 6 223 . The number of sequences within distance 16 of any given 
codeword is equal to the number of sequences of weight less than 17, which is 

X !)* = -302 x 10 64 

for the code in question. The total number of sequences is q = 256 . Thus, the ratio of 

codewords and neighbors to the total calculates out to 0.26 x 10' 13 . This is the worst case 
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probability of landing within 16 symbols of another codeword. 

A more rigorous expression which relates undetected codeword error to channel 
error rate can be derived based on a formula from Forney [Ref 3]. The number N d w of 
codewords of weight w for a code of minimum distance d is 
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and the probability of error, i.e., the probability of receiving a sequence within distance 16 
of another codeword is given as 


N 16 


.-XU X 


(N-w) ! w ! (q-1 ) k ~ w ~ i ~ g, (q-2) 1 P 


N« 


e- l "» 
w=d k=0 i=0 £=0 


\i r>w+2jl+i— k 


q-p) 


N-w-2i,-i+k 


!(N-w-£)!i!(k-i -l)\ (w-k+£)! 


with m = min (k-i, N-w), and p the channel symbol error rate. This expression is 
corrected from [Ref 3], which has an error in the upper limit of the inner sum. This 
expression was evaluated using double precision Fortran, and the results are plotted in 
Figure 2.0-2. Note that, as expected, the maximum undetected error rate is equal to 
0.26 x 10~ 13 . 

The use of a decoder that corrects a fixed maximum number of errors, 16 in this 
case, will result in rejection of transfer frames with more than 16 errors. The probability of 
accepting a correct message, which is the probability of 16 or fewer errors, is a function of 
the channel error rate. The correct message rate is plotted in Figure 2.0-3. The use of a 
decoder that corrects to some codeword regardless of the number of errors will result in a 
slighdy higher probability of accepting a correct message since there are some error patterns 
of 17 or greater that will still be closer to the true codeword than any other codeword. 
However, the probability of undetected error goes to one as the channel error rate 
deteriorates since all detection capability is lost. Since the frame header code does not 
provide the required error rate, as we will show in the next section, it would be advisable to 
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use algebraic decoding with a maximum error correction of 16 symbols rather than 
maximum likelihood decoding with arbitrary error correction. 

21 Reed-Solomon Code Performance over a Binary Channel 

The previous results cover Reed-Solomon code performance over a non-binary 
channel, or a channel with equi-probable symbol errors, as depicted in Figure 2.0-2. The 
worst case error detection performance using this channel can be heuristically bounded by 
noting that, regardless of the error pattern affecting the k information symbols, the errors 
affecting the N-k check symbols must cause perfect agreement between every computed 
check symbol and every received check symbol. Note that with a Reed-Solomon code, any 
change in one information symbol affects each and every computed check symbol. The 
worst case probability of receiving a check symbol that agrees with an erroneous computed 
check symbol is P e / (q-1), which is less than (q-1) -1 . The probability that all check 
symbols agree is (q-l)" (N_k) with error detection only, and no error correction. With 
q = 256 and N-k = 32, this gives a P e < 10~ 77 . 

By comparison, with a binary channel, low binary weight error patterns are more 
likely than error patterns with a large number of bit changes. The concern is that a 
codeword differing from the transmitted codeword by N-k+1 symbols may only differ by 
N-k+1 bits, one bit per symbol. In a manner analogues to the non-binary channel, we can 
calculate the probability of match for a set of N-k check symbols consisting of n bits each, 
and each with one bit error. The probability of one check symbol matching is pq n_1 with 

one bit error in the symbol. This is maximized with p = ^, so that the probability of N-k 

lYn-l)(N-k) 

check symbols matching with one error each is P c ^-J [} ~ nj * t0T 

our example, this is P c < (JJ 2 • ^JJ* 32 = 1CT 42 for one specific sequence of check 

symbols given that the appropriate information symbols have been received. The overall 
probability of error is much less since the overall probability of receiving the information 


symbols as well as the check symbols must be factored in. The question that we can't 
answer at this point is how many low binary weight codewords exist relative to the total 
number of codewords. In all likelihood, the density of binary weight check sequences is 
uniform for a given binary weight of the information sequence, in which case we can argue 
that the error detection performance over a binary channel is equivalent to the performance 
over a non-binary channel. Recall that a Reed-Solomon code has the property that a change 
in any one information symbol changes each and every check symbol. Thus, each check 
symbol can be written as a linear combination of the information symbols as shown below 

C 0 = a 0>0 1 0 + «o,i h + * * * + a o,k-i ! k-i 
Cl = <Xi,o !o + a i,i + * * * + a i,k-i h-i 

C31 = «31,0 lo + a 31,l + •** + a 31,k-l I k-l 

with all the coefficients a y non-zero. Now, the next step is to argue that is a column of 
the coefficient matrix of the a's is not a constant, each change of an information symbol 
will result in every check symbol changing by a different amount corresponding to a 
different value of the coefficient. For example, if the all-zero codeword is transmitted and 
one information symbol is received in error, all the check symbols will be non-zero. At 
most eight can have binary weight one if all the coefficients in one column are unique, and 
the remaining check symbols have minimum binary weight two. Thus, under this 
assumption, the minimum binary weight non-zero codeword is 1 + 8 + 2 x 24 = 57. In 
order to pursue this, the coefficient matrix for a specific code needs to be generated, and the 
properties examined. Looking at the code from this perspective, however, intimates that a 
change in one or more information symbols "randomizes" the check symbol sequence. 


To lend some credence to these arguments, we examine the performance of a Reed- 
Solomon code in the next section, and compare it to a binary BCH code and a Hamming 
code of the same rate. The simulation results will show that the Reed-Solomon code 
outperforms the other two. While we can't run a full-up simulation of the (255,223) code 
at all error rates, we can examine areas of concern. That is, we can induce specific 
numbers of errors, 17, 18, and so on, to bias the simulation to a high error rate and 
accelerate the statistical reliability of the results. 

We feel that this effort will be justified since this code has the potential for 
simultaneously satisfying both the error correction desired and the false message 
requirements. 

3.0 FRAME HEADER SIMULATION RESULTS 

A frame header code to provide protection against acceptance of erroneous headers 
has been prepared. [Ref 1] specifies a header BER of 10~ 12 This bit error rate is not 
possible at all input error rates using the candidate error detection code proposed. The code 
proposed in [Ref 2] is a (64,48) binary Hamming code. In fact, the best any code with 16 
check bits can do is guarantee an error rate of 2 -16 unless the input, or channel, error rate is 
small. Fortunately, the results of section 2 indicate that this BER requirement can be met 
by the channel code. However, should a header code be required, we have simulated the 
performance of the proposed Hamming code as well as a comparable binary BCH code and 
a Reed-Solomon code. A1 codes were tested with the same binary error pattern to provide 
consistency. 

The Hamming code has a generator polynominal given as 
g(x) = 1 + x 5 + X 12 + X 16 . 

The BCH code had a generator polynominal given as 

g(x) = (1 + X 2 + x 3 + x 4 + x 8 ) (1 + X + x 2 + x 4 + x 5 + x 6 + x 8 ) 

= 1 + X + x 5 + x 6 + x 8 + x 9 + x 10 + x 11 + x 13 + x 14 + x 16 


The Reed-Solomon code check symbols consist of the first three symbols, C 0 , Cj, and C 2 , 
generated by constraining the first three syndrome equations to zero. The fourth symbol is 
an overall parity check on the codeword. Header error rates for these three codewords are 
plotted in Figure 3.0-1. Ultimately, the Reed-Solomon code outperforms the other two. 

A listing of the frame header simulation follows. A copy of this program was given 
to LOCKHEED personnel at JSC for their use. It compiles under LAHEY F77L 
FORTRAN using a PC with the 8087 coprocessor. 


[Ref 1] "Space Station Specification," Section 4.3.2.2, January 5, 1987 

[Ref 2] "Channel Coding Performance Analysis," (Private Communication), Warner 
Miller and James Morakis, February, 1987. 

"Concatenated Codes," Forney, G. David, The MIT Press, 1966. 


[Ref 3] 
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PROGRAM RS 

INTEGER*2 IC ( 0 : 1 5 ) , IFR , 
DIMENSION IETAB( 0 : 31 ) 


IMUL(0: 15) ,ISUM, IETAR(0:3) 


Hamming generator 

DATA IFH/Z' 10210000' / 

Reed-Solomon symbol field generator polynomial 
DATA IFR/Z ' 13 ' / 


* BCH generator 

DATA IFB/Z'6F630000'/ 

* 8087 Initialization 

CALL INVAL0P( 12 ) 

0PEN(2,FILE="LPT1") 

* Y=RRAND( ) 

I0NE= 1 
CWCNT=0 . 

CWEH=0 

CWEB=0 

CWER=0 

WRITE(*,*) "INPUT BER" 

READ( * , *) PE 

* Initialize PR generator 

CALL DRANSET 
IETAB(0)=1 
DO 10 1=1,31 

* 8087 word shift 

* IETAB and IETAR are used to add errors to input 

IETAB( I )=ISHFT(IETAB( 1-1 ) , IONE) 

10 CONTINUE 

IETAR(0)=1 

DO 101 1=1,3 

IETAR(I)=2*IETAR(I-1) 

101 CONTINUE 

* IMUL is GF(16) multiplication table 

IMUL(0)=0 
IMUL ( 1 ) =2 
DO 15 1=2,15 
IMUL( I )=I*2 

* IEOR is 8087 XOR 

IF (IMUL(I) .GE.16) IMUL( I )=IE0R( IMUL ( I ) , IFR ) 

15 CONTINUE 

* Print BER on screen @ 10K, on printer @ 100k cwords 

30 DO 21 1=1,10 

DO 20 IPR-1, 10000 

* Clear R-S input words 

DO 33 J=0 , 1 5 
33 IC(J)=0 

* Input error counter 

BTMP=0 

* Clear BCH input words— Note B is BCH, H is Hamming 

IR1 B=0 
IR0B=0 

* Clear Hamming input words 

IR1 H=0 
IR0H=0 

* Clear decoder syndrome registers (Ref. pg. 150 Vol 1 Peterson) 

ICDH=0 

ICDB=0 

11=0 

* Generate first 32 input bits packed into IC and IR1B 

DO 50 12=0,7 
DO 60 13=0,3 


* Generate PR# 

Y=DRAN( 1 . ) 

* Y=RND ( ) 

IF(Y.LE.PE) THEN. 

* Set bit to 1 if PR# less than input error rate 

IC(I2) = IOR(IC(I2) , IETAR ( 13 ) ) 
BTMP=BTMP+1 

IR1B=I0R(IR1B, IETAB(Il) ) 

ENDIF 

11=11+1 

60 CONTINUE 

50 CONTINUE 

11=0 

* Generate second 32 input bits 

DO 501 12=8,15 
DO 601 13=0,3 
Y=DRAN ( 1 . ) 

IF(Y.LE.PE) THEN 

IC(I2)=I0R(IC(I2) , IETAR (13)) 
BTMP=BTMP+1 

IR0B=I0R(IR0B,IETAB(I1)) 

ENDIF 

11=11+1 

601 CONTINUE 

501 CONTINUE 

* Copy BCH input to Hamming input 

IR1H-IR1B 

IR0H=IR0B 

♦.Increment cword counter 
CWCNT=CWCNT+1 . 

* Skip if no input errors 

IF(BTMP.EQ.O) GOTO 20 
BERR=BERR+BTMP 
DO 110 16=0,63 

* Sign bits hold current bit to be processed 

* XOR syndrome , input for BCH 

IFLGB=IE0R(ICDB,IR1B) 

ICDB=ISHFT(ICDB,IONE) 

* 64-bit left shift 

IR1B=ISHFT(IR1B,I0NE) 

IF(IROB.LT.O) IR1B=I0R(IR1B,I0NE) 
IR0B=ISHFT(IR0B,I0NE) 

* Feedback g(x) if high order bit is 1 

* Next do Hamming 

IF(IFLGB.LT.O) ICDB=IE0R( ICDB , IFB) 
IFLGH=IEOR( ICDH.IR1H) 

ICDH=ISHFT( ICDH, IONE) 

IR1H=ISHFT(IR1H, IONE) 

IF(IROH.LT.O) IR1H=I0R(IR1H,I0NE) 
IR0H=ISHFT(IR0H,I0NE) 

IF(IFLGH.LT.O) ICDH=IE0R( ICDH , IFH) 

110 CONTINUE 

* Increment error count if syndrome not zero 

IF (ICDB.EQ.O) THEN 
CWEB=CWEB+1 

ENDIF 

IF (ICDH.EQ.O) THEN 
CWEH=CWEH+1 

ENDIF 

* R-S decoder 

1004 ISUM=0 


* First do overall parity check 

DO 80 15=0,15 

80 ISUM=IEOR(ISUM, IC(I5) ) 

* Detect error if nonzero, i.e. not UNDETECTED error 

IF(ISUM.NE.O) GOTO 20 

* Calc first syndr orae=sum( IsubI4*Alpha**I4 ) 

ISUM=IC(0) 

DO 70 14=1,14 
ISUM=IMUL(ISUM) 

ISUM=IE0R( ISUM,IC(I4)) 

70 CONTINUE 

* Detect error if nonzero 

IF(ISUM.NE.O) GOTO 20 
ISUM=IC(0) 

* Second syndr ome=sura( IsubI6*Alpha**(2*I6) ) 

DO 90 16=1,14 
ISUM=IMUL(ISUM) 

ISUM=IMUL(ISUM) 

ISUM=IE0R(ISUM,IC(I6)) 

90 CONTINUE 

* Detect error if nonzero 

IF(ISUM.NE.O) GOTO 20 
ISUM=IC(0) 

* Third syndrome 

DO 100 17=1,14 
ISUM=IMUL( I SUM) 

ISUM=IMUL( ISUM) 

ISUM=IMUL(ISUM) 

ISUM=IE0R( ISUM , IC( 17 ) ) 

100 CONTINUE 

* If all syndromes zero, we have undetected error, inc 

IF (ISUM.EQ.O) THEN 
CWER=CWER+1 

ENDIF 

20 CONTINUE 

* Calc error rates, PB is check on input rate 

CWOUTB=CWEB/CWCNT 
CWOUTH=CWEH/CWCNT 
PB=BERR/CWCNT/64 
PUniTTR-PWFR /rwrNT 

WRITE ( * , * ) PB , CWCNT , CWOUTH , CWOUTB , CWOUTR 

21 CONTINUE 

WRITE ( 2 , * ) PB , CWCNT , CWOUTH , CWOUTB , CWOUTR 

GOTO 30 

END 

* Big mother PR# generator 

FUNCTION DRAN(X) 

DOUBLE PRECISION M , A , C , RM , Y , XX 
DIMENSION RANT(97 ) 

J=1 .DO+(97.DO*Y)/M 
Y=RANT( J) 

DRAN=Y*RM 
XX=DM0D( A*XX+C , M) 

RANT( J)=XX 
RETURN 

* Initialization for PR# gen. 

ENTRY DRANSET 

DRANSET=1 

M=1771875. 

A=24 1 6 . 

C=374441 . 


count 



RM=1 .DO/M 

XX=DMOD( A*1 2345 . ,M) 
DO 10,1=1,97 
XX=DMOD( XX*A+C , M) 
RANT ( I ) =XX 
CONTINUE 

XX=DMOD( XX*A+C,M) 

Y=XX 

RETURN 

END 
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TECHNICAL MEMORANDUM TM8706-4 

TO: Sid Novosad, Henry Chen, NASA/JSC 

FROM: James Dodds 

SUBJECT: SSIS Header Code 

DATE: June 19, 1987 


Warner Miller and Henry Chen recently sent me some new material on the use of a 
binary Kasami code for burst error control with the SSIS frame header. I understand that 
there is a CCSDS meeting next week (6/22), and that this is one of the topics of 
consideration. I have a few quick comments prior to a more detailed analysis. The 
assumptions made by Goddard regarding the error rate into the (85,64) binary code are: 

• Probability of burst length greater than 10 bits is zero. 

• Input error rate less than 10 - 5. 

With these assumptions, they bound the probability of header error as the probability of 2 
or more channel errors with no burst errors, or 1 channel error with a burst of length 10 or 
less. This is certainly a conservative bound, given the assumption. We might want to take 
a harder look at the burst error statistics, however. A related issue is the performance of 
the error trapping decoder. Historically (if 20 years in the past can be considered history) 
the error trapping decoder was considered because of its relative simplicity. However, 
modem high speed processors allow us to consider the use of other decoding techniques. 
The reason for bringing this up is that I don’t know what the residual error detection 
capability is with this code in the event that a burst length of 1 1 or greater occurs, or if it is 
overloaded with random errors. 


I took a quick look at an alternate, a (17,13) Reed Solomon code with 5-bit 
symbols. This code is 85 bits long, with 20 check bits, 64 address bits and 1 "freebie." It 
can correct bursts in the range of 6 to 10 bits, with six guaranteed. In addition, the residual 
error detection capability is estimated to be as follows: 


BERin 

P 

e, header 

0.5 

0.125 

0.1 

0.25 x 10" 1 

0.05 

0.50 x 10- 2 

0.01 

0.59 x 10- 4 


A (19,13) Reed-Solomon code, with a length of 95 bits, can correct bursts in the range of 
1 1 to 15 bits, with 1 1 guaranteed. The residual error detection capability is also better. 

I'll call you Monday and we can discuss pursuing this further. 
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TECHNICAL MEMORANDUM 

TO: Sid Novasad, NASA-JSC 

FROM: James Dodds 

SUBJECT: S S IS SM Codes 

DATE: September 22, 1987 


TM8709-1 


1.0 INTRODUCTION 

Axiomatix has implemented a nominal rate 1/2 (8 data rows + 8 check rows) SM-8 
decoder in order to verify the Schilling/Manela results. This memo describes details of the 
specific algorithm implemented. We are getting close to the results claimed by Manela in 
his thesis for the SM-8 code, and in fact, one variation of the decoding algorithm is 
exceeding his results. The outstanding question is, of course, the complexity of the 
hardware implementation vis-a-vis a Reed-Solomon code of comparable performance. 

The Axiomatix decoder has a checkered history; the initial decoder had performance 
several orders of magnitude worse than the Manela results. However, after working with 
the codes, and gleaning incidental pieces of information from David Manela and Don 
Schilling, we were able to approach, but not duplicate, their results. We have tried several 
sets of parity check slopes, which are not necessarily optimum. Schilling and Manela are 
not willing to provide us with the "optimum" slopes, but it appears that we are close 
enough for practical purposes. 

2.0 SM-8 CODE CHARACTERISTICS 

The specific code implemented is depicted in Figure 2.0-1. The maximum block 
size is limited by the maximum allowable storage space of the Fortran compiler object 


SSIS SM Codes 


2 


September 22, 1987 



FIGURE 2.0-1. SM-8 CODE BLOCK INCLUDING EXTRA PARITY 
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modules, 64 K, The code of Figure 2.0-1 has parity slopes of 1/0, 1/1, 1/2 1/7. A 

code with slopes consisting of the reciprocal of prime numbers was also implemented. The 
code performance is close to the 1 /n slope code, but the block size is necessarily shorter 
due to the steeper slope. The nominal code rate is 1/2, but the flush bits lower the code rate 
by a significant factor. The number of flush bits can be calculated as follows. The slope 
1/0 parity line requires no flush bits, the slope 1/1 requires 7 flush bits, the 1/2 slope 
requires 14 flush bits, and the 1/n slope requires 7n flush bits. Thus, the total number of 
flush bits is 


F = ^ 7 * i = 7 ^ n ~2 ^ n ~ = 196 bits 
i=l 

The overall code rate for the 1 /n code is R = 95 = 0 - 477 

, . 1376 

The rate of the prime number code is 2752 + 294 - U ‘ 43 

One of the ARQ strategys described in the next section lowers the effective code 
rate. The code rate loss is, of course, a function of input error rate, as the number of 
retransmissions increases over the noisy channel. 

3.0 SM-8 DECODER IMPLEMENTATION 

Two variations of an SM -8 decoder have been implemented. The basic decoder 
decodes a block of data using a single "pass," while the ARQ version decodes a block, 
checks the cumulative syndrome register against a threshold, and makes an additional pass 
if any cell is above the threshold. This variation allows us to use a lower threshold and 
pick up more errors. The number of passes must be limited, both to prevent error 
propagation, and to prevent infinite decoder loops. 

The simulation uses a pseudo-random number (PRN) generator which computes a 
uniformly distributed PRN, compares the number with the selected input error rate, and 
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stores a 1 in those bits which have a PRN less than the input error rate. A block diagram of 
the simulation is shown in Figure 3.0-1. Referencing this figure, one of eight arrays, R1 
through R8, are dedicated to each of the eight parity lines. If a parity line has odd parity, 
all bits on the diagonal line of the appropriate array are set to one. After all array have been 
set, the contents of each array are added, with the result stored in array ICNT. For 
example, a single error at row J, column I in the information block will result in 1 s being 
stored at all diagonals of R1 through R8 that intersect at I, J. Thus ICNT (J,I) will contain 
an 8 after all the R arrays are added together. 

The cell corrections are made from bottom (Row 8) to top (Row 1), and the 
threshold level for correction is initially set to 8, then 7, and so on. Using a higher 
threshold first, as per a hint by David Manela, corrects the most likely (single) errors first. 
This strategy gives a significant improvement in decoder performance. Figure 3.0-2, a 
reproduction of a plot from the Manela thesis, shown the results of the Axiomatix 
simulation relative to the Manela results. 

3.1 Modified SM-8 Decoder 

In the course of trying to reconcile the Axiomatix results with the published 
SM-8 results, a variation of the Manela decoder was implemented. The basic decoding 
algorithm is unchanged, but the threshold is lowered to a value of 4. If no other changes 
were made, this would increase the error rate due to a propagation effect. However, once 
the block is initially decoded, we interate on the block. That is, the cumulative weights are 
recomputed, and the decoding process repeats. This procedure can lead to an infinite loop 
at low SNR, so a maximum count of 2 iteration per block is imposed. If, after two 
interations some cells still contain values equal to or greater than 4, a retransmission is 
requested. This strategy reduces the code rate, but is a very effective error detection 
strategy. At an input error rate of 0.02, the unmodified decoder has an output error rate of 

3.2 x 10 -6 , while the modified decoder has an error rate of 1 x 10“^. The net code rate is 
lowered by a factor of 0.98, equivalent to a loss of 0.07 dB. 
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FIGURE 3.0-1 
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AXIOMATIX SM-8 DECODER 


(FIGURE 3.2-5 P(out) versus P(in ) R=1/2) 

FIGURE 3.0-2. AXIOMATIX SIMULATION RESULTS. 
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TECHNICAL MEMORANDUM TM8710-1 

TO: Sid Novosad, NASA/JSC 

FROM: James Dodds 

SU B J EOT: A Critique of the Schilling-Manela Codes 

DATE: October 14, 1987 

Axiomatix has been evaluating the performance of a subset of the SM codes 
presented in Dr. Manela’s thesis (Ref. 1). In addition, we have reviewed a critique of the 
thesis by Lin and Costello (Ref. 2) provided by NASA/JSC. In general, we agree with the 
Lin and Costello review of the codes: the codes are threshold decodable convolutional 
codes, and the decoding complexity engendered by the iterative decoding technique can 
cause throughput delay. One advantage claimed for the code proposed for Space Station 
use, the partially autoconcatenated (PASM) code, is the ARQ strategy permitting 
transmission of additional parity only. However, this is not unique to the SM codes. Any 
systematic code can be re-decoded with additional transmitted parity bits. This technique is 
of dubious advantage in decoding a high rate code as proposed for the Space Station. For 
example, a rate 0.9 code will, on the average, have 9 out of 10 errors in the information 
bits. Thus, it may be more efficient to retransmit the whole block rather than just parity for 
high rate codes. This is particularly true considering that each additional parity line has a 
shallower slope than the prior line, and the overall rate decreases faster than just the ratio of 
data rows to the total number of rows. 

The probabilistic nature of the decoding process was evidenced in the Axiomatix 
implementation of the decoder. Each cell of the composite matrix is compared with a 
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threshold. In the case of the SM-8 code, the threshold starts at 8, since this would be 
induced by a single error (single errors are more likely than double errors, and so on). If 
the cell equals the threshold, the appropriate input bit is toggled, and all affected cells are 
recomputed. This process repeats for each cell. Next, the threshold is lowered to 7 and the 
whole threshold comparison repeats. Note that the decoding rate will vary as a function of 
the error patterns. Two widely separated errors will cause two recomputations at the 
threshold = 8 level, whereas two errors affecting the same composite cell will require only 
1 recomputation at the threshold = 7 level. This means that the computational speed of the 
decoder must be fast enough to handle worst-case error patterns, or data will either be lost 
or output out-of-order. 

The actual computational requirements of the decoder will be governed by a 
probability distribution function which is a function of channel noise, much like the 
probability of buffer overflow of a sequential decoder. That is, the probability of 
computational failure will increase as channel noise increases, and the decoder clock rate 
must be set fast enough to provide an acceptably low probability of failure. The nature of 
the computational distribution function may have to be estimated using simulations, since 
the iterative nature of the decoding process is not deterministic. 

The Axiomatix implementation of the proposed rate 0.9 PASM code has been 
frustrated by the FORTRAN compiler 64K byte memory constraint If we were to use the 
same programming shortcuts that allowed us to rapidly implement the SM-8 code, the 
memory requirement for a single composite matrix would exceed this limit In lieu of using 
straightforward programming techniques, we would have to invoke complex bit shuffling 
algorithm which will slow down the decoder (and the programmer). An alternative is to 
program the decoder in "C" which is more suitable for this type of problem, but which will 
cause a delay due to our " learning curve" with the C language, or tie into the TCD VAX 
with a modem and remote terminal. 
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In conclusion, barring any new revelations, the SM codes do not appear to offer 
any significant advantage over other threshold decodable convolutional codes, and the 
dynamic behavior of the decoder may preclude its use over high-speed links. The 
proposed (255,223) Reed-Solomon code can provide the required error protection. This 
code can correct a minimum of 16 bit errors per block, with a residual error detection 
capability guaranteeing a probability of undetected error less then 10 . Since the chip set 

exists to implement this decoder, we feel that it is logical to use the Reed-Solomon code 
rather than the less mature SM code. On the positive side, the algorithm is simple and easy 
to understand; this class of codes presents an excellent introduction to multiple error 
correcting codes without resorting to the use of higher order algebraic fields and could well 
serve as a first example in a coding course. 
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FROM: James Dodds 

SUBJECT: WANDERER TRACKER PERFORMANCE GOALS 


DATE: February 17, 1987 


1.0 INTRODUCTION 

This memo addresses two of the significant concerns regarding candidate systems 
to detect wandering behavior and track a wanderer-cost and effectiveness. We feel that 
minimum cost can be attained by using devices based wholly or predominantly on 
commercial products and technology wherever possible. The use of a low cost intrusion 
alarm as part of the perimeter detection system is discussed in a subsequent paragraph. 
Effective wanderer tracking forces us to consider conflicting requirements, for example, 
long range implies high power transmitters, which in turn may require excessive battery 
power and short battery life. Again, some suggested tradeoffs and possible solutions are 
discussed in subsequent paragraphs. 

2.0 TECHNICAL CONSIDERATION FOR EFFECTIVE WANDERER 

TRACKING 

The following discussion assumes that a potential wanderer will be equiped with a 
radio-frequency transmitter, a perimeter alarm will detect the presence (or absence) of the 


February 17, 1987 

Wanderer Tracker Performance Goals 
Page 2 


RF signal, and tracking capability will be provided by RF direction finding. Topics 
considered include: 

• Minimization of battery power requirements 

• Directionality of signal 

• Signal penetration of structures 

• Range 

• Available technology 

The first topic, minimization of power, has lead us to consider a system with a 
wanderer's transmitter in the "normally off" condition. With this system, depicted in 
Figure 2.0-1, the transmitter is held in the off state by the presence of an acoustic or 
electronic field provided, for example, by a commercial intrusion alarm. Such alarm 
systems can be purchased for under $200 per unit (Ref. 1) and provide an inexpensive 
source of microwave or acoustic energy. A receiver on the wanderers unit detects the 
presence of the energy to suppress the wearers transmitter. The wearers unit can also be 
equipped with a "dead-mans" switch which will trigger the transmitter if the unit is 
removed by the wanderer. Battery life should be extended substantially, since power is 
only needed for the receiver during periods of normal activity. This system has the added 
advantage that only a single, inexpensive transmitter is needed for each area to be 
monitored, regardless of the number of doors and windows, thus eliminating antennas, 
loops, or receivers at each portal. The same system can be used to monitor outdoor activity 
by careful placement of transmitters. 

The wearer units can be tested by periodically shutting down the area transmitters. 
The system can be considered fail-safe, since the failure of any component of the system 
can be detected automatically. Failure of an area transmitter will cause all wearer units to 
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respond, failure of one wearer unit (assuming each has a unique ID) will be evident during 
periodic checkout, as will failure of the base receiver when all units fail to respond during 
checkout. However, to distinguish individual wearer units during checkout, we must use a 
modulation which cannot be jammed by alternate wearers. A baseline modulation for 
consideration is a pulse-position modulator with time division multiplex pulse timing 
unique to each wearer. The wearer-unique code provides identification as well as low 
mutual interference. 

In addition to the three components disscussed, the energy source, the wearer unit 
• and the base receiver, a fourth component is needed for direction finding. For simplicity, 
we assume the direction finding unit is not required to identify a specific wearer unit once 
identification is made at the base receiver, but rather it is sufficient to track the energy of the 
wearer unit. Thus, the portable unit can be consideraly simplified. At this point we need to 
examine another tradeoff area, that of the frequency band utilized by the wearers unit. As 
depicted in Figure 2.0-2; two types of direction finding can be employed; (1) locate peak 
of pattern and (2) locate null of pattern. The requirement for a portable, hand-held receiver 
dictates that a hom antenna at EHF be used for the peak-of-pattem method, while a loop or 
ferrite antenna can be used at HF to locate with a null-of-pattem. Direction finding on a 
null is much less satisfactory for an inexpenrienced operator, particularly with the broad 
nulls at HF using a loop or ferrite. In addition, it is difficult to distinguish the 180° 

directional ambiguity, as shown in Figure 2.0-2; at EHF, a highly directional beam can be 
formed with a moderately sized antenna. Thus, it would seem that frequencies at K a band 

and above would be prime candidates for the wanderer system if the low cost technology is 
available. However, these frequencies are easily blocked by structures and foliage, 
whereas HF energy is not attenuated as readily by obstacles. A hybrid approach would be 
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to use a dual frequency system with coarse direction finding at HF in the presence of 
obstacles, and accurate direction finding at EHF. One frequency band could be 
implemented as a simple CW source, with modulation on the other band to provide the 
wearer ID. The complexity of the wearer unit would be slightly increased, with the worst 
impact on the portable direction-finding receiver which would probably need to 
accommodate both peak and null detection capability. 

2.1 Wearer Identification 

The problem of discriminating between different wanderers is exacerbated when, 
during self test, all wearer untis are turned on as the area transmitter is turned off. Some 
form of orthogonal modulation must be used to prevent mutual interference, overloading of 
the area receiver's front end, and suppression of weak signals from more distant wearers. 

An example of a modulation scheme with a relatively simple demodulation, not 
requiring sophisticated spread spectrum demodulation or correlation receivers, would have 
the wearer units time synchronized to the area receiver clock. This system is depicted in 
Figure 2.1-1. Pulses from the area transmitter periodically reset counters in the wearer 
units and in the area receiver. Each wearer unit clock then counts to a preset value, Tj, 

different for each unit, before commencing transmission of alarm pulses. The minimum 
counter period of each wearer unit is greater than the reset pulse period T r from the area 

transmitter, thus preventing transmission by wearer units within range of the area 
transmitter. Once activated, the wearer units transmit at a pulse interval T r The stability of 

the free-running oscillator in the wearer unit need only be good enough to provide adequate 
time to identify the unit by preventing the pulse epochs from drifting into the time slot 
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corresponding to another wearer. For example, a system capable of accommodating 50 
wearer units might have the following parameters: 


Pulse duration : 

10 jisec 

Guard time : 

10 jisec 

Repetition : 

1 KHz 

Number of time slots : 

50 


If we assume that the receiver can tolerate a 5 |i.sec drift (1/2 a guard time), and lsecond 
(1 thousand pulses) is an adequate time to identify a wanderer, then the transmitter stability 
requirement is about 5 x 10 ^, which is not difficut to achieve. 

Installation having more than one area transmitter will require that the transmitters 
be synchronized, since it is assumed that only one area receiver will be needed. The 
transmitters can be hardwired together, or the synchronization signal can be superimposed 
on the AC powerline. Note that simplification can be effected for small installations not 
requiring wearer identification, but the economy of mass production would probably dictate 
a common wearer unit, whether or not identification is desired. 

3.0 CONCLUSIONS 

We have proposed a candidate wanderer tracking system with features which 
provide fail-safe operation, minimizes battery power requirements, and can provide good 
tracking capability. A more detailed investigation will provide insight into the cost 
effectiveness of considering a dual-band wearer unit or whether a compromise single 
frequency system would provide a better cost/performance ratio. 
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1. INTRODUCTION 


The BOSS Lightwave Module Library contains modules which are useful for 
analyzing certain single mode fiber digital and analog communication links. 
The library features a general Single Mode Fiber module, an Avalanche Photo- 
detector module, a PIN Photodetector module, a Semi— Analytic Error Rate Esti- 
mator module for digital On-Off Keyed (00K) systems, and a Power Series User 
module, plus various other lower-level and internal modules (listed below). 

In addition, two example systems are included to demonstrate how to use 
the modules. 

The following is a list of the modules in the Lightwave Module Library 
arranged according to group name: 


ANALOG MODULATORS 

LASER (POWER SERIES) 

BASIC BULDING BLOCKS *TYPE/UNITS CONVERSION* 
REAL TO DBM 

CALIBRATION DEVICES /METERS 
DELAY METER (REAL) 

OPTICAL AVERAGE POWER 
PRINT AVERAGE OPTICAL POWER 

CHANNELS 

SINGLE MODE FIBER (LINEAR) 

DIGITAL SOURCES 

OOK OPTICAL SOURCE 

ESTIMATORS 

NOISE BW .IMPULSE INJECT (REAL) 

OOK_ERROR RATE ESTIMATOR 

ESTIMATORS INTERNALS* 

COUNT ERRORS_& GENERATE STOP 

DELAY TX_SIG 

NOISE_BW COMPUTE (REAL) 

OOK ERROR PROB CALCULATE 
OOK PRINT 

OOK_POINT TO_DISTANCE 
FILTERS 

BUTWTH FILTER (REAL) 

FILTERS * INTERNALS* 
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2ND ORDER HR SECTION (REAL) 
TAPPED DELAY LINE CELL (REAL) 

NOISE AND INTERFERENCE 

AVALANCHE PHOTO-DETECTOR 
AVALANCHE PHOTODIODE RAN_GEN 
PIN PHOTO-DETECTOR 


The top-level modules are discussed in more detail below. Many of the modules 
in the lightwave database are lower-level internal modules and are not expli- 
citly discussed in detail here. Of course, on-line documentation is available 
for these modules as well as all other modules In BOSS. On-line documentation 
for the modules discussed below is provided in the Appendix. 

It should be pointed out that in general, modeling of lightwave systems 
is rather problem dependent. It is difficult to separate the functionality of 
the various components of lightwave systems into modules which are totally 
independent of other modules in the system. This is the reason why modules 
for other lightwave system degradations are not currently available in the 
Lightwave Module Library. 


2 . LASER MODULE 

This module models the input current to optical power output characteris- 
tic for a laser diode above threshold using a power series nonlinearity. The 
module is intended for modeling analog lightwave communication systems where 
the nonlinearity of the laser is the predominate degredation (in the laser). 
The maximum order of the polynomial is five, and is specified either by har 
monic levels (in decibels relative to the fundamental) or by direct entry of 
the polynomial coefficients themselves. 

The parameters of the module are: 


HARM/COEF 

TEST TONE MAGNITUDE 
5TH HARMONIC LEVEL 
4TH HARNOMIC LEVEL 
3RD HARMONIC LEVEL 
2ND HARMONIC LEVEL 
EFFICIENCY 
BIAS POWER 
BIAS CURRENT 


2 - 


3. OOK OPTICAL SOURCE 


This module outputs an On-Off Keyed (OOK) non-return-to-zero signal used 
to drive a lightwave fiber module* The random binary data on-amplitude is 
controlled by the value of the parameter "AVG OPTICAL POWER", which has units 
of decibel milliwatts (dbm). Initialization code reads in this value and 
computes the on-amplitude assuming that the of f -amplitude is zero* 


4. SINGLE MODE FIBER MODULE 

The Single Mode Fiber modules provided in the BOSS lightwave Module 
Library is an implementation of a model given by Duff [1]. The baseband 
transfer function of a single mode fiber used in this module is given by the 
following equation: 

00 

H(f ) - J S( A)L( X) e" jaiLT(X) dX ( 1 ) 

— 00 

where 

a) ** 2irf 

L ■ fiber length 

S( A) = optical source spectrum vs wavelength 

L(A) - l/(fiber loss) vs wavelength 

T(X) * fiber group delay vs wavelength 

(antiderivative of dispersion vs wavelength) 

In the BOSS module, an initialization subroutine calculates the transfer 
function given in Equation 1, takes the inverse FFT, and truncates the result- 
ing impulse response. The impulse response is then convolved with the incom- 
ing signal using a tapped delay-line, implemented in block diagrams. 

The integral in Equation 1 is implemented using Simpson's rule* The 
interval specified by the lower and upper limits of integration should be 
large enough to include insignificant portions of the source spectral 
envelope, S(A). 
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The module has three "modes of operation: 


1. Both source spectral envelope and fiber chromatic dispersion speci 
fied by closed form approximate functions; 

2. Source spectral envelope specified by tabular data file, and closed 
form approximation for chromatic dispersion; 

3. Both source spectral envelope and chromatic dispersion specified by 
tabular data files- 


In all three modes, the fiber loss per kilometer is specified by tabular data 
files. 


The closed form approximation for the source spectrum as a function of 


wavelength is given by 
1 


S(X) 


a 2 / 2 tt 


(x - x s ) 

exp( t ) 

2 a 


( 2 ) 


where 


o g •» rms spectral width 

X g “ source center wavelength 


The chromatic dispersion closed form approximation as a function of wave 


length is (assuming a silica fiber) 


dT( X) 
dX 


X - X r 


(s/c) (- 


(3) 


where 


s -a unitless constant 

c * speed of light 

Xq « zero-dispersion wavelength 

The group delay is then 

x x o" \ 

T( X) - (s/c)[£n(4-) + j—] 

Xq ^ 


(M 
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. n1v to the closed form approximation, while 

Some of the parameters only apply Although BOSS 

others only apply to functions specified By ^ ^ che ^ameter 

always requires a valid entry for each param , module mode of 

mo, he Ignored by the module If It does not apply 

operation- 

. „ hon the closed form approximation 
The following parameters are read on y 

for chromatic dispersion is used: 

dispersion adjustment constant 

ZERO DISPERSION WAVELENGTH 

Parameters read only when the closed for. approximation for the source 
spectrum is used are: 

RMS SOURCE SPECTRAL WIDTH 
SOURCE CENTER WAVELENGTH 

aoolv only when chromatic dispersion as a func- 
The following parameters PP y 
tion of wavelength Is specified In a tabular data «• 

GROUP DELAY POINT (WAVELENGTH) 

GROUP DELAY POINT (DELAY) 

DISPERSION FILENAME 

, A onlv when the source spectral envelope is 

The following parameter is read only 

specified by a tabular data file: 

SOURCE SPECTRAL ENVELOPE FILENAME 

The remainder of the parameters are used regardless of the mode 
tion of the fiber module, and are .8 follows: 


FIBER LENGTH 

# OF POINTS IN FIBER H(F) 

# OF PARTITIONS FOR INTEGRATION 
INTEGRATION UPPER LIMIT 
INTEGRATION LOWER LIMIT 

FIBER H(T) SAVE FILENAME 
FIBER LOSS FILENAME 
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recalculate fiber H(T) 

MODE OF OPERATION 


Computing the fiber Impulse response for a given set of parameters can 
require from five to forty-five minutes of cpu time. In many applications, 
the same fiber characteristics may be used from one simulation no the next. 
To prevent wasting time recalculating the fiber Impulse response in these 
cases, the module writes the calculated impulse response to the "FIBER R(T) 
SAVE FILENAME" when the value of the parameter "RECALCULATE FIBER H(T) is set 
to TRUE. In subsequent simulation runs, the value would be set to FALSE, 
which causes the program to read the impulse response from the file specified 
by the parameter "FIBER H(T) SAVE FILENAME." 


To Insure that the parameters passed to the fiber module In a given 
simulation tun are the same as the ones used to calculate the fiber Impulse 
response in the save file, the module writes out a "verification record 
(explained In the BOSS user's .annual) at the bottom of the file. This ver 
cation record contains the numeric values for all parameters which were re- 
quired to generate the fiber Impulse response. Before the fiber module win 
successfully read the verification record, all of the parameters of the vert 
f ication record must agree with the parameters passed Into the fiber module. 
This prevents reading fro. a fiber impulse response file which was generated 
with parameters different fro. those values shown on the simulation parameter 

list. 


5. AVALANCHE PHOTODETECTOR 

This module Is an Implementation of an Avalanche PhotoDetector (APD) based 
on the Poisson and McIntyre distributions [2]. 

The number of primary charge carries (hole-electron pairs) in an AP 
modeled by a Poisson process, with the average number per At given by 

* - h p mc (t) • “ + ”d ’ “ <5> 


where 


i 
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P inc(0 

n 

h 

u 

n d 


- incident optical power 
■ quantum efficiency 
= Plank’s constant 
» optical frequency 

= number of dark charge carriers/second 


The exact number of primary carriers, np, is drawn from this Poisson* Each 
primary charge carrier pair generates secondary charge carriers (providing the 
gain) in a random way. The discrete probability density of the number of 
secondary charge carriers plus the original primary charge carrier was found 
by McIntyre to be 


PCrij* m) 


(I -*)"- 1 rOjqj) 


l+k(m-l) 

1 1-k J 


[ l+k(m-l ) ] (m-1)! T (- 
Hence for a given At, the output current of the APD is 


G_i ( m_1 > 

1V1 


i(t) - 


C np 

— 1 
At i=i i 


( 6 ) 


( 7 ) 


where 


C e ■ the charge of an electron 

n i - number of secondary charge carries plus the primary generated in At 
seconds 

The actual random number generation is done by the primitive module 
AVALANCHE PHOTODIODE RANJ3EN . The McIntyre distribution is implemented using 
the table lookup method. The following parameters to the APD module are 
related to the distribution: 


AVG APD GAIN 
APD "K" PARAMETER 
MAX CDF PROBABILITY 
CDF TABLE SIZE 
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The cumulative distribution function (cdf) is generated once at the beginning 
of the simulation from values given in the above parameters. A binary search 
is used to produce variates from the table during the simulation. "MAX CDF 
PROBABILITY" and "CDF TABLE SIZE" control the accuracy as explained in the on- 
line documentation in the Appendix* 

Note that usage of this module in shot noise limited cases where the 
error rate however, is low (e.g., 10 -9 ) require to much simulation time to be 

practical. It does accurately model the APD, and could be useful for waveform 

analysis. 

The Poisson random number generator uses the acceptance-rejection method 
to output a deviate for the given mean [3]. 

The remaining module parameters are: 

APD INITIAL SEED 
DARK CURRENT (NANOAMPS) 

QUANTUM EFFICIENCY 
SOURCE CENTER WAVELENGTH 

These parameters are further discussed in the on-line documentation given in 
the Appendix. 

6. PIN PHOTODETECTOR 

A PIN photodetector module is also provided, and is the same as the APD 
module, except that there is no secondary charge carrier generation. Hence, 
for the PIN, the number of primary charge carriers is the same as in Equation 
(5). The output current Is given by 

C • np 

w “-V- 

where 

C e « the charge of an electron 

np * the number of primary charge carriers generated in At seconds 
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The parameters are the same as the APD, except there are no parameters per- 
taining to the McIntyre distribution. 

7. 00K SEMI-ANALYTIC ERROR ESTIMATOR 

This module uses the combined analytic-simulation technique for estimat- 
ing the bit error rate of digital 00K systems [4]. This module is used in 
conjunction with the NOISE BW IMPULSE INJECT (REAL) module which is inserted 
in the block diagram at the point where the additive white Gaussian noise 
enters the system. The Error Estimator can be used only if the modules down- 
stream of the point where additive white Gaussian noise enters the system are 
linear. 

When set up properly, this module automatically determines the equivalent 
noise bandwidth of the cascade of modules downstream from the additive white 
Gaussian noise source, and estimates the amount of delay between the transmit- 
ted and received bitstreams. Results from the calibration phase are available 
for viewing in the Postprocessor using the View Printed Output option. 

After the module has finished the calibration phase, it begins the error 
estimation phase. For each received sample, the distance to the decision 
threshold is measured (the decision threshold is a parameter to the module). 
The estimated noise bandwidth and the user specified value of noise power 
spectral density are used to determine the variance of the Gaussian noise at 
the input to the estimator. The measured distance to the threshold is divided 
by the square root of this variance and a Gaussian tail probability (Q 
function) lookup is performed. The resulting error probability for this 
received sample is averaged with the error probability for the corresponding 
sample in all the bits in the estimation phase of the simulation. 

The results are plotted in a multi-dimensional plot for viewing in the 
Postprocessor. If the BOSS Parameter Iteration option is NOT used, the only 
independent axis for the plot is NORMALIZED BIT TIME (MICROSECONDS). This 
plot is useful for determining the sensitivity of the error rate to timing 
errors in the sample and decide portion of the receiver- 

If Parameter Iteration is used, the extra independent axis can be used to 
plot the variation of a parameter (where each iterated value is used in a 
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separate simulation) on the same multi-dimensional error rate plot (see the 
BOSS user's manual for a more detailed description of Parameter Iteration). 


This Estimator module has an option which allows it to count actual 
errors which occur due to shot noise or severe ISI and store these in a multi- 
dimensional plot with NORMALIZED BIT TIME (MICROSECONDS) as the independent 
axis. The user can input the minimum number of errors which must be counted 
for every sample in the bit interval. When this number of errors has been 
exceeded for all samples in the bit interval, the Estimator module outputs a 
signal which can be used to drive a TERMINATE SIMULATION module. This module 
shuts down the simulation by setting time equal to STOP TIME. This feature 
allows the Estimator to be used in situations where there is a significant 
amount of shot noise (enough to cause the sample to cross the decision thres- 
hold) in addition to thermal noise. 


- From 10 to 100 errors should be counted to obtain reasonable confidence 
in the error probability estimate, so for extremely low error rates direct 
simulation is impractical. In these cases, an alternate simulation technique 
such as Importance Sampling should be used [4]. 


The parameters for the Error Estimator module are 


CALIBRATION START TIME 
NOISE BW CALC DURATION 

# OF SAMPLES IN CORR SUB SEQ 
MAX DELAY TO CALCULATE 

# OF SYMBOLS FOR ERROR EST 

# OF SAMPLES /SYMBOL 
THERMAL NOISE PSD 
DECISION THRESHOLD 
MINIMUM # OF ERRORS 
PROB OF ERROR AXIS LABEL 
INDEPENDENT AXIS #2 LABEL 
INDEPENDENT AXIS it 2 VALUE 


8. EXAMPLES 

Two examples are presented in this section which use most of the mod 
in the Lightwave Module Library. The examples are mainly for illustration, 
and are not intended to show the only valid usage of the modules. 
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The first example is an end to end single mode fiber system, with on-off 
keying (block diagram shown in Figure 1). The simulation parameters are shown 
in Figure 2. In this example, the fiber module uses the closed form 
approximations for both source spectral envelope and dispersion given in 
Equation 2 and Equation 3 (hence the parameter value CLOCLO). Since the 

parameter “RECALCULATE FIBER H(T)" is FALSE, the module will read the fiber 
impulse response from the file HIGHRES.SAV. The decision threshold value of 
0.125 was determined from a previous simulation run. 

The system level parameter values for the error estimator, PIN 
photodetector, fiber, and source modules are shown respectively In Figure 3. 
The INDEPENDENT AXIS #2 parameters are not used in this example, but are 
available for Parameter Iteration if desired. Notice that the three tabular 
filenames in the fiber module are all exported to the same filename. Since 
the fiber module is using closed form approximations for source spectrum and 
dispersion, it ignores the file name given for the respective tabular file 
parameters. 

Data written out into the Postprocessor by the error rate estimator and 
PRINT AVERAGE OPTICAL POWER modules is shown in Figure 4 (this data can be 
read by selecting the View Printed Output option). The estimator writes out 
the beginning and ending times for calculating the noise bandwidth of the 
receiver and the start time for measuring the delay between the transmitted 
and received bitstreams. Also, the time which the actual error rate estima- 
tion begins is printed along with the minimum value of STOP TIME required to 
process the requested number of bits. The calculated delay In samples is 
displayed along with the single sided noise bandwidth in Megahertz. This 
information Is useful when initially setting up a run with the error rate 
estimator module. The PRINT AVERAGE OPTICAL TOWER module writes out the 
optical power In both dBm and Watts, as seen in Figure 4. 

A plot of Log (Probability of Bit Error) versus normalized symbol time 
for this example is shown in Figure 5. From the plot, sensitivity to timing 
errors in the receiver can be easily seen. If the independent axis # 2 
parameters had been used with Parameter Iteration, this plot would have had an 
additional dimension (the additional dimension would be orthogonal to the 
plane of the page). 
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The block diagram for the second example is shown In Figure 6. In this 
example, the three QPSK channels are Frequency Division Multiplexed and the 
composite signal used to intensity modulate a laser. The QPSK portion of the 
system is simulated using complex envelope representation. A Pseudocarrier is 
introduced to interface this complex envelope signal to the lightwave portion 
of the link, which operates on real signals. The pseudocarrier is the removed 

at the output of the lightwave link. From this point on, the signal is again 
in complex envelope form. 

Figure 7 shows the simulation parameters and the system level parameters 
for the Laser module. In the example, a normalized baud rate of one is used, 

with a pseudocarrier frequency of 10. The laser operates with a third 
harmonic level of 14 dB below the fundamental. Thus when the laser has an 

input tone with a magnitude of 10, (the test tone magnitude) the output will 

have a third harmonic value 14 dB below the output at the input tone 

frequency* 

In this particular simulation, the nonlinearity of the laser is the 
predominate degradation, therefore the fiber module is not used (since it is 

linear). A gain module is used for scaling in place of the fiber module. 

In Figure 8, the magnitude spectrum of the laser output in shown. Figure 
9 displays a scatter plot and a time plot of the receive filter output. 

It is hoped that the examples provide some assistance in using the 
Lightwave Module Library. The Block diagrams and parameter listings should 

provide a starting point for using the modules to analyze other lightwave 
system configurations. 
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Figure 3. Parameter values at the system level -or: 
Error Rate Estimator, PIN Photodetector, 
Single Mode Fiber, and OOK Source. 
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Your choice of parameters for the 
30K._ ERROR ESTIWTCR result in the 
following comoutsd times : 
noise bw start time, noise ow stop time 
: 7. 1400002E-02 0. 1429000 

delay est start time, error *st start time 
s 0.2142000 0.2213432 

error eet minimum stoo time • 3.071352 

On-Off keyed error eet tx to rx delay 25 (sameles; 

On-Off keyed error eat noise bv 104.2952 (single sided, her zzi 

There were 394 bits used in ook error estimation. 

Incident average ootical power -38.90030 Decibel milliwatts 
Incident average optical power 1 . 2646481E-07 watts 


Figure 4. Postprocessor printed output for digital lightwave 
example . 



Figure 5. Bit error rate plot for digital example. 
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Figure 7. Simulation parameters for FDMA example and the 
system level parameters for the Laser module. 



Figure 8. Laser module output spectrum for the FDMA example. 
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11. APPENDIX 


This Appendix contains the on-line documentation for the following top- 
level lightwave modules: 

AVALANCHE PHOTO-DETECTOR 
LASER (POWER SERIES) 

00K OPTICAL SOURCE 
OOK_ERROR RATE ESTIMATOR 
PIN PHOTO-DETECTOR 
SINGLE MODE FIBER (LINEAR) 
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AVALANCHE PHOTO-DETECTOR 


Page 1 
8 June 1987 


MODULE NAME: AVALANCHE PHOTO-DETECTOR 

GROUP: NOISE AND INTERFERENCE 

DATABASE: DUAO: [ BOSS. ADD_DB. LIGHT] 

AUTHOR: TOWNSEND 

CREATION DATE: 18-May-1987 16:41:03 


DESCRIPTION: 

This module implements an avalanche diode photo detector. 
The input signal should be an optical power signal (Uatts). 
Any small amount of negative power on this input will be 
ignored by. the module (this may occur in an FFT filter, 
for example). Shot noise effects are simulated using 
an AVALANCE PHOTODIODE RAN_GEN module.. The output current 
is in units of microamps. 


REVISIONS: 


Author : SYSADMIN 

Date : 6-Jun-1987 1:50:40 

Description: 

Module RECOMPILED by CONVERTER. 

Author : SYSADMIN 

Date : 18-May-1987 16:41:03 

Description: 

Module RECOMPILED by CONVERTER. 


INPUT SIGNALS: 

OPTICAL POWER IN Type: REAL 

Lower Limit: -1.7e38 

Upper Limit: 1.7e38 

The incident optical power in watts. 


OUTPUT SIGNALS: 

OUTPUT CURRENT Type: REAL 

Lover Limit: 0-0 

Upper Limit: 1.7e38 

The avalanche photo diode output current in MICROAMPS. 


PARAMETERS: 

APD INITIAL SEED 
Lover Limit: 1 


Type : INTEGER 
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AVALANCHE PHOTO-DETECTOR 


Upper Limit: 2147483647 

Initial seed for random number generator. Should be 
a large odd integer. 

AVG~APD GAIN Type: REAL 

Lower Limit: 3.0e-39 

Upper Limit: 1.7e38 

The gain of the apd (average number of total charge 
carriers for each primary charge carrier). 

APD "K" PARAMETER Type: REAL 

Lower Limit: 0 

Upper Limit: 0.99999 

The apd "k" parameter (ionization ratio). 

DARK CURRENT (NANOAMPS) Type: REAL 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

The PIN dark current in nanoamps. 

QUANTUM EFFICIENCY Type: REAL 

Lover Limit: 0.0 

Upper Limit: 1.0 

The quantum efficiency of the pin diode. 

SOURCE CENTER WAVELENGTH Type: REAL 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

Center wavelength of the source in micrometers. This is 
used to determine the optical frequency. 

MAX CDF PROBABILITY Type: REAL 

Lover Limit: 0.8 

Upper Limit: 0.999999 

This parameter specifies the maximum probability used 
when generating the discrete cumulative distribution 
function. Ie, this is the maximum probability accumulated 
during the generation of the cdf (the probability is set to 
1 at this point in the cdf). A typical number^ 
for this parameter is 0.9999. If this number is set 
too close to l f the routine may never finish generating 
the cdf due to roundoff error, hence it is usually 
best not to exceed 0.9999. 


Page 2 
June 1987 
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AVALANCHE PHOTO-DETECTOR 


Page 3 
8 June 1987 


CDF TABLE SIZE Type: INTEGER 

Lover Limit: 1 

Upper Limit: 5000 

This parameter defines the number of elements in the 
cumulative distribution function (cdf) for the McIntyre 
(internal) random number gen. This number should 
be as large as the largest number in the internally 
generated cdf lookup table, up to a maximum of around 
1000. Note: 1000 is always a reasonable value. 

So, if for example, the average apd gain is 10, an 
the "k" parameter is equal to .2, the largest 
number used in the cdf table is around 20 to get "max 
probability" equal to 0.9999. So a cdf table size of 20 
would be large enough in this example. The largest number 
used in generating the cdf is read using the POST PROCESSOR 
VIEW PRINTED OUTPUT option. If the cdf table size is 
larger than necessary, the random number generator will 
be slower (because of the larger table to search). 

For situations where the cdf table requires values >1000 
or so, the size should be fixed at 1000, since the table 
will do a linear interpolation if necessary. 


COMPUTED PARAMETERS: 

ELECTRON CONST Type: REAL 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

Electron constant, computed. 


DARK CONST Type: REAL 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

Dark constant, computed. 


POWER CONST Type: REAL 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

Power constant, computed. 


MODULES USED IN BLOCK DIAGRAM: 
SELECT 
>= REAL 
ADDER 

MULTIPLIER 
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AVALANCHE PHOTO-DETECTOR 


Page 4 
8 June 1987 


CONST GEN 

AVALANCHE PHOTODIODE RAN_GEN 

INITIALIZATION CODE: 

Subroutine: avalanche_photo_detector 

Arguments: 

QUANTUM EFFICIENCY 
SOURCE CENTER WAVELENGTH 
DARK CURRENT (NANOAMPS) 

POWER CONST 
DARK CONST 
ELECTRON CONST 

Description: 

This routine calculates 3 constants: "power const", 

"dark const", and "electron const". "Power const" equals 
the product of "source center wavelength", "quantum 
efficiency", and "dt", divided by 1.989e-13. "Power 
const" converts the input optical power units to average 
number of charge carriers. "Dark const" is the number 
of dark charge carriers in delta t seconds. "Electron 
const" convert:' charge into current in microamps. 
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LASER (POWER SERIES) 


MODULE NAME: LASER (POWER SERIES) 

GROUP: ANALOG MODULATORS 

DATABASE : DUAO : ( BOSS . ADD_DB . LIGHT 1 

AUTHOR: TOWNSEND 

CREATION DATE: 26-Mar-1987 9:56:13 

DESCRIPTION: 

This module models the optical power output vs current 
input characteristic for a laser diode above threshold. 

The laser is represented by an AM-AM, memoryless non- 
linearity which is approximated by a fifth order poly- 
nomial. Laser diodes are biased for analog applications 
above the threshold (lasing region). The input current, 
output optical power, and slope of the characteristic (the 
output efficiency) at the bias point must be given. The 
coefficients of the remaining terms in the polynomial 
approximation may be entered directly (if "harm/coef" is 
.false.) or as the levels of the 2nd through 5th harmonics 
(if ""harm/coef" is .true.), measured relative to the 
fundamental with a test tone of peak magnitude "test tone 
magnitude". 

The polynomial approximation assumes that the input 
current remains above the threshold current. If the input 
drops below this, the calculated output power may be <0. 

In this case, the output power is set to 0, since power 
output cannot be negative. In this case, the results will 
not be accurate for samples near or below threshold. 

REVISIONS: 


Author : SYSADMIN 

Date : 6-Jun-1987 2:04:28 

Description: 

Module RECOMPILED by CONVERTER. 

Author : SYSADMIN 

Date : 26-Mar-1987 9:56:13 

Description: 

Module RECOMPILED by CONVERTER. 


INPUT SIGNALS: 

INPUT CURRENT Type: REAL 

Lover Limit: 0.0 

Upper Limit: 1.7e38 

Laser input current. Should always be greater than the 
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LASER (POVER SERIES) 


Page 2 
8 June 1987 


threshold current. 


OUTPUT SIGNALS: 

OPTICAL PVR OUT Type: REAL 

Lover Limit: 0.0 

Upper Limit: 1.7e38 

Laser optical power output. 


PARAMETERS: 

HARM/COEF Type: LOGICAL 

Lover Limit: NIL 

Upper Limit: NIL 

When this logical parameter is .true., the values of 
parameters "XXX harmonic level' 1 are interpreted as harmonic 
levels. When .false., they are interpreted as coefficents 
of the 5th order polynomial (efficiency is the first order 
term). 


TEST TONE MAGNITUDE Type: REAL 

Lover Limit: 3.0e-39 

Upper Limit: 1.7e38 

Peak magnitude of sinusoidal input current for which the 
harmonic levels were determined. 

This parameter value is not used if M harm/coef H is .false. 

5TH HARMONIC LEVEL Type: REAL 

Lover Limit: -1.7e38 

Upper Limit: 1.7e38 

Level of fifth harmonic (in dB) relative to the 
fundamental. If the parameter "harm/coef ,f is .false., 
this parameter is the coefficient of the fifth order 
term in the polynomial. 


4TH HARMONIC LEVEL Type: REAL 

Lover Limit: -1.7e38 

Upper Limit: 1.7e38 

Level of fourth harmonic (in dB) relative to the 
fundamental. If the parameter "harm/coef " is .false., 
this parameter is the coefficient of the fourth order term 
in the polynomial. 


3RD HARMONIC LEVEL Type: REAL 

Lover Limit: -1.7e38 
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LASER (POWER SERIES) 
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8 June 1987 


Upper Limit: 1.7e38 

Level of third harmonic (in dB) relative to the 
fundamental. If the parameter "harm/coef " is .false., 
this parameter is the coefficient of the cubed term in 
the polynomial. 

2ND HARMONIC LEVEL Type: REAL 

Lover Limit: -1.7e38 

Upper Limit: 1.7e38 

Level of second harmonic (in dB) relative to the 
fundamental. 

If the parameter "harm/coef" is .false., this parameter 
is the coefficient of the squared term in the polynomial. 


EFFICIENCY Type: REAL 

Lover Limit: 3.0e-39 

Upper Limit: 1.7e38 

Slope of the output pover vs input current characteristic 
at the bias point. The units of this parameter must 
correspond to the units of the output pover/input current. 

BIAS POVER Type: REAL 

Lover Limit: 3.0e-39 

Upper Limit: 1.7e38 

Pover output at bias point (corresponding current specified 
by the parameter "bias current"). Units should be same as 
the output power signal. 

BIAS CURRENT Type: REAL 

Lover Limit: 3.0e-39 

Upper Limit: 1.7e38 

The units of this parameter should be the same as the 
input current, and the parameter "test tone magnitude". 


COMPUTED PARAMETERS: 



A2 

Lower Limit: 
Upper Limit: 

-1 . 7e38 
1.7e38 

Type: 

REAL 

Computed. 

A3 

Lower Limit: 

-1 . 7e38 

Type: 

REAL 

Upper Limit: 

1.7e38 
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LASER (POWER SERIES) 
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Computed . 


A4 

Lower Limit: -1.7e38 

Upper Limit: 1.7e38 

Computed. 


A5 

Lower Limit: -1.7e38 

Upper Limit: 1.7e38 

Computed. 


Type: REAL 


Type: REAL 


MODULES USED IN BLOCK DIAGRAM: 

SUBTRACTOR 

GAIN 

MULTIPLIER 
ADDER 
>= REAL 
CONST GEN 
SELECT 

INITIALIZATION CODE: 

Subroutine: laser_pwr_series 

Arguments: 

EFFICIENCY 

2ND HARMONIC LEVEL 

3RD HARMONIC LEVEL 

4TH HARMONIC LEVEL 

5TH HARMONIC LEVEL 

TEST TONE MAGNITUDE 

HARM/COEF 

BIAS POWER 

A2 

A3 

A4 

A5 

Description: 

This subroutine calculates the coefficients of the second 
through fifth order terms of the polynomial approximation 
of the power output vs current input characteristic of the 
laser diode. These values are determined from the levels 
of the 2nd through 5th harmonics (relative to the 
fundamental), the magnitude of the test tone, the 
output efficiency (ie, slope of the characteristic at the 
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bias point). If "harm/coef" is .false., then the entered 
harmonic level parameters are used as the polynomial coef- 
ficients . 
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OOK OPTICAL SOURCE 


Page 1 
8 June 1987 


MODULE NAME: OOK OPTICAL SOURCE 

GROUP: DIGITAL SOURCES 

DATABASE: DUAO: [ BOSS. ADD_DB. LIGHT] 

AUTHOR: TOWNSEND 

CREATION DATE: 3-Mar-1987 8:01:29 

DESCRIPTION: 


This module outputs an on-off keyed signal used to drive 
a lightwave fiber module. The user specifies the on level 
via the value of parameter "avg optical power", which is 
in units of decibel milliwatts. 

The logical transmitted bitstream is also output for 
useage in an error estimator module. 

REVISIONS: 


Author : SYSADMIN 

Date : 4-Jun-1987 17:43:54 

Description: 

Module RECOMPILED by CONVERTER. 

Author : SYSADMIN 

Date : 3-Mar-1987 8:01:29 

Description: 

Module RECOMPILED by CONVERTER. 


INPUT SIGNALS: (none) 

OUTPUT SIGNALS: 

TRANSMITTED BITSTREAM Type: LOGICAL 

Lower Limit: NIL 

Upper Limit: NIL 

The (logical) transmitted bitstream. 


OPTICAL OUTPUT Type: REAL 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

Signal representing the optical power out of a digital 
lightwave source (for use with the SINGLE MODE FIBER 
(LINEAR) module). 


PARAMETERS: 

AVG OPTICAL POWER 


Type: REAL 
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Lower Limit: -300 

Upper Limit: 300 

The average optical power in decibel milliwatts. 


BIT RATE Type: REAL 

Lover Limit: 3.0e-39 

Upper Limit: 1.7e38 

Bit rate in Hertz. 


I SEED Type: INTEGER 

Lover Limit: 1 

Upper Limit: 2147483647 

Initial seed for the uniform random number generator. Should 
be a large odd integer. 


COMPUTED PARAMETERS: 

ON LEVEL Type: REAL 

Lover Limit: 0.0 

Upper Limit: 1.7e38 

Computed by init code, the optical power for the on state. 


MODULES USED IN BLOCK DIAGRAM: 

RANDOM DATA 
LOGICAL TO NUMERIC 

INITIALIZATION CODE: 

[Subroutine: ook_opt ical_source 

Arguments: 

AVG OPTICAL POWER 
ON LEVEL 

Description: 

This routine reads in the average optical power in decibel 
milliwatts, and outputs the on level necessary to produce 
this amount of optical power. 
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MODULE NAME: OOK_ERROR RATE ESTIMATOR 

GROUP: ESTIMATORS 

DATABASE: DUAO: [ BOSS. ADD_DB . LIGHT ] 

AUTHOR: TOVNSEND 

CREATION DATE: 6-Mar-1987 8:21:22 

DESCRIPTION: 

This module uses the combined analytic-simulation 
technique for estimating the probability of error for 
On-Off Keyed (00K) systems. This module is used in con- 
junction with the NOISE BW IMPULSE INJECT (REAL) module, 
which is used during calibration. 

During calibration time, the module calculates the 
equivalent noise bandwidth of the linear portion of 
system, estimates the number of samples of delay 
between the transmitted and received 
signals. (These numbers are printed out for viewing 
in the post processor (VIEW PRINTED OUTPUT)). 

The simulation stop time must be greater than or equal 
to est_stop, where est_stop is given by the following 
formula: 

est stop = caltime + 3*nbwdur + (nsseq + 2*maxdel)*dt 
~ + symtime*(numsym+l) 

where: 

caltime = "calibration start time" 
nbwdur = "noise bw calc duration" 
nsseq * "# samples in corr subseq" 
maxdel = "max delay to calculate" 
dt = simulation dt 

symtime = the symbol time in seconds 
numsym = "# symbols for error est" 

This module also counts the number of bit errors 
which occured during simulation and outputs the number 
of errors vs normalized symbol time as a multi-dimensional 
plot in the post processor. 

The module provides an additional independent axis for 
the multi-dimensional error rate plot. This is for use 
with the BOSS PARAMETER ITERATION option. 

Also, this module outputs a signal which can be used to 
drive a TERMINATE SIMULATION module when a specified number 
of errors have been counted in the simulation for ALL 
samples in a bit interval. For further details refer to 
external documentation. 

REVISIONS: 

Author : KOMP 

Date : 5-Jun-1987 17:17:11 

Description: 
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Module RECOMPILED by CONVERTER. 

Author : SYSADMIN 

Date : 6-Mar-1987 8:21:22 

Description: 

Module RECOMPILED by CONVERTER. 


INPUT SIGNALS: 

TRANSMITTED BITSTREAM Type: LOGICAL 

Lower Limit: NIL 

Upper Limit: NIL 

The (logical) transmitted bitstream input. 


RECEIVED SIGNAL Type: REAL 

Lower Limit: -1.7e38 

Upper Limit: 1.7e38 

The received (and degraded) on-off keyed signal. 


OUTPUT SIGNALS: 

TERMINATE SIMULATION Type: LOGICAL 

Lower Limit: NIL 

Upper Limit: NIL 

This signal becomes .true, when all samples in the symbol 
have exceeded the user input "max # of errors". This 
signal is typically connected to a TERMINATE SIMULATION 
module. 


PARAMETERS: 

CALIBRATION START TIME Type: REAL 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

Start time (absolute) of error rate estimator calibration 
in seconds. Nothing is happening in the error estimator 
before this time. 


NOISE BV CALC DURATION Type: REAL 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

Duration of the noise bandwidth calculation time interval 
in seconds. 
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# OF SAMPLES IN CORR SUB SEQ Type: INTEGER 

Lover Limit: 1 

Upper Limit: 4096 

This parameter specifies the number of samples used to 
perform the correlation between the 2 input data streams. 
The trade-off is accuracy of the delay estimate vs the 
amount of time needed to perform the correlation. 


MAX DELAY TO CALCULATE Type: INTEGER 

Lover Limit: 1 

Upper Limit: 2147483647 

The maximum delay which the module can calculate. (Units 
are in samples). 


# OF SYMBOLS FOR ERROR EST Type: INTEGER 

Lover Limit: 1 

Upper Limit: 2147483647 

This is the number of symbols which are used in the actual 
calculation of the probability of error. 


# OF SAMPLES/SYMBOL Type: INTEGER 

Lover Limit: 1 

Upper Limit: 10000 

The number of samples in a symbol interval. 


THERMAL NOISE PSD Type: REAL 

Lover Limit: 3.0e-39 

Upper Limit: 1.7e38 

The power spectral density (SINGLE SIDED) of the 
analytically handled thermal noise. 


DECISION THRESHOLD Type: REAL 

Lover Limit: 0.0 

Upper Limit: 1.7e38 

The decision threshold used to determine whether a 1 or 
zero was sent. 


MINIMUM # OF ERRORS Type: REAL 

Lover Limit: 0.0 

Upper Limit: 1.7e38 

The minimum number of errors which must occur in each 
sample of the symbol interval before the "terminate 
simulation" signal toggles .true. 
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PROB OF ERROR AXIS LABEL Type: CHARACTER 

Lower Limit: 1 

Upper Limit: 255 

Label which will appear on the dependent axis of the plot. 
This parameter can also be thought of as the title of the 
plot, and must be unique for multiple instances of 
multi-dimensional plots in the same simulation. 


INDEPENDENT AXIS #2 LABEL Type: CHARACTER 

Lower Limit: 1 

Upper Limit: 255 

(If used) label which will appear on the independent axis 
of the multi-dimensional error rate plot corresponding to 
the independent axis #2 values. This parameter and the 
parameter "independent axis #2 value" provide an additional 
independent axis for the error rate plot in the post 
processor, and are typeically used with the BOSS 
PARAMETER ITERATION option. If not used, the parameter 
should be set to the default value *N0T USED*. 


INDEPENDENT AXIS #2 VALUE Type: REAL 

Lover Limit: -1.7e38 

Upper Limit: 1.7e38 

Value which will appear on the independent axis 
corresponding to the independent axis #2 label. The value 
entered for this parameter is ignored if the independent 
axis #2 label has the value *N0T USED*. This parameter 
and the parameter "independent axis #2 label" provide 
an additional independent axis for the error rate plot 
in the post processor, and are typically used with the 
BOSS PARAMETER ITERATION option. 


COMPUTED PARAMETERS: 

SIM STOP TIME Type: REAL 

Lover Limit: 0.0 

Upper Limit: 1.7e38 

The simulation stop time (computed by init code). 


DELAY EST START TIME Type: REAL 

Lover Limit: 0.0 

Upper Limit: 1.7e38 

Absolute time (seconds) in which the estimation of 
delay begins. 
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NOISE BV START TIME Type: REAL 

Lover Limit: 0.0 

Upper Limit: 1.0e37 

Time (seconds) in which noise bandwidth calculation begins. 
Computed by init code. 


NOISE BU STOP TIME Type: REAL 

Lower Limit: 0.0 

Upper Limit: 1.0e37 

Time (seconds) after which the noise bandwidth calculation 
ends. Computed by init code. 

ERROR EST START TIME Type: REAL 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

Time (seconds) after which actual calculation of the 
probability of error begins. Computed by init code. 


MODULES USED IN BLOCK DIAGRAM: 
ERROR ENABLE_& SYMBOL COUNTER 
SYMBOL SAMPLE- NUMBER 
T_>= T_0N &_T_<= T_0FF 
NOI SE_BW COMPUTE (REAL) 

DELAY METER (REAL) 

OOK PRINT 
DELAY TX_SIG 
SAMPLE &_H0LD 
NOT 

PRINT SIGNAL 

T _ >= T _ 0N 
OR 

LOGICAL TO NUMERIC 

COUNT ERR0RS_& GENERATE STOP 

OOK ERROR PROB CALCULATE 

INITIALIZATION CODE: 

Subroutine: ook_error_est 

Arguments: 

CALIBRATION START TIME 
NOISE BU CALC DURATION 

# OF SAMPLES IN CORR SUB SEQ 
MAX DELAY TO CALCULATE 

# OF SYMBOLS FOR ERROR EST 

# OF SAMPLES/ SYMBOL 
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NOISE BU START TIME 
NOISE BW STOP TIME 
DELAY EST START TIME 
ERROR EST START TIME 
SIM STOP TIME 

Description: 

This is rather complicated. The routine computes the 
various times required by modules in the simulation from 
more user friendly parameters. Many of the parameters this 
routine calculates are written into the PRINTED DATA option 
in the post processor. Refer to external documentation 
for further details. 
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MODULE NAME: PIN PHOTO-DETECTOR 

GROUP: NOISE AND INTERFERENCE 

DATABASE : DUAO : ( BOSS . ADD_DB . LIGHT ] 

AUTHOR: TOWNSEND 

CREATION DATE: 19-Feb-1987 14:17:40 

DESCRIPTION: 

This module implements a PIN diode photo detector. The 
input signal should be an optical power signal (Watts). 
Any small amount of negative power on this input will be 
ignored by the module (this may occur in an FFT filter, 
for example). Shot noise effects are simulated using 
a Poisson random number generator. The output current 
is in units of microamps. 

REVISIONS: 

Author : SYSADMIN 

Date : 6-Jun-1987 1:53:33 

Description: 

Module RECOMPILED by CONVERTER. 

Author : SYSADMIN 

Date : 19-Feb-1987 14:17:40 

Description: 

Module RECOMPILED by CONVERTER. 


INPUT SIGNALS: 

OPTICAL POWER IN Type: REAL 

Lower Limit: -1.7e38 

Upper Limit: 1.7e38 

The incident optical power in watts. 


OUTPUT SIGNALS: 

OUTPUT CURRENT Type: REAL 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

The pin diode output current in MICROAMPS. 


PARAMETERS: 

DARK CURRENT (NANOAMPS) Type: REAL 

Lower Limit: 0.0 
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Upper Limit: 1.7e38 

The PIN dark current in nanoamps. 


INITIAL SEED Type: 

Lower Limit: 1 

Upper Limit: 2147483647 

Initial seed for Poisson random number generator 
be large and odd. 

QUANTUM EFFICIENCY Type: 

Lower Limit: 0.0 

Upper Limit: 1.0 

The quantum efficiency of the pin diode. 

SOURCE CENTER WAVELENGTH Type: 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

Center wavelength of the source in micrometers, 
used to determine the optical frequency. 


COMPUTED PARAMETERS: 

ELECTRON CONST Type: 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

Electron constant, computed. 

DARK CONST Type: 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

Dark constant, computed. 

POWER CONST TyP e: 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

Power constant, computed. 


MODULES USED IN BLOCK DIAGRAM: 
SELECT 
>= REAL 
ADDER 

MULTIPLIER 


INTEGER 

Should 

REAL 

REAL 
This is 

REAL 

REAL 

REAL 
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CONST GEN 
POISSON RAN_GEN 

INITIALIZATION CODE: 

Subroutine: pin_photo_detector 

Arguments : 

QUANTUM EFFICIENCY 
SOURCE CENTER WAVELENGTH 
DARK CURRENT (NANOAMPS) 

POWER CONST 
DARK CONST 
ELECTRON CONST 

Description: 

This routine calculates 3 constants: "power const", 

"dark const", and "electron const". "Power const" equals 
the product of "source center wavelength", "quantum 
efficiency", and "dt", divided by 1.989e-13. "Power 
const" converts the input optical power units to average 
number of charge carriers. "Dark const" is the number 
of dark charge carriers in delta t seconds. "Electron 
const" converts charge into current in microamps. 
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MODULE NAME: SINGLE MODE FIBER (LINEAR) 

GROUP: CHANNELS * 

DATABASE: DUAO : [ BOSS . ADD_DB . LIGHT ] 

AUTHOR: TOWNSEND 

CREATION DATE: 10-Feb-1987 10:33:43 

DESCRIPTION: 

This module implements one model of a single mode fiber. 

The initialization code reads in the parameters which 
describe the fiber and source characteristics, and 
generates from these the transfer function H(f) of the 
fiber. The init code then takes the inverse FFT of H(f), 
and truncates the insignificant portions from the resulting 
impulse response. This impulse response is written to 
a save file so that in subsequent runs it doesn't have to 
be recalculated (this is controlled by a parameter). 

The model is linear since it relys of FFT relationships. 

The model is discussed in "IEEE Journal on Selected Areas 
in Communications", January 1984, VOL SAC-2, NO 1, on 
pages 171 thru 184, by Don G. Duff. The details of this 
model are to complicated to describe here, please refer 
to external documentation for a more complete description 
of this module. 

REVISIONS: 

Author : SYSADMIN 

Date : 4-Jun-1987 17:44:54 

Description: 

Module RECOMPILED by CONVERTER. 

Author : SYSADMIN 

Date : 10-Feb-1987 10:33:43 

Description: 

Module RECOMPILED by CONVERTER. 


INPUT SIGNALS: 

JNPUT Type: REAL 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

Fiber input signal. This signal is proportional to the 
optical power of the fiber input. 


OUTPUT SIGNALS: 
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OUTPUT Type: REAL 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

Fiber output signal (optical power signal). 


PAD Q * 

DISPERSION ADJUSTMENT CONSTANT 
Lower Limit: 3.0e-39 

Upper Limit: 1.7e38 


Type : REAL 


When the internal closed form approximation for dispersion 
is used, this constant is multiplied by the dispersion 
function. A typical value is 0.0412. 


Type: REAL 


RMS SOURCE SPECTRAL WIDTH 
Lower Limit: 3.0e-39 

Upper Limit: 1.7e38 

The rms source spectral envelope width (standard deviation) 
of the closed form internal Gaussian approximation to 
source spectral envelope. This parameter is ignored if 
tabular source spectral envelope is used. Units are 
micrometers . 


Type : REAL 


SOURCE CENTER WAVELENGTH 
Lower Limit: 3.0e-39 

Upper Limit: 1.7e38 

The center wavelength of the source (micrometers). This 
parameter is used only when the closed fo ™ Gaussian 
approximation to the source spectral envelope is used. 

The parameter is ignored in all other cases. 

ZERC>”DISPERSION WAVELENGTH Type: REAL 

Lower Limit: 3.0e-39 

Upper Limit: 1.7e38 

The fiber zero dispersion wavelength (micrometers). This 
value is used only when the "mode of operation is set 
to require closed form internal approximation. of 
dispersion, and is otherwise ignored. 

GROUP DELAY POINT (DELAY) T YP e: REAL 

Lower Limit: -1.7e38 

Upper Limit: 1.7e38 

The group delay value (picoseconds/Kilometer) the R ° int 

which the tabular group delay finction is forced to pass 
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thru. (This function is obtained by integrating the 
tabular dispersion). This value is not used if the 
internal closed form function for dispersion is used. In 
this case, the group delay equals zero at the zero 
dispersion wavelength. 


GROUP DELAY POINT (WAVELENGTH) Type: REAL 

Lower Limit: 3.0e-39 

Upper Limit: 1.7e38 

The wavelength value (micrometers) for the point which the 
tabular group delay function is forced to pass thru. (This 
function is obtained by integrating the tabular 
dispersion). This value is not used if the internal closed 
form function for dispersion is used. In this case, the 
delay is zero at the zero dispersion wavelength, 
form function for dispersion is used. 


FIBER LENGTH Type: REAL 

Lower Limit: 0.0 

Upper Limit: 1.7e38 

The length of the fiber in Kilometers. 


# OF POINTS IN FIBER H(F) Type: INTEGER 

Lower Limit: 2 

Upper Limit: 10000 

The number of points in the fiber transfer function. This 
number sets the resolution of the transfer function and 
specifies the number of times that the integration over 
wavelength is performed. Must be EVEN. 


# OF PARTITIONS FOR INTEGRATION Type: INTEGER 

Lower Limit: 4 

Upper Limit: 10000 

The number of partitions used in the Simpson's rule 
numerical integration over wavelength. There are several 
considerations in choosing this number. First, it must 
be an EVEN number. The larger this number, the more 
accurate the integration but the longer it will take to 
generate the fiber impulse response. Also, if this number 
is no chosen large enough, the integrand will not be 
sampled often enough and errors due to undersampling will 
result. Many times these errors will show up as a 
transfer function with many large bumps (it won't look 
anything like a Gaussian shaped spectrum). 


INTEGRATION UPPER LIMIT Type: REAL 
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Lover Limit: 3.0e-39 

Upper Limit: 1.7e38 

Upper limit of integration which generates a sample of the 
fiber transfer function. Units are in micrometers. 

See external documentation for further info. 


INTEGRATION LOWER LIMIT Type: REAL 

Lover Limit: 3.0e-39 

Upper Limit: 1.7e38 

Lover limit of integration which generates one sample of 
the fiber transfer function. The units of this value are 
micrometers (this integration is with respect to wave- 
length). See external documentation for further info. 


FIBER H(T) SAVE FILENAME Type: NEV-FILE 

Lover Limit: NIL 

Upper Limit: NIL 

Filename of file containing the impulse response of the 
fiber. A new version of this file is written out each 
time the parameter "recalculate fiber h(t)" is .true. 

When this parameter is set to .false., the module READS 
the impulse response from this file. 

This file is required because calculation of the 
fiber impulse response can require a reletively large 
amount of cpu time. This file allows the user to generate 
an impulse response for a given set of fiber parameters, 
then for subsequent simulation runs, read from this file 
(assuming the fiber parameters are not changed in these 
subsequent runs). 


DISPERSION FILENAME Type: OLD-FILE 

Lover Limit: NIL 

Upper Limit: NIL 

Filename of file containing a tabular specification of 
fiber dispersion (picoseconds/kilometer/micrometer) vs 
wavelength (micrometers). This file must be in the BOSS 
preferred file format. Data from this file is only used 
if the mode is either CLOTAB or TABTAB, but a valid old 
filename must be entered in any case. This function must 
span the range of integration specified by the parameters 
"integration lover limit" and "integration upper limit". 

The antiderivi tive of this function is calculated to 
get the group delay function for the fiber. This group 
delay function is then, raised to the power e and multiplied 
by the fiber length in kilometers. 

The resulting group delay function is required to pass 
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thru the point specified by the parameters "group delay 
point (wavelength)" and "group delay point (delay)" (this 
evaluates the constant of integration). 

The maximum number of entries is 10,000, 


SOURCE SPECTRAL ENVELOPE FILENAME Type: OLD-FILE 

Lover Limit: NIL 

Upper Limit: NIL 

Filename of file containing the tabular specification 
of source spectral envelope vs wavelength (micrometers). 
For the overall gain to be correct, the source spectral 
envelope function MUST be normalized so that the area 
under this function equals unity. The BOSS preferred file 
format must be used, and wavelength values must be in the 
first column, with source spectral envelope values in the 
second. This function does not have to span the range 
of integration specified by the parameters "integration 
lover limit" and "integration upper limit", but the 
function will be assumed to be zero outside the values 
given. Since this function is a factor in the integrand, 
it should approach zero within the range of integration. 

This file is not required if the mode is either 
CL0CL0 or CLOTAB, BUT a valid old filename must still be 
entered in any case (the fiber loss filename could be used 
in these cases). 

The maximum number of entries is 10,000. 


FIBER LOSS FILENAME Type: OLD-FILE 

Lover Limit: NIL 

Upper Limit: NIL 

Name of file containing the tabular values of fiber loss 
(dB/Kilometer ) vs wavelength (micrometers). (Note that the 
loss is 1/gain). The file MUST be in the BOSS prefered file 
format, with wavelength in the first column and loss in the 
second column. This file is required in ALL modes of 
operation. The minimum and maximum values of wavelength 
MUST span the range of integration specified by the 
parameters "integration lover limit" and "integration 
upper limit". 

The maximum number of entries is 10,000. 


RECALCULATE FIBER H(T) Type: LOGICAL 

Lover Limit: NIL 

Upper Limit: NIL 

If this logical variable is .true., the current parameter 
values are used to recalculate a new fiber impulse 
response and write it to the save file. If .false., the 
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impulse response is read from the save file. This 
parameter is needed because the calculation of the fiber 
impulse response takes a long time (up to 1 cpu hour). 

So the first time a simulation is run with new parameter 
values, this parameter should be set to .true. For all 
subsequent runs using the same parameters, the value 
should be .false. THEREFORE, THIS PARAMETER SHOULD 
BE EXPORTED TO THE HIGHEST (SIMULATION) LEVEL! 

MODE OF OPERATION Type: CHARACTER 

Lower Limit: 1 

Upper Limit: 6 

This parameter specifies how the module will compute the 
fiber impulse response. There are 4 legal character 
entries (case INsensitive) : 

cloclo - Closed form internal approximation for source 
spectral envelope and dispersion, 
clotab - Closed form internal approximation for source 
spectral envelope and tabular dispersion data 
supplied from a file. 

tabclo - Tabular source spectral envelope data from a file 
and closed form internal approximation for 
dispersion. 

tabtab - Both source spectral envelope AND dispersion are 
supplied as tabular inputs from 2 separate files. 
For further information about the tabular data files 
or the approximations see the respective help info. 


COMPUTED PARAMETERS: 

TAP GAIN Type: REAL 

Lower Limit: -1.7e38 

Upper Limit: 1.7e38 

Vector Length: (# OF TAPS) 

Tap gain for this cell. 

# OF TAPS Type: INTEGER 

Lover Limit: 1 

Upper Limit: 4096 

The number of taps in the fiber impulse response. 
(Computed by initialization code). 


MODULES USED IN BLOCK DIAGRAM: 
TAPPED DELAY LINE CELL (REAL) 
SINK 

CONST GEN 
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INITIALIZATION CODE: 


Subroutine: single_mode_f iber 

Arguments: 

MODE OF OPERATION 
RECALCULATE FIBER H(T) 

FIBER LOSS FILENAME 
SOURCE SPECTRAL ENVELOPE FILENAME 
DISPERSION FILENAME 
FIBER H(T) SAVE FILENAME 
INTEGRATION LOVER LIMIT 
INTEGRATION UPPER LIMIT 

# OF PARTITIONS FOR INTEGRATION 

# OF POINTS IN FIBER H(F) 

ZERO DISPERSION WAVELENGTH 
FIBER LENGTH 

SOURCE CENTER WAVELENGTH 
RMS SOURCE SPECTRAL WIDTH 
DISPERSION ADJUSTMENT CONSTANT 
GROUP DELAY POINT (DELAY) 

GROUP DELAY POINT (WAVELENGTH) 

# OF TAPS 
TAP GAIN 

Description: 


Most of the work in this module is performed here. The 
routine reads in all parameters required to calculate the 
fiber impulse response. Description of the module 
describes this routine, so we will not be redundant here. 
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LASER DIODE MODELS 


Fiber optic system sources and detectors are generally composed of direct 
bandgap semiconductor materials such as Gallium Arsenic which have an energy 
gap approximately equal to the energy in a photon of light of the desired 
wavelength. In photodetectors, most received photons cause ionization of a 
lattice atom creating a hole-electron pair (which can then create more hole- 
electron pairs in an avalanche photodiode). In optical sources, a forward 
biased junction is used to inject minority carriers, which create photons when 
they recombine. 

Optical sources are of two types; light emitting diodes (LED ’ s) and laser 
diodes. In LED’s, recombination of injected minority carriers causes emission 
of light with a relatively broad wavelength range. The optical power emitted 
is approximately proportional to the LED current, but relatively little power 
may be coupled into an optical fiber. Also, the frequency response of LED’s 
limit the frequencies at which they can be modulated. 

Laser diodes have a better frequency response and can be modulated at 
frequencies up to several gigahertz, and more power can be coupled into a 
fiber over a narrower wavelength range. These properties make lasers 
preferable to LED’s for applications which require high power or high 
frequency modulation. 

However, lasers also have disadvantages, one of which is their nonlinear- 
ity. At low Input currents, laser diodes act essentially as LED’s, producing 
light by spontaneous emission. When the current is Increased to a value 
called the threshold current, stimulated emission (lasing) begins and causes 
the power output to increase much faster with increasing current than it did 
below threshold. In stimulated emission, a photon striking a lattice atom in 
the active region can either be absorbed In ionizing the atom, or can cause a 
photon of the same wavelength to be released. Conditions are created so that 
the latter occurs more frequently, so that the light is amplified. Reflecting 
facets are placed at the ends of the active region so that light may continue 
to be amplified. It is this positive feedback which causes the power to be 
increased over that due to spontaneous emission. The initial photons to start 
this process come from spontaneous emission. 


The power output versus current input curve (the laser’s L-I characteris- 
tic) rises sharply above threshold due to the higher power output with stimu- 
lated emission than with spontaneous emission alone* This makes the laser 
diode an inherently nonlinear device, and therefore it is essential in analog 
applications that the laser remain above threshold for all values of the 
input. The laser is thus biased with a current which is about halfway between 
the threshold current and the highest laser current (in the middle of the 
lasing region). Because of the steepness of the L-I characteristic above 
threshold, the bias current is generally a small fraction above the threshold 
current. This makes temperature stability of the laser crucial, since the 
threshold current varies with temperature. 


Since lasers are generally kept above threshold, it is the output in this 
region that it is most important to model. In this region, the slope of the 
L-I characteristic for all early and many current laser diodes exhibit promi- 
nent nonlinearities called "kinks," making them unsuitable for analog use. 
These kinks have several proposed explanations, such as the excitation of 
higher order transverse modes [1] (the active region of a laser diode is a 
resonant cavity). Newer lasers, such as narrow stripe lasers, have decreased 
the active region width so that these modes do not occur at normal current 
levels, and so these lasers can have a very linear characteristic above thres- 
hold. ^ Third order intermodulation signals (the most important if signal 
frequencies are kept within an octave) can be as much as 60 to 80 dB below the 
fundamental {2], so for some lasers a linear model (requiring only the bias 
point and slope of the L-I characteristic at the bias point) may be fairly 
accurate. 


If the L-I relationship is not considered to be linear, the most 
straightforward way to model it is as a polynomial 


L * al + bl^ + cl 3 + • • • 


where L is the power output and I is the current input relative to their bias 
values. The coefficient of the I term (a) is the slope of the L-I character- 
istic at the bias point (the efficiency), and the other coefficients may be 
determined from measured harmonic magnitudes, along with the input currents 
for which they were measured. 
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Nonlinearities are also sometimes described by measured values of AM to 
AM and AM to PM characteristics, which indicate the effect of variations in 
the input signal amplitude on the output magnitude and phase, respectively. 
Measurements on some lasers have shown [3] that below resonance, most of the 
intermodulation terms are due to AM to AM characteristics, while at the la 
ser T s resonant frequency, AM to PM characteristics may become significant due 
to increased reflection noise. 

There are many models based on the physics of laser diodes , some of which 
are reviewed by Buus [A] . Some of these [5, 6, 7] allow the calculation of 
points on the L-I characteristic, which can then be used directly, or to 
determine the coefficients of a polynomial approximation. However, these 
models require values of many laser parameters, many of which would not be 
accurately known, and they require iterative solution. Although these models 
agree fairly well with measured threshold currents and efficiencies, the 
degree to which they model the nonlinearity is unknown. 

There is another source of nonlinearity which may be as significant as 
the nonlinearity of the L— I characteristic. Most laser diodes do not emit a 
single wavelength of light, but rather a narrow spectrum of discrete wave- 
lengths. The active regions of many lasers (such as stripe geometry lasers) 
have one dimension (the length) which is much greater than the other two. The 
active region is a resonant cavity, and the output wavelengths correspond to 
different longitudinal modes (normally, only the lowest order transverse modes 
are permitted). The spacing between the output wavelengths thus depends on 
the cavity length, and the relative magnitudes depend on the spontaneous 
emission and the optical gain of the different wavelengths. The shape of the 
output spectrum is generally modeled as Gaussian or Lorentzian [8] , but the 
center wavelength and width of the spectrum are both functions of the output 
power level. 

The shift in the center wavelength, when combined with a fiber with 
dispersion (different wavelengths propagate at different speeds) give rise to 
harmonics for a sinusoidal input current, and thus is a source of nonlinearity 
[9, 10]. This shift in center frequency may be on the order of the spacing 

between modes (or less) and still give rise to a possibly significant effect 
[9]. Published spectra [6, 11, 12] sometimes show a greater shift than 
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this. The width of the output spectrum has been shown to be approximately 
inversely proportional to the output power, so that as the output power in- 
creases, the spectrum shifts to longer wavelengths and narrows [6] • 

The model of the laser output spectrum may also need to account for modal 
partition noise. Although for a given constant input current, the total power 
output is nearly constant, the power in each mode may vary greatly. If the 
light then travels through a dispersive fiber, different wavelengths will 
travel at different velocities, so the power reaching the receiver will vary 
with time. Since the partitioning is random, this will have the effect of 
producing noise (unlike the shift in wavelengths above, which occurs synchron- 
ously with modulation and appears as a nonlinearity). Given the center wave- 
length, modal spacing, a measure of the spectral width and the total output 
power, the time averaged spectrum can be calculated using a Gaussian or 
Lorentzian envelope. Then, a random component can be added to each mode since 
the exact partitioning of the output power is random. Since the total output 
power is nearly constant, the power in the different modes must be negatively 
correlated, and the covariance between modes can be approximated by [8] 

C(i, j) = -K 2 A(i) A(j) 

C(i> i) “ K 2 A(i) [1 - A(i) ] 

where C(i, j) is the covariance between modes i and j, C(i, i) is the variance 
of the power in mode i, A(i) is the fraction of the total power in mode i, and 
K is a partition coefficient which has been measured and determined to be 
between 0.3 or 0.4 and 0.6 [8, 13]. 

So a complete model of the output spectrum would require that first the 
total output power be determined based on the input current and information 
about the L-I characteristic, and from this average power the spectral width 
and center wavelength could be calculated. From this and the mode spacing, 
the number of modes with significant power could be determined along with the 
wavelength and average power for each. Then, a random component with the 
above covariances could be added to each. Since modal partition noise is more 
of a problem when there are few modes, it may be reasonable to ignore this 
factor when there are many modes since the effect Is reduced while the calcu- 
lations required are increased. 
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This would be required for each sample of the input current, which may 
make the model impractical to use without further simplification- The compli- 
cation is that the nonlinearity associated with the spectral shift and the 
modal partition noise are both effects due to both the laser source and dis- 
persion in the optical fiber, so to model the laser alone so that when com- 
bined with any fiber model these effects will be observed seems to require 
calculating an output spectrum for each current level. 

The discussion so far has assumed that the laser is memoryless, so that 
the output at any time depends only on the input at that time. The measured 
frequency response of laser diodes is generally flat a low frequencies, then 
has a resonant peak, after which the power output decreases with increasing 
frequency at 40 dB/decade[2] . As long as the modulation frequency is in the 
flat region of the frequency response (below the resonant frequency) this is a 
reasonable assumption. If it is modulated above this region, then the fre- 
quency response might be inserted as a filter placed after the input current 
but before the output power is calculated. However, the frequency response is 
often flat up to several gigahertz, and it is generally desirable to use the 
laser in this region, so the frequency response can probably be neglected most 
of the time . 


- 59 - 


REFERENCES 


1. Shen, T. M., ’’Effects of TE-TM Kinks on the Performance of Semiconductor 
Lasers in Lightwave Communication Systems,” Journal of Lightwave Technol- 
ogy, Vol. LT-4, No. 9, September 1986, pp. 1420-1424. 

2. Channin, D. J., ’’Optoelectronic Performance Issues in Fiber-Optic Communi- 
cations,” RCA Review , Vol. 46, December 1985. 

3. Way, W. I., and A. Afrashteh, ’’Linearity Characterization of Connectorized 
Laser Diodes Under Microwave Intensity Modulation by AM/AM and AM/PM 
Measurements , ” 1986 IEEE MTT-S International Microwave Symposium Digest , 
pp. 659-662. 

4. Buus, J., ’’Principles of Semiconductor Laser Modelling,” IEE Proceedings , 

Vol. 132, Pt. J., No. 1, February 1985, pp. 42-51. 

5. Streifer, W. D. , R. Scifres and R. D. Burnham, ’’Analysis of Diode Laser 

Properties,” IEEE Journal of Quantum Electronics , Vol. QE-18, No. 11, 
November 1982, pp. 1918-1929. 

6. Streifer, W. D., R. Scifres and R. D. Burnham, ’’Longitudinal Mode Spectra 
of Diode Lasers,” Applied Physics Letters, Vol. 40, No. 4, February 15, 
1982, pp. 305-307. 

7. Streifer, W. D., R. Scifres and R. D. Burnham, "Analysis of Spontaneous 

Emission Effects on Spectra and L v$. I Characteristics of Diode Lasers,” 

Japanese Journal of Applied Physics, Vol. 21, No. 5, May 1982 , pp. L282- 
L284 . 

8. Laughton, A. K. , "Mode Partition Noise in Gain-Guided Lasers and Its 
Effect on a Multimode Fibre-Optic System,” IEE Proceeding s, Vol. 132, Pt. 
J., No. 6, December 1985, pp. 359-363. 

9. Lyons, P. B., S. D. Personick and G. A. Hayward, "Applications of Optical 

Fibers to Analog Telemetry Delay Lines and Sensing Systems," IEEE Journal 

on Selected Areas in Communications, Vol. SAC-1, No. 3, April 1983, pp. 
555-561. 

10. Personick, S. D. , Fiber Optics , Plenum Press, 1985, pp. 123-125. 

11. Wada, M. , K. Hamada, T. Shibutani, H. Shimizu, M. Kums, K. Itoh, G. Kano 

and I. Teramoto, "A New Chemical Etching Technique for Formation of Cavity 
Facets of (GaAl) As Lasers," IEEE Journal of Quantum Electronics, Vol. QE- 
21, No. 6, June 1985, pp. 658-662“. " 

12. Burkhard, H. , and E. Kuphal, "Three- and Four-Layer LPE InGaAs(P) Mushroom 

Stripe Lasers for A * 1.38, 1.54 and 1.66 pm," IEEE Journal of Quantum 

Electronics , Vol. QE-21, No. 6, June 1985, pp. 650-657. 

13. Duff, D. G., "Computer-Aided Design of Digital Lightwave Systems," IEEE 
Journal on Selected Topics in Communications, Vol. SAC-2, No. 1, January 
1984, pp. 171-185. 


- 60 - 



BIBLIOGRAPHY 


Multiple Topics 


The following books and articles each discuss most of the important 
characteristics of laser diodes, such as their nonlinearity, frequency re- 
sponse, modal noise and mode partition noise. Duff also discusses modeling of 
many laser characteristics that show up in digital applications, such as turn- 
on and turn-off response, pattern dependence and oscillations (also discussed 
by Channin) . The linearity of some current laser diodes are quite good 
(Channin) and the frequency response of some allows them to be modulated above 
1 GHz. 

!• Channin, D. J., "Optoelectronic Performance Issues in Fiber-Optic Communi- 
cations, RCA Review , Fol. 46, December 1985. 

2. Daly, J. C., (Ed.), Fiber Optics , CRC Press, 1984. 

3. Duff, D. G., "Computer Aided Design of Digital Lightwave Systems," IEEE 
Journal on Selected Topics in Communications , Vol. SAC-2, No. 1, January" 
1984, pp. 171-185. 

4. Gowar, J., Optical Communications Systems , Prentice-Hall International 
Series in Optoelectronics, 1984 . 

5. Personick, S. D. , Fiber Optics , Plenum Press, 1985. 

Quantum Physical Models 


The following articles and books are concerned with modeling the physical 
operation of laser diodes. Buus is a review of the relevant equations used 
and of the various published models based on these. The articles by Streifer, 
Scifres and Burnham give a model which can be used to calculate the L-I char- 
acteristic of a laser diode given a number of its physical parameters. Many 
of these would be hard to determine, the model requires an iterative solution, 
and the results do not necessarily accurately describe the characteristics we 
are interested In. However, these articles also model the spectral width, 
which appears to be Inversely proportional to the output power, which may be 
useful in modeling the output spectrum. 


- 61 - 



6 . 


Buus, J., "Principles of Semiconductor Laser Modeling," IEE Proceedings 
Vol. 142, Pt. J., No. 1, February 1985, pp. 42-51. 

7. Streifer, W. D. , R. Scifres and R. D. Burnham, "Analysis of Diode Laser 

Properties, IEEE Journal of Quantum Electro nics. Vol. QE-18 No. 11 
November 1982, pp. 1918-1929'.' ** ’ “ ’ 

8. Streifer, W. D.,^R. Scifres and R. D. Burnham, "Longitudinal Mode Spectra 

of Diode Lasers," Applied Physics Letters, Vol. 40. No. 4. February 15 
1982, pp. 305-307. * 3 

Streifer, W. D., R. Scifres and R. D. Burnham, "Analysis of Spontaneous 
Emission Effects on Spectra and L vs. I Characteristics of Diode Lasers,” 
Japanese Journal of Applied Physics . Vol. 21, No. 5, May 1982, pp. L282- 


10. Thompson, G. H. B. , Physics of Semiconductor Laser Devices, Wiley, 1980. 


Nonlinearity 

These articles are concerned with the nonlinearities of laser diodes. 
Shen discusses the kinks found in the above threshold region of many laser 
diodes and the effects of these on digital communication systems. Way and 
Afrashteh measured the AM/AM and AM/PM characteristics of laser diodes for 
modulation frequencies above 1 GHz, and conclude that below the resonant 
frequency, the AM/AM characteristic is most important. Lyons, Personick and 
Hayward discuss a nonlinearity due to the variation of the center wavelength 
synchronously with modulation and dispersion in the optical fiber. This is 
also discussed in Reference 5 above* 


11. Lyons, P. B., S. D. Personick, and G. A. Hayward, "Applications of Optical 

® r ® ^ nalog Telemetry Delay Lines and Sensing Systems,". IEEE Journal 
on Sel ected Areas in Communications , Vol. SAC-1, No. 3, April 1983 , ppT 
555-561* 

12. Shen, T. M. , "Effects of TE-TM Kinks on the performance of Semiconductor 

Lasers in Lightwave Communication Systems," Journal of Lightwave Technol- 
ogy^, Vol. LT-4 , No. 9, September 1986, pp. 1420-1424"! """ 

13. Way, W. I., and A. Afrashteh, "Linearity Characterization of Connectorized 

Laser Diodes ^ Under Microwave Intensity Modulation by AM/AM and AM/PM 
Measurements, 1986 IEEE MTT-S International Micr owave Symposium Digest, 
pp. 659-662. ■ ‘ — — -*® * 


- 62 - 



Some Published L-I Characteristics and Spectra 


Although there seems to be little published quantitative data on the 
nonlinearity of laser diodes or on the magnitude of shift of the center wave- 
length described above, there are very many published L-I characteristics and 
output spectra at different power levels which show the shift in the center 
wavelength (and the narrowing described in an earlier section). Below are 
listed a few articles containing one or both of these types of data. 


14. Burkhard, H., and E. Kuphal, "Three- and Four-Layer LPE InGaAs(P) Mushroom 

Stripe Lasers for X = 1.38, 1.54 and 1.66 ym,” IEEE Journal of Quantum 

Electronics , Vol. QE-21, No. 6, June 1985, pp. 650-657. 

15. Dolginov, L. M., A. E. Drakin, P. G. Eliseer, B. N. Sverdlov and E. G. 
Shevchenko, ”CW InGaAsP/InP Injection Lasers with Very Low Threshold 
Current Density at Room Temperature, ” IEEE Journal of Quantum Electronics , 
Vol. QE-21, No. 6, June 1985, pp. 646-649. 

16. Dutta, N. D., R. B. Wilson, D. P. Wilt, P. Besoni, R. L. Brown, R. J. 

Nelson and R. W. Dixon, "Performance Comparison of InGaAsP Lasers Emitting 
at 1.3 ym and .1.55 ym for Lightwave Systems Applications,” AT&T Technical 
Journal, Vol. 64, No. 8, October 1985, pp. 1857-1884. 

17. Haupt, H. and 0. Hildebrand, "Lasers and Photodetectors in Europe,” IEEE 

Journal on Selected Areas in Communications , Vol. SAC-4, No. 4, July lU£5~, 
pp. 444-455. 

18. Wada, M. , K. Hamada, H. Shimizu, M. Kume, F. Tajiri, K. Itoh and G. Kano, 

"High Power Lasers of the Twin-Ridge-Substrate Type,” IEE Proceedings , 

Vol. 132, Pt. J.. No. 1, February 1985, pp. 3-8. 


Noise 


The power output of a laser diode is almost constant for a constant 
current input, but there can be significant noise due to two mechanisms. One 
is reflections, in which the laser’s output characteristic changes when light 
is reflected back into the laser, which is discussed in Garvey and Quinn, and 
which would probably be difficult to model since it depends on the connections 
outside the laser source. The other type of noise is mode partition noise, 
which is caused by random variations in the power output of each mode (wave- 
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probably be necessary to model the output spectrum. Laughton discusses the 
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ABSTRACT 

Because of the presence of intersymbol in- 
terference and non-Gaussian noise sources, the 
evaluation of the performance of digital light- 
wave communication links is difficult. Direct 
Monte Cafrlo type simulations are impractical 
because of the long run times required to esti- 
mate low bit error rates. In this paper, tech- 
niques that combine analysis and simulation in 
modeling and implementing the functional blocks 
of a single-mode lightwave system using the 
SYSTID simulation package are discussed. A simu- 
lation example is presented which illustrates the 
capabilities of the semi-analytic approach to 
digital single-mode lightwave system analysis. 

I. INTRODUCTION 

The rapid pace in the research and develop- 
ment of improved and economical fiber cable, 
devices, and systems technology has resulted in 
large-scale deployment of lightwave communication 
systems. Early lightwave communication systems 
utilized multimode fibers which allowed only 
limited information bandwidth and repeater spac- 
ing [1]. As technology advanced, emphasis shift- 
ed away from multi-mode lightwave systems, in 
favor of the higher bandwidths and lower losses 
available with single-mode systems. This growth 
has increased the range of potential applications 
and led to lightwave systems being used more and 
more in highly competitive commercial applica- 
tions. As a result, lightwave systems engineers 
require more detailed information on the possible 
tradeoffs available when designing or deploying a 
system. Analytical methods alone are becoming 
insufficient because of many necessary simpli- 
fying assumptions. The lightwave system designer 
la often forced to evaluate design alternatives 
by constructing prototypes? an approach which can 
be costly - especially in a research or early 
design phase* 

In a recent paper (23, Ikiff presented models 
and analysis methods for computer-aided design of 
digital lightwave systems. Duff's approach was 
based primarily on analytical approximations to 
compute performance measures such as error rates, 
eye diagrams, and power penalties. In contrast, 
the approach presented here is based on detailed 
simulation models for functional blocks in light- 
wave communication links. Computer simulation of 
digital lightwave systems is an altsrnatlvs to 


both hardware evaluation and pure analytical 
methods of analysis. An important advantage of 
computer simulation is the flexibility available 
to the designer. Simulation packages such as 
SYSTID (3) offer the user the capability to study 
the effects various parameters have on total link 
performance via bit error rate (BER) estimates, 
eye diagrams, spectral plots, time plots, etc. 
In general, computer simulation involves perform- 
ing a series of transformations on a data struc- 
ture (i.e., a signal) and analyzing the effects 
of these transformations using the performance 
measures mentioned above. The transformations 
are between and within the time and frequency 
domains, and model the behavior of cascaded or 
single components in the system being simulated. 

Three important approaches to computer simu- 
lation are as follows: a) direct or brute-force 

Monte Carlo simulation, b) modified Monte Carlo 
simulation (importance sampling), and c) semi- 
analytic Monte Carlo simulation. Digital light- 
wave systems of interest typically operate with 
bit error probabilities from to 10’^. For 

the brute- force Monte Carlo approach to be accu- 
rate, a large number of errors must be observed 
which at these low error rates require on the 
order of several million simulated bits. Gener- 
ating and processing this number of bits requires 
too much computer time to be practical. Impor- 
tance sampling and semi-analytic approaches re- 
duce the number of bits simulated to a more rea- 
sonable value. 

Modeling lightwave system functional blocks 
and using these in a semi-analytic simulation 
framework are discussed in this paper. Lightwave 
systems are ideally suited for the semi-analytic 
approach. The effects of intersymbol interfer- 
ence (XSI) introduced by the optical source and 
the fiber, in addition to the quantum noise in- 
troduced by the photodetector are simulated. 
Thermal noise (assumed to be white Gaussian) 
enters the simulation only in the receiver which 
is downstream of any non-linearities introduced 
by the optical source or photodetector. Hence 
the effects of thermal noise are handled analyti- 
cally. Typically, short VN sequences (on the 
order of hundreds of bits) can be used to evalu- 
ate the effects of the non-Gaussian degradations, 
which greatly reduce the computer time required 
as compared to the brute- force approach. 

The contents of this paper are divided into 
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four sections. Section II is a brief discussion 
of the SYSTID simulation package. Section III 
presents the models for the lightwave functional 
blocks considered here. Discussed in this sec- 
tion are the single-mode optical fiber models 
featuring versions which allow tabular input of 
fiber parameters, and a choice between time do- 
main or frequency domain signal processing. A 
PIN diode model and two BER estimator models are 
also discussed. 

Section IV presents an example lightwave 
system simulation which illustrates the capabil- 
ities of SYSTID as a tool for lightwave systems 
analysis. Results from more detailed simulation 
studies are presented in [4]. The simulated 
system is a single-mode 140 Mb/s, 1.28 m LED 
link vrith a fiber length of 15 Km. Waveform 
plots are shown at various points throughout the 
system, along with a plot of probability of bit 
€tror versus received optical power. Concluding 
remarks are presented in Section V. 

II* THE SYSTID SIMULATION PACKAGE 

The SYSTID simulation package (developed 
originally by the Hughes Aircraft Corporation) is 
a powerful tool for analyzing and designing com- 
munication systems [3]. SYSTID features a flexi- 
ble topology format with an input language which 
corresponds closely to system level block dia- 
grams. This minimizes the amount of time a user 
spends translating the block diagram of the sys- 
tem under analysis Into the simulation source 
program. Although the user can easily construct 
new blocks, the SYSTID Model Library contains 
many commonly used functional blocks such as 
signal generators, filters, modems, and source 
encoders. Results of simulations are analyzed 
using the SYSTID postprocessor RIP. RIP is menu- 
driven with a large number of commands available 
such as spectral plot3 r time plots, eye diagrams, 
histograms, correlation, covariance and others. 
To decrease the amount of time a user spends 
answering questions In a given menu, standard 
default values are provided when possible. 

The structure of SYSTID makes it easy for 
the user to construct and use models which re- 
quire tabular data input. This feature was used 
in the single-mode fiber models to allow labora- 
tory hardware measurements to be used in simula- 
ting the fiber, as discussed in Section III. 

IH # MODELING FUNCTIONAL BLOCKS IN A LIGHTWAVE 
COMMON ICATIONS LINK 


are finding increasing application with single- 
mode fibers for the shorter distance distribution 
plant (7). 

The spectrum of the optical source plays an 
important role in determining the effects ja f 
dispersion on a single-mode fiber communication 
system. When using one of the four single-mode 
fiber models discussed below, the user can either 
input an arbitrary source spectral density curve 
in tabular form, or supply the needed parameters 
to specify a built-in Gaussian approximation to 
the source spectrum. If the source spectrum were 
very narrow, it could be approximated as an im- 
pulse and the effects of the single-mode fiber 
could be handled analytically. This is not a 
realistic assumption, however, and simulation 
techniques are required to evaluate the effects 
of a non-ideal source spectrum on system perform- 
ance. 

Single-Mode Fiber 

The baseband transfer function of a single- 
mode fiber can be modeled by the following equa- 
tion (8): 

H(f> - / S(X)L(X)e' :i “ I '' r( X) dX (1) 

where 

w m 2Ttf 

L “ fiber length 

S(A) * optical source spectrum vs wavelength 

L (A) - 1/(fiber loss) vs wavelength 

T(A) * fiber group delay vs wavelength 

The two major physical characteristics which 
influence the performance of the single-mode 
fiber are fiber loss and chromatic dispersion. 
In eq. (1), L( A ) is the reciprocal of loss (gain) 
as a function of wavelength. The fiber group 
delay term, T(A), is actually the antiderivative 
of the dispersion characteristic. The degrading 
effects of chromatic dispersion are related to 
the spectral purity of the optical source. Real 
optical sources contain a distribution of energy 
around the source center wavelength A , An opti- 
cal source with a broad spectrum of optical ener- 
gy results in a narrowing of the information 
bandwidth of the fiber. This can be easily noted 
frcm eq. (1) if we approximate L(X) by a con- 
stant, and T(A) by a linear function of A. Now 
the source spectrum S(A) and the fiber transfer 
function H(f) can be visualized as Fourier trans- 
form pairs (except for a scale constant). Hence, 
a broad 5(A) versus A corresponds to a narrow 
H(f) versus f, and vice versa. 


Optical Sources 

Depending upon the application, optical 
sources under consideration for lightwave commu- 
nication systems are either semiconductor light 
emitting diodes (LED) or laser diodes (5, 6). 
Both types of optical sources tend to degrade the 
input bit stream, although incoherent LEDs are 
more well behaved (and therefore easier to model) 
than semi-coherent laser sources. Lasers, how- 
ever, are capable of providing higher values of 
launched optical power and are widely used in 
long distance single-mode fiber systems. LEDs 


If the optical source spectrum is narrow 
enough to be approximated by an impulse, the 
fiber transfer function is given by 


H(f) 


L(A^) e 


-juLT(A s ) 


( 2 ) 


which results in a delayed and attenuated version 
of the input signal. Although this case is not 
interesting in Itself, this approach was used to 
verify the fiber models implemented in SYSTID. 

The source spectrum can sometimes be approx- 
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imated by a Gaussian as follows [2, 8], 


s(A) 


where 



(A - X ) 2 

~T7~' 

a 


a - nas spectral width 
A m source center wavelength 


(3) 


Chromatic dispersion for silica 
sometimes approximated by 1 2, 8], 

X - 


dT(X) 

dA 


(s/c)<- 


-) 


where 


s » a unitless constant 

c * speed of light 

Ag * zero-dispersion wavelength 


fibers is 


(4) 


The group delay is then 

x *0 - x 

T(X) - (s/c)[ln(f) + ] (5) 

A 0 A 


These closed-form approximations to the 
source spectrum and chromatic dispersion are not 
adequate when it is desired to simulate a single- 
mode fiber and optical source which have charac- 
teristics that deviate drastically from those in 
eqs. (3) and (5). In addition, no reasonable 
closed-form expression exists for fiber loss. 
Hence, some means for arbitrary entry of the 
various functions is required for a single-mode 
fiber model to be general enough to handle arbi- 
trary dispersion, loss and source spectrum char- 
acteristics. 


Four different models of the single-mode 
fiber have been implemented in SYSTID. All four 
models use eq. (1 ) to calculate the fiber trans- 
fer function. Two of the models use frequency 
domain signal processing methods, while the other 
two perform signal processing in the time domain 
via the convolution sum. One of the frequency 
and one of the time domain models use the closed 
form approximations to the source spectrum and 
chromatic dispersion curves given in eqs. (3) and 
(4) respectively. The other set of models allows 
the user to input either the source spectrum or 
dispersion, (or both) as tabular functions. The 
advantage of using the "closed-form" models is 
that they are easier to use and require less 
computer time to calculate the transfer func- 
tion. All models require fiber loss data to be 
input in tabular form. 

The integral in eq. (1) is implemented using 
Simpson's rule. Euler's identity is used to 
break the integrand into real and imaginary parts 
yielding 

H(f) - / S ( A ) L ( A ) cos ( uiLT ( A ) ) d A 

- j / S( A )L( A)sin( wLT( A) )d A (6) 

—m 

To save computer time when calculating H(f), the 
values of S(A)*L(A) and 2*ir*L-T(A) are calcu- 
lated during the first integration (i.e., the 
first value of baseband frequency, f) and stored 
in arrays for future integrations. The upper and 


lower limits of integration are parameters to the 
model, and are chosen to include the significant 
non-zero values of S(A). One additional con- 
straint on the function S(A) is that /s(A)dA - 1, 
otherwise the attenuation of the fiber will be 
incorrectly biased. 

It is important to use sufficient resolution 
in the discrete integration to provide an accu- 
rate approximation to the continuous integra- 
tion. It was found that if the resolution of the 
integrations in eq. (6) was inadequate, the fiber 
transfer function at higher frequencies would be 
incorrect. It was determined that double preci- 
sion arithmetic is required to eliminate roundoff 
errors which occur under some circumstances. 
Double precision is only required when calculat- 
ing the fiber transfer function and not when 
performing the signal processing, so the speed 
penalty for using double precision arithmetic is 
not proportional to the number of samples pro- 
cessed by the simulation. 

To help reduce the amount of computer time 
required during execution of a lightwave simula- 
tion, all four fiber models can be configured to 
either write the calculated transfer function to 
a save file, or read from this save file. This 
feature is useful when it is desirable to hold 
the single-mode fiber parameters constant and 
study the effects of varying parameters of other 
models in the simulation. 

The group delay function, T(A) is obtained 
as follows. The antiderivative of the user-input 
tabular dispersion data is automatically calcu- 
lated yielding group delay plus an arbitrary 
constant versus wavelength. Next, the arbitrary 
constant of integration is evaluated by biasing 
the group delay function to force it to pass 
through a user-specified point. There is no 
restriction on what type of function the user 
specifies for dispersion, loss, or source spec- 
tral density, as long as adequate resolution is 
used. 


The calculated transfer function of the 
single-mode fiber is used directly in the two 
models which perform the signal processing in the 
frequency domain. The FFT of the fiber input 
signal is multiplied by the transfer function and 
the inverse FFT of the result is calculated. To 
preserve a modular program structure, these fre- 
quency domain signal processing operations are 
performed by a SYSTID library model. Use of 
frequency domain signal processing requires the 
input signal to be periodic, and introduces an 
artificial delay equal to the size of the FFT. 
This approach offers the advantage of very fast 
execution time when compared to the time domain 
convolution approach. 

The two siodels which use time domain signal 
processing implement the discrete convolution sum 
using a tapped delay line algorithm. In these 
models the calculated transfer function is in- 
verse transformed to obtain the sampled impulse 
response of the single-mode fiber. Each sample 
of the impulse response is a tap gain in the 
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tapped delay line. It is not uncommon to have on 
the order of 2000 points to accurately character- 
ize the impulse response. Hence, this method is 
slower than the FFT approach but has the advan- 
tages of neither requiring periodic data nor 
introducing an artificial delay. 

PIN Diode Photodetector 

The PIN diode photodetector is a semicon- 
ductor device that produces an electrical current 
which is approximately proportional to the inci- 
dent optical power and is a potential source of 
degradation in a lightwave communication system 
link. In practical PIN receivers the effects of 
quantum or shot noise are overshadowed by the 
thermal noise of succeeding receiver stages. In 
this analysis, however, quantum noise has been 
added because of the need to eventually model an 
APD receiver, where the photodetector noise be- 
comes more important. 

The number of charge carriers (electron-hole 
pairs) in a PIN diode can be modeled mathemati- 
cally as a Poisson process, with a mean (and 
variance) proportional to the incident optical 
power _[2, 91. The average number of charge car- 
riers N generated in a time At is given by 


exponentially distributed. The mean of the in- 
terarrival time random variable is equal to the 
reciprocal of the mean of the Poisson process 
(for an observation time interval equal to uni- 
ty). For each sample in the simulation, the PIN 
diode model produces the Poisson distributed 
random number using the above relationship. If 
the average number of carriers, N, given in eq. 
(8) is large, then a Gaussian random number is 
generated with mean and variance equal to N . 
This is an extremely good approximation for TT 
larger than 50, and requires much less computer 
time. 

Bit Error Estimator Models 

There are currently two bit error rate (BER) 
estimation models implemented in SYSTID which are 
designed for lightwave communication system simu- 
lation. Both employ the semi-analytic approach 
to reduce the amount of simulation time required 
to obtain accurate estimates of link BER. The 
primary difference between the two is that one 
assumes thermal noise (additive white Gaussian) 
is added at a single point in the receiver, while 
the other model assumes two different thermal 
noise sources. 


N m P. (t) - 
hv inc 


At + n -At 


(8) 


where 


P, (t) 
inc 
n 

h 


incident optical power 

quantum efficiency 

Plank's constant 

optical frequency 

number of dark charge carriers/ 

second 


The probability that exactly k electron-hole 
pairs is generated in a time At can be written as 


P(n - k) 


(N)* (e)‘ N 
kl 


(9) 


which is simply the Poisson probability mass 
function. The PIN diode output current is pro- 
portional to the number of carriers, and is given 
by 


i(t) 



( 10 ) 


where 

- the charge of an electron 
n • the number of charge carriers genera- 
ted in At seconds 


Implementation of the mathematical model of 
the PIN detector shot noise in SYSTID is rela- 
tively straight-forward. The majority of the 
programming effort was involved in generating the 
Poisson random numbers , and providing the para- 
meters necessary to allow the user some flex- 
ibility in controlling the trade-off between 
simulation run-time and accuracy. 


It was necessary to modify a BER estimator 
model which already existed in SYSTID library so 
that it could be used with the on-off keying 
modulation scheme used in lightwave systems. 
These models automatically calculate the noise 
equivalent bandwidth of the portion of the system 
which lies downstream of the point where thermal 
noise enters the system. This is accomplished by 
injecting an impulse into this portion of the 
system (which is assumed linear) and summing the 
magnitude squared of the impulse response. The 
signal power at the point where noise enters the 
system i3 also measured and input to the BER 
estimator model. 

The BER estimator models measure the dis- 
tance from the received sample to the threshold, 
and calculate the probability of error for the 
values of signal power to noise power ratio re- 
quested by the user. This calculation is re- 
peated for the corresponding sample in each bit, 
and the results averaged. For example, if there 
are 16 samples per bit then there are 16 dif- 
ferent average error probabilities calculated. 
If sampling were assumed to occur at the optimum 
instant during each bit interval, the overall BER 
would be the minimum of the 16 (for example) 
different BER estimates. Because it is not al- 
ways reasonable to assume that the receiver tim- 
ing recovery is accurate enough to provide opti- 
mum sampling# the BER models discussed here pro- 
vide the user with tables of the probability of 
error estimate for each of the samples in the bit 
interval, maki"9 #4 * l * r evaluate the ef- 

fects of timing errors on system performance. 

IV . rrcuTVIXVE SYSTEM SIMULATION EXAMPLE 


Several methods for generating the Poisson 
random numbers are available, but the model we 
implemented utilizes the fact that in a Poisson 
process the interarrival time between events is 


To 

computer 

systems. 


help demonstrate the capabilities of 

ai.ul.tlon ot Uqhtw.v. communication 

a simple SYSTID simulation aseepie i. 


- 69 - 




Flgur* 1. Olgltitl lightwave common I cat Ion* system block dl 
of source output, fiber output, PIN diode output, equal tzetl 
filter output; equalization filter output eye diagram; and f 
power Is -kO dBm. 

presented which employs the lightwave models 
previously discussed (LED source, single-mode 
fiber, PIN photodiode, BER estimator)* The exam- 
ple is not Intended to show any new advances in 
lightwave system technology. The purpose of the 
example is to demonstrate the potential useful- 
ness of simulation In analyzing lightwave sys- 
tems* 

Figure 1 shows a digital lightwave communi- 
cations system block diagram with RIP postpro- 
cessor plots of the signals at various points* 
In Figure 1, the actual thermal noise samples are 
generated and shown for illustration. In the 
actual semi-analytic simulation, the effects of 
thermal noise are determined analytically - no 
thermal noise samples are generated, and the 
sample and decide block is replaced with a BER 
estimator model. 

The information source is modeled using an 
KR2 PN sequence generator, which in these ex- 
amples produces 120 bit PN sequences with 16 
samples per bit. The bit rate is 140 Hb/s. The 
FN generator selected for this example allows the 
user to specify the rise and fall times inde- 
pendently, and is therefore useful for modeling 
the fact that most optical sources have turn-on 
and turn-off characteristics which are not 
equal. The Gaussian approximation to the optical 
source spectral density curve is used in the 


igram with waveforms generated by SYST10. Shown are time ploti 
on filter Input (Including thermal noise), and equalization 
Iber transfer function magnitude (HHz). The received optical 


example* The source operates with a center wave- 
length of 1*28 micrometers and an RMS spectral 
width of 30 nanometers. 

The single-mode fiber in the example simula- 
tion uses the dispersion curve given in eq. (5), 
with a zero-dispersion wavelength equal to 1.31 
micrometers. A tabular fiber loss function was 
used with a fiber loss of approximately 0.5 dB /Km 
at 1.31 micrometers. The single-mode fiber model 
also includes optical source power level scaling, 
which in this example corresponds to a launched 
optical power of -28 dBm. Four thousand subin- 
tervals were used in the discrete implementation 
of the integral in eq. <1 ) to evaluate the fiber 
transfer function. The resulting transfer func- 
tion was calculated over the entire simulation 
bandwidth (-1120 to +1120 MHz) and contained 2240 
frequency samples. 

The remainder of the functional blocks in 
the example simulation model the various mecha- 
nisms of a lightwave receiver. The number of 
dark charge carriers generated by the PIN diode 
photodetector was set equal to 1.34 x 10 11 car- 
riers/second (21 nanoamperes) . The PIN diode 
photodetector represents the only non-linearity 
in the lightwave communications link. 

Thermal noise enters the receiver in the 
circuitry following the PIN diode as shown in 
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Figure- 1. When using the serai-analytic simula- 
tion approach in SYSTID, the ENBSIG model is 
inserted into the simulation at the point where 
the Gaussian noise enters the system. For the 
BER estimates to be correct, the portion of the 
system from the ENBSIG model to the end of the 
signal path must be linear. This turns out to be 
a very realistic assumption - this section of 
actual lightwave receivers is always linear. For 
purposes of illustration, the third order, 100 
MHz Butter worth filter is used to simulate the 
equalization filter in a lightwave system receiv- 
er* However, a variety of SYSTID library models 
are available to simulate the equalization fil- 
ter* 

To keep the example simple, a single thermal 
noise source was used. Inputs to the BER esti- 
mator model include the receiver filter output 
signal, the input bitstream, the decision thres- 
hold value, and the thermal noise power spectral 
density* For each value of received optical 

power, the BER estimator model produces a table 
which lists probability of error versus S/N for 
each of the 16 samples in a simulation bit. The 
table includes probability of error estimates for 
each sample in the bit interval so that the ef- 
fects of non-optimum sampling due to timing er- 
rors can be evaluated* 

A plot of probability of bit error versus 
received optical power is given for this example 
simulation in Figure 2. The plot in Figure 2 
assumes ideal sampling* 



V* CONCLUSIONS 

Techniques for modeling lightwave system 
functional blocks in SYSTID have been presented, 
along with examples of their usage in an end-to- 
end lightwave system simulation* Flexible-tcpo- 
logy simulation packages such as SYSTID can be 
used to analyze different types of communications 
links (i.e., satellite and lightwave). Future 
trends in simulation package design are leading 
toward combining this flexibility with more user- 
friendly graphics editing capabilities. 
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figure 2. Problblltiy of bit trror vi received 
averse optical power for the SYSTI0 example. 
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